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ABSTRACT 


\ 

■  This  work  focuses  on  a  method  for  two-dimensional  pattern  recognition. 
The  method  includes  a  global  search  scheme  for  candidate  windows  of  interest, 
based  on  Fourier  domain  cross-correlation.  A  method  to  normalize  the  input 
scene  by  local  rectangular  regions,  in  an  attempt  to  efficiently  approximate 
search  window  normalization,  is  presented.  Also  developed  is  a  candidate 
window  (potential  target)  similarity  measure,  based  on  the  normalized  Ll  and 
Euclidean  distances,  which  is  independent  of  the  template  DC  value  and  its 
energy.  Observations  on  the  performance  of  the  algorithm  applied  to  visual 
spectrum  photographs  of  tanks  in  a  realistic  environment  are  included.  Also 
included  is  the  software  needed  to  implement  the  algorithm  on  a  Data  General 
Eclipse  S/250  minicomputer. 


SCENE  ANALYSIS:  NON-LINEAR  SPATIAL 
FILTERING  FOR  AUTOMATIC  TARGET  DETECTION 

I.  INTRODUCTION 

GENERAL 

Target  detection  is  the  area  of  pattern  recognition 
concerned  with  locating  a  given  class  of  objects  embedded 
within  a  background  scene.  The  distinction  to  be  made  is 
whether  the  object  in  question  belongs  in  the  target  class  or 
in  the  non-target  class.  The  challenge  of  this  type  of  auto¬ 
matic  analysis  of  complex  visual  data  by  machine  has  proven 
to  be  surprisingly  difficult.  The  problem  remains  largely 
unsolved  despite  considerable  research  effort  [1:  28]. 

JUSTIFICATION 

Automated  systems  capable  of  identifying  objects  in  a 
cluttered  background,  irrespective  of  size,  orientation, 
illumination,  or  position  would  have  a  nearly  unlimited  range 
of  applications.  The  following  list  suggests  a  few  of  the 
major  areas  [2:  596]: 

1)  Document  processing:  recognition  of  unformatted,  non- 
segmented,  multi-font  characters; 

2)  Industrial  automation:  robot  assembly  and  inspection; 

3)  Military  applications:  analysis  of  reconnaissance  im¬ 
agery,  enhancement  of  weapon  delivery  display  systems,  and 
realization  of  the  autonomous  missile. 


BACKGROUND 


The  target  detection  process  can  be  broken  down  into 
three  steps  as  follows: 

1)  Determine  characteristic  features  of  the  target  tem¬ 
plate  that  are  invariant  to  size,  orientation,  energy,  and 
background  changes; 

2)  Perform  a  global  search  of  the  cluttered  input  scene 
for  the  features  of  step  (1)  and  identify  local  "windows"  of 
interest  (target  candidates) ;  and 

3)  Classify  the  information  within  the  candidate  windows 
into  one  of  the  two  disjoint  classes  of  either  targets  or 
non-targets. 

The  problems  associated  with  this  target  detection  pro¬ 
cess  are  numerous.  One  problem  is  in  analytically  describing 
the  quintessence  of  a  class  of  objects,  given  the  physical 
characteristics  of  only  a  finite  number  of  template  objects. 
Once  a  suitable  set  of  target  characteristics  (or  features) 
is  selected,  the  isolation  of  those  features  from  a 
background  scene  which  may  posess  some  of  the  same  features 
as  the  target  may  be  impossible.  For  these  and  other 
reasons,  no  general  solution  exists  to  detect  objects 
embedded  in  "real-life"  cluttered  scenes. 


Previous  studies  performed  in  this  area  at  the  Air  Force 
Institute  of  Technology  have  shown  promising,  although 
limited,  results.  One  of  these  studies  is  a  thesis  by 


"Picture  Correlation 


Israeli  Air  Force  Major  Moshe  Horev, 

Model  for  Automatic  Machine  Recognition"  [3],  In  it,  Horev 
describes  a  series  of  transformations  performed  in  the 
Fourier  domain  to  determine  the  size  and  orientation  of 
targets  within  a  cluttered  scene.  He  then  suggests  to 
perform  a  non-linear  (and,  hence  potentially  unpredictable) 
operation  of  combining  the  modified  phase  of  the  scene  with 
the  magnitude  of  the  template  in  hopes  of  enhancing  the 
target  objects;  this  procedure  is  known  as  the  "Phase  of 
the  Image,  Magnitude  of  the  Template"  (PIMT)  process.  An 
immediate  observation  is  that  the  success  of  the  process  may 
be  both  scene  and  template  dependent.  The  PIMT  process  has 
debatable  merit,  but  the  existence  of  a  scale-rotation 
transformation  does  suggest  an  area  for  further  research. 
With  the  premise  that  the  size  and  orientation  of  a  candidate 
target  within  a  scene  is  known  or  can  be  determined,  can  a 
process  be  developed  to  then  locate  and  accurately  classify 
the  target? 

PROBLEM/SCOPE 

This  study  includes: 

1) The  initial  development  of  a  target  detection 
algorithm  that  is  invariant  to  background  scene  composition 
or  energy; 

2)  The  implementation  on  a  digital  minicomputer  of  the 
software  routines  needed  to  process  the  input  imagery, 
perform  global  scene  searches  through  high-speed  correlation, 
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and  discriminate  local  candidate  windows  by  using  Ll  and  L2 
distance  metrics;  and 

3)The  preliminary  test  results.  A  complete  statistical 
performance  analysis  of  the  process  was  not  conducted  due  to 
time  constraints.  Suggestions  for  improving  the  performance 
of  the  process  are  included. 

ASSUMPTIONS 

The  following  assumptions  were  made  concerning  the  input 
test  scenes: 

1) The  size  and  orientation  of  targets  in  a  scene  is 
known,  or  can  be  determined  by  existing  methods,  three  of 
which  are  diffraction  pattern  sampling,  cross-correlation 
with  a  bank  of  scaled  and  rotated  templates,  or  by  performing 
a  scale-rotation  transformation; 

2)  The  digitized  scene  images  are  accurate  represen¬ 
tations  of  the  continuous  scenes  from  which  they  were 
obtained; 

3)  The  digitized  scene  images  may  have  been  corrupted  by 
additive  uncorrelated  noise;  and 

4)  Illumination  over  the  continuous  scene  is  varying 
slowly. 

OVERVIEW  OF  PRESENTATION 

The  next  chapter  discusses  methods  of  detecting  objects 
in  scenes  through  the  classic  technique  of  template  matching. 
Distance  factors  based  on  the  Ll  and  L2  distance  metrics  are 


derived;  these  factors  will  be  used  to  classify  candidate 
windows,  or  potential  targets. 


Chapter  three  briefly  discusses  the  discrete  Fourier 
transform  and  some  of  its  properties.  A  method  for 
performing  high-speed  correlation  by  multiplication  in  the 
Fourier  domain  is  given;  cross-correlation  will  be  used  in 
the  detection  process  to  perform  searches  for  potential 
targets  in  the  input  scene. 

In  chapter  four  a  scene  normalization  scheme  is 
presented.  The  normalization  is  necessary  to  improve  the 
ability  of  the  cross-correlation  to  locate  candidate 
targets. 

The  software  needed  to  implement  the  detection  process 
is  described  in  chapter  five.  The  source  code  has  been 
included  in  the  appendix. 

In  chapter  six,  some  observations  of  detection  process 


are  made.  Explanations  are  given  for  the  weaknesses  of  the 
algorithm,  and  suggestions  for  improving  the  performance  are 
discussed. 


a 


II.  TEMPLATE  MATCHING 


Often  in  scene  analysis  problems  a  simple  question  is  to 
be  answered:  Does  the  input  scene  contain  a  previously 
specified  object?  A  technique  classically  employed  to 
determine  the  presence  of  an  object  is  the  fundamental  method 
of  template  matching,  in  which  the  template  brightness 
function  is  compared  point-by-point  with  the  scene  brightness 
function.  In  most  cases,  a  perfect  template  match  will  not 
be  found,  so  some  realistic  distance  measure  D(m,n)  indi¬ 
cating  the  degree  of  similarity  between  the  template  window 
and  the  scene  needs  to  be  computed  for  all  possible  points  in 
the  scene. 

LI  and  L2  DISTANCE  MEASURES 

Let  the  array  (or  vector)  t(x,y)  represent  in  some  sense 
the  template  pattern,  and  let  the  array  s(x,y)  represent  the 
scene  to  be  searched.  For  our  purposes  of  discussion,  it  is 
immaterial  how  the  arrays  are  obtained,  whether  from  digi¬ 
tizing  (sampling  and  quantizing)  the  continuous  brightness 
distributions,  infrared  distributions,  or  some  function  of 
these  distributions  (for  example  only  the  low-frequency 
Fourier  components) .  Two  common  definitions  of  distance 
measures  used  are  given  by  equations  (1)  and  (3)  [4:  279]. 

Dl(m,n)  *  ^  £  ls(Xry)  “  t(x-m,y-n)|  (1) 
x  *  y 


s  £  2^V[s(x,y)  -  t  (x-m,y-n)  ] 
x  *  y 
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D2(m, n) 


(3) 


D2(m,n)  2  Is (xry)  -  t(x-m,y-n)J 

f  x  *  y 

*  —  for  all  x,y  such  that  (x-m,y-n)  is  within  the  area  of 
overlap  of  the  scene  and  template  windows  (0  <  m  <  M+Jr 
0  <  n  <  N+K  for  a  JxK  template  and  an  MxN  search  area) .  See 
Figure  1  for  an  illustration  of  the  labelling  convention 
used.  Note  that  m  and  n  represent  a  specific  translation 
between  s(x,y)  and  t(x,y). 
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The  definitions  for  Dl(m,n)  and  D2(m,n)  are  known  as  the 
metrics  based  on  the  LI  and  L2  norms  respectively.  The 
D2(mfn)  measure  is  also  known  as  the  standard  Euclidean 
distance  between  two  vectors,  that  is 


A  1  A 

D2(m,n)  =  (s  -t)  (s  -  t) 


(4) 


where  t  = 


t (xl-m,yl-n) 


t (x2-m,y2-n) 


t (xn-m,yn-n) 

m 

NORMALIZED  CROSS-CORRELATION 

Insight  can  be  gained  from  equation  (3)  by  expanding  it 
as  follows: 

2  "  fs2 

x  *  y 

2  -1 

(5 


D2(m 


Equivalently, 


,n)  =  EE  [s  (x,y)  -  2s(x,y)t (x-m,y-n) 
x  *  v  L 

-m,y-n)J 


2 

+  t  (x- 


2  2 

D2(m,n)  *  EE  s  (x,y)  -  2  EE  s ( x , y ) t ( x-m, y~n ) 

x  *  y  x  y 


+  (x-m,y-n) 

x  y 


(6 


Let  Es,  Et,  and  Rst  be  defined  as  follows: 


Es(m,n)  = 

x  *  y 


(7 


Et(m,n)  =  EE  t  (x-m,y-n) 


(8 


(9) 


Rst(m,n)  =  y*!  s  (x,  v)  t  tx-m.v-n) 

x  y 

(There  is  no  restriction  on  the  range  of  x  and  y  in  equations 
(8)  and  (9)  as  the  template  function  is  considered  to  be  zero 
outside  the  area  of  interest.) 

The  term  Es(m,n)  represents  the  scene  energy  (or 
equivalently  the  vector  length)  of  the  search  window ,  which 
will  vary  ever  the  search  area.  The  term  Et(m,n)  represents 
the  template  energy,  which  is  constant  for  all  values  of  m 
and  n.  The  term  Rst(m,n)  is  the  cross-correlation  between 
the  scene  and  t  template,  and  generally  is  largest  when 
the  distance  D2(m,n)  i '  smallest.  The  cross-correlation  term 
is  not  an  absolute  measure  of  the  template  difference, 
however,  since  the  scene  window  energy  Es(m,n)  is  position 
variant.  For  this  reason,  Rst(m,n)  is  normalized  to  achieve 
invariance  to  input  energy  [5:  553].  The  normalized  cross¬ 
correlation,  denoted  by  Nst(m,n),  is  defined  as  follows: 

Rst (m,n) 

Nst  (m,n)  =  -  (10) 

y  Et(m,n)  y Es (m,n) 

where  the  usual  restriction  is  placed  on  x  and  y  for  the 
computation  of  the  scene  energy  term. 

The  significance  of  this  normalization  can  be  realized 
by  appealing  to  the  Schwarz  inequality  [6:  159].  The  Schwarz 


inequality  is  stated  as  follows: 


For  two  real  functions  f(x)  and  g(x)  defined  on  a<x<b, 

£f(x)g(x)  <  VElf(x)!2  ^  5^[g  (x)  1  (11) 

with  equality  when  f(x)=kg(x),  where  k  is  a  constant  scale 
factor.  As  applied  to  the  cross-correlation  terms. 


E  Es<x  ,y) t (x-m,y-n)  < 
x  y  - 

ViE  Es2(*'y>J 

*  x  *  y 

or,  by  rearranging  terms 


[E  E  1 2(x-m,y-n)  ] 

x  y 


(12) 


Nst(m,n)  <  1  (13) 

with  equality  only  when  the  scene  area  under  consideration 
exactly  matches  the  template.  Thus  the  normalized  cross¬ 
correlation  can  be  used  as  a  decision  criterion  regardless 
of  the  distribution  of  the  non-normal ized  scene  energy, 
assuming  that  the  scene  energy  can  be  recomputed  for  each 
shift  of  the  search  window.  A  simple  decision  rule  would 
classify  the  scene  window  information  into  the  target  class 
only  when  Nst  exceeded  some  preset  upper  threshold  value, 
into  the  non-target  class  when  Nst  was  below  a  lower 
threshold,  and  would  not  make  a  decision  when  Nst  was  between 
the  thresholds. 


With  the  constraint  that  the  scene  and  template  can  take 
on  only  non-negative  values,  a  lower  bound  for  the  threshold 
value  is  zero.  Fortunately,  a  tighter  bound  closer  to  unity 


can  be  determined  for  cross-correlations  with  scenes  of 
interest  (ones  that  approach  the  template  in  form).  Consider 
a  scene  of  constant  value  c,  c>0.  Then  the  cross-correlation 
value  may  be  determined  by  applying  equations  (7)  through 
(10)  : 

Es  =  Z  L  c  (14) 

x  y 


2 

Es  =  c  JK 


(15) 


[Nst | s=c] 


<=£  £  t (x-m,y-n) 
x  y 


(16) 


£  £  t (x-m,y-n) 
x  y 

Nst  |  c  =  - -  — - 

V  EtJK 


(17) 


with  J  =  width  of  search  window 
K  =  length  of  search  window 


Note  that  the  normalized  cross-correlation  between  the  tem¬ 
plate  and  a  scene  of  constant  non-zero  value,  Nst|c,  is 
independent  of  the  scene  value.  Scene  windows  which  yield  a 
Nst  less  than  Nst|c  need  not  be  considered  for  further 
discrimination. 

NORMALIZED  LI  and  L2  DISTANCES 

Two  other  distance  measures  which  take  into  account  the 
array  energies  are  the  normalized  LI  and  L2  distance 
measures,  defined  in  Eqs.  (18)  and  (19). 
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NL1 (m, n) 


x  *  y 


s (x,y) 
SL1E 


t (x-m.y-n) 
TL1E 


NL2(m,n)  =  ' 


x  *  y 


(s (x,y)  t (x-m,y-n)\: 

Es(m,n)  Et  / 


with 


slie  =53 

x  *  y 

TL1E  =  53  53  t(x-m.y-n) 
x  y 


The  normalized  L2  distance  can  be  determined  more  effi¬ 
ciently  by  Eq.  (22) : 


NL2(m,n)  =i/2  [1  -  Nst(m,n)] 


One  way  of  "visualizing"  the  normalized  L2  distance  is 
to  think  of  it  as  the  Euclidean  distance  between  the  points 
where  the  scene  vector  and  the  template  vector  intersect  the 
unit  hypersphere.  Thus  NL2  is  dependent  only  upon  the  angle 
between  the  vectors,  and  not  on  the  vector  lengths. 


The  maximum  normalized  distances  to  be  considered  as 
possibly  identifying  a  target  location  will  be  those 
corresponding  to  the  distances  computed  between  a  template 
and  a  scene  with  a  constant  value.  These  distances  are  given 
in  Table  I.  for  typical  tank  template. 


TABLE  I.  DISTANCES  FROH  TEMPLATE  H3  TO  A  CONSTANT  VALUED  SCENE 
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DISTANCE  FACTORS 

At  this  point  the  distance  factors  used  to  classify 
candidate  target  windows  will  be  introduced.  The  corre¬ 
lation,  LI,  and  L2  factors  correspond  to  the  normalized  cor¬ 
relation,  LI,  and  L2  distances  linearly  scaled  into  a  0-100 
range,  with  100  corresponding  to  a  exact  match  and  0 
corresponding  to  the  distance  to  a  constant  valued  scene. 
Consider  the  mappings 


Nst  — >  0 
NL2  0 
NL1  ->  0 

and 

Nst  *1—5 
NL2  *0—5 


for  Nst 
for  NL2 
for  NL1 

100 

100 


<  Nst  |  c 

>  NL2  |  c  (23) 

>  NL1  |  c 

(24) 


NL1  *  0  — >  100 


The  functions  to  achieve  these  mappings  are 


CFACTOR  = 


L2FACTOR  = 


Ll FACTOR  = 


{100 [ (NSt  -  NSt I C ) /  (1  -  Nst | C) ] 
0 

(100 [1  -  NL2/NL2 | c] 

0 

!100[1  -  NLl/NLl | c] 

0 


Nst  >  Nst 
else 


c 

(25) 


NL2  <  NL2  |  c 
else 


(26) 


NLl  <  NLl  |  C 
else 


(27) 


A  score  is  computed  to  take  into  account  all  three 
distance  measures  as  follows: 

SCORE  =  y  (CFACTOR  *  L2FACTOR  *  LlFACTOR)  (28) 

Note  that  SCORE  intentionally  favors  the  Euclidean 
measure  over  the  Ll  metric  (recall  that  the  normalized 
correlation  is  an  invertible  function  of  the  normalized 
Euclidean  distance) .  The  SCORE  will  always  be  a  number  from 
0  to  100  inclusive.  The  behavior  of  the  SCORE  of  a  template 
measured  against  itself  for  various  window  shifts  is  given  in 
Table  II. 


One  of  the  problems  in  determining  the  distance  func¬ 
tions  is  that  they  are  computationally  expensive,  often 
infeasible,  because  of  the  large  size  arrays  required  for 
most  applications.  This  is  the  case  with  many  linear 


TABLE  II.  BEHAVIOR  OF  DISTANCE  FACTORS 
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processing  algorithms.  Indirect  computational  technique 
based  on  unitary  transforms  permit  more  efficient  linea 
processing  than  conventional  methods.  Tht  next  ^apte 
introduces  the  Fourier  transform  as  a  linear  processing  tool 


III.  LINEAR  PROCESSING 


An  efficient  method  of  linear  processing  is  through  the 
use  of  unitary  transforms.  A  unitary  transform  meets  the 
following  three  conditions  [5:  232]: 

1) It  is  a  linear  transformation; 

2)  Its  operation  is  exactly  invertible;  and 

3)  Its  operating  kernel  satisfies  certain  orthogonality 

conditions. 

A  unitary  transform  of  particular  importance  in  the  field  of 
image  processing  is  the  two-dimensional  Fourier  transform. 
In  addition  to  its  use  as  a  linear  processing  tool,  the 
Fourier  transform  provides  a  means  of  extracting  features 
from  images.  For  instance,  the  center  or  DC  term  is  propor¬ 
tional  to  the  average  image  brightness.  The  low-frequency 
terms  contain  the  gross  form  information,  while  the 
high-frequency  terms  indicate  the  amplitude  and  orientation 
of  the  edges  (the  detail)  . 


TWO-DIMENSIONAL  DISCRETE  FOURIER  TRANSFORM 

Consider  a  two-dimensional  periodic  sequence 

?(x,y)  =  g (x+qM,y+rN)  (29) 

where  q  and  r  are  integers,  and  M  and  N  are  the  periods  in 
the  x  and  y  direction.  Such  a  sequence  can  be  represented  by 
a  finite  sum  of  exponentials  in  the  form 


M-l  N-l 


G(fx,fy)exp[ j2  Tf  (xf x/M+yfy/N)  ]  (30) 


fx»0  fy=0 
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1 


where 


G(fx,fy)  2 ?(x,y)exp[-j27<  (xf x/M+yfy/N)  ] 


(31) 


x  y 


and  j=  v “I  . 


i 


l 


Note  that  G(fx,fy)  will  have  the  same  periodicity  as  the 
sequence  g(x,y) .  If  a  finite  area  sequence  g(x,y)  is  con¬ 
sidered  to  be  one  period  of  g(x,y),  and  G(fxrfy)  is  taken  to 
be  one  period  of  G(fx,fy),  then  g(x,y)  and  G(fx,fy)  will  form 
a  discrete  Fourier  transform  pair.  In  equation  form  [7: 
117]  , 


9(x,y)  = 


G(fx,fy)  = 


—  ^^G(fx,fy)exp[  j27t(xf  x/M+yfy/N)  ]  0<x^M-l 

MN  Tx  fy  (Ky£N-l 

(32) 

0  otherwise 


'S*'  g  (x,y)  exp  [-j27T(xf  x/M+yfy/N)  ]  O^fx^M-1 

“  O^fy^N-1 

(33) 

0  otherwise 


The  notation  to  be  used  to  indicate  a  Fourier  Transform  pair 
is 

g(x,y)  -« - ►  G (fx, fy)  (34) 

Equivalently, 

F{g(x,y)}  =  G (f xr  fy)  (35) 

and 


-1 

F  {G ( f x, fy ) }  *  g (x,y) 


(36) 


17 


Having  defined  the  Fourier  transform,  two  theorems  for 
use  in  later  developments  will  be  stated  without  proof  [7: 
110]  . 

SHIFT  THEOREM 

If  t(x,y)« - ►  T(fx,fy)  ,  (37) 

then  t(x-m  ,y-n  )  * - ►  exp[-j2  Tf  (mf x/M+nfy/N)  ]T(fx,fy) 

REVERSAL  THEOREM 

If  t  (x,y)  < - ►  T(fx,  fy)  , 

then  t(-x,-y)  * - ►T*(fx,fy), 

where  T*  is  the  complex  conjugate  of  T.  (38) 


The  convolution  theorem  suggests  a  method  for  performing 
correlation  in  the  Fourier  domain.  It  will  now  be  stated 
with  its  proof,  modeled  after  a  proof  for  continuous  signals 
[4s  307]  . 

CONVOLUTION  THEOREM 

For  s(x,y)  + - ►Stfx^y)  and  t(x,y)  <« - ►  T(fx,fy)  , 

F{s(x,y) *t (x,y) }  =  S(fx,fy)T(fx,fy)  (39) 

Proof : 

By  definition, 

F(J]Zs(x»y)t(m-x»n-yi}  =  (40) 

x  y 

££[  Z)  £  s  (x,y)  t  (m-x,  n-y)  J  exp  [- j2  71  (mf  x/M+nfy/N)  ] 
m  n  x  y 

Interchanging  the  summation  order, 

■  T,  T,  s(x*y)  {  £  L  t (m-x,n-y)exp[-j27t  (mf x/M+nfy/N)  ]  } 
x  y  m  n 
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scenes,  so  that  T*(fx,fy)  needs  to  be  computed  just  once  and 
stored. 

ajrray  extension  for  linear  correlation 

Care  must  be  taken  in  choosing  sequence  and  transform 
lengths.  Consider  the  linear  correlation  between  two  N-point 
sequences  R(m)  =  ^s(x)t(x-m),  where  R(m)  will  have  up  to 

2N-1  non-zero  points.  The  indirect  correlation  method  using 
N-point  discrete  Fourier  transforms  will  result  in  an  N-point 
sequence,  which  is  the  circular  correlation  of  the  input 
sequences.  To  obtain  the  linear  correlation,  the  discrete 
Fourier  transforms  must  be  computed  on  the  basis  of  2N-1  or 
more  points,  with  the  input  sequences  extended  with  at  least 
N-l  zeros.  In  general,  for  s(x)  of  length  Si  and  t(x)  of 
length  Tl,  the  indirect  linear  correlation  method  requires 
that  discrete  Fourier  transform  be  computed  on  the  basis  of 
at  least  Sl+Tl-1  points. 

For  the  two-dimensional  case,  the  sequence  arrays  are 
extended  as  follows  (See  Figure  3)  [5s  288] : 

1)  Imbed  the  TlxT2  template  image  sequence  in  the  lower 
right  quadrant  of  an  all  zero  MlxM2  matrix; 

2)  Imbed  the  SlxS2  input  scene  in  the  upper  left  quadrant 
of  an  all  zero  MlxM2  matrix; 

3) Compute  all  discrete  Fourier  transforms  on  the  basis 
of  M1>S1+T1-1  and  M2>S2+T2-1  to  avoid  wrap-around  error. 


By  taking  advantage  of  the  speed  gained  by  implementing 
the  indirect  correlation  method,  the  ability  to  normalize  the 
scene  search  window  during  the  correlation  process  is  lost. 
This  inability  is  a  major  shortcoming  of  the  indirect 
correlation  method.  An  alternate  normalization  scheme  must 
be  implemented  to  approximate  the  window-by-window 
normalization  method. 


IV.  NORMALIZATION  SCHEME 


Consider  dividing  the  scene  to  be  normalized  into  a  grid 
of  rectangles.  A  compromise  between  global  normalization  and 
search  window  normalization  would  be  to  divide  each  of  the 


scene  values  within  a  given  rectangle  by  some  constant  that 
is  proportional  to  the  square  root  of  the  energy  of  that 
rectangle.  This  normalization  scheme  may  be  implemented  as 
in  Figure  4. 


(H  .V  )  <h  >  v  > 


(H  .V  > 


(H  .  V  ) 
1  2 


(H  .V  > 
i  J 


(H  .V  ) 
i  k 


(H  .V  ) 

1  k 


FIGURE  4:  NORMALIZATION  GRID  SETUP 


(Hi,Vj)  —  grid  rectangle  in  grid  row  j,  grid  column  i 
N  —  of  pixel  rows 

M  —  of  pixel  columns 


(i, j,M,N,N/k,  and  M/1  are  all  integers). 


Define 


Eij  =  ^2^2  s  (x»y) 

x  y 

(i-1) M/1  +  1  <  x  <  iM/1 
(j-l)N/k  +  1  <  y  *  jN/k 


Then  for 


s' (x,y)  —  normalized  array  value 

Nij  —  normalization  coefficient 

Eij  —  energy  in  grid  rectangle  (Hi,Vj) 


where 


s' (x,y)  =  s(x,y)/Nij 


Nij  =  a  y  Eij 


Some  thought  must  be  given  to  the  size  of  the  rectangles 
chosen,  for  as  the  size  is  increased,  the  clutter  energy-to- 
target  energy  ratio  is  also  increased  (ideally  this  ratio 
should  be  zero) .  As  the  size  of  the  rectangles  is  decreased, 
the  scene  begins  to  lose  contrast  as  the  normalized  values 
asymptotically  approach  a  constant  value  (namely  1  when  the 
rectangle  size  is  lxl)  .  Another  problem  accompanying  the 
decrease  in  rectangle  size  is  the  possibility  of  sectoring 
part  of  a  target  into  separate  rectangles  increases;  this 
could  have  a  deleterious  effect  on  the  cross-correlation 
function.  See  Figure  5  for  an  illustration  of  a  normalized 
scene. 


c 
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Consider  the  following  two  cases,  shown  in  Figure  6,  in 
which  a  3x€  normalization  grid  has  been  chosen. 


•7  ■  —  LCr-CLPCY 

[__  ,  —  ;iip:h30cy 


FIGURE  6i  EXAMPLE  OF  TARGET  SECTORING  PROBLEM. 

In  case  I,  the  target  will  be  satisfactorily  normalized 
against  the  high  energy  background,  allowing  a  successful 
global  search  by  correlation.  In  case  II,  the 
clutter-to-target  area  is  large  for  all  grid  rectangles. 
The  high  energy  areas  dominate,  yielding  false  peaks  in  the 
correlation  function. 

In  the  next  chapter,  the  software  necessary  to  process 
the  images  and  test  the  target  detetection  process  is 
described.  Also,  the  detection  process  is  further  discussed. 


V.  SOFTWARE  DESCRIPTION 
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The  software  that  was  used  in  this  study  is  described  in 
this  chapter.  The  programs  (indicated  by  capitalized  titles) 
have  been  grouped  into  the  following  categories: 

1)  Image  input  and  output; 

2)  Scene  and  template  synthesis; 

3)  Correlation  implementation; 

4)  Process  evaluation;  and 

5)  Support  subroutines. 

All  source  code  that  was  generated  during  this  thesis 
effort  and  added  to  the  AFIT  Signal  Processing  Laboratory 
software  archives  is  included  in  the  appendix.  All  programs 
are  written  in  Data  General  (DG)  Fortran  5  (except  for  VIDE07 
and  NMOVE,  which  are  written  in  Fortran  IV) .  All  programs 
were  written  by  James  Cromer,  with  the  exception  of  PLTTRNS, 
INVERSE,  and  DIRECT  (by  Ronald  Schafer)  ,  and  IOF,  UNPACK,  and 
REPACK  (by  Robin  Simmons) . 


IMAGE  I/O 


Before  a  digital  computer  can  be  used  to  analyze  an 
image,  the  image  must  first  be  converted  to  a  form  usable  by 
the  computer.  Specifically,  the  image  must  be  represented  in 
some  sense  by  an  array  of  numbers.  The  process  used  to 
obtain  this  array  is  known  as  digitization,  in  which  some 
image  parameter  is  sampled  and  quantized  at  points  throughout 
the  scene.  In  this  study,  the  achromatic  brightness,  or  gray 
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level ,  is  sampled  in  a  256  column  by  256  row  quadruled  grid 
format,  and  quantized  into  one  of  16  levels  (4-bit 
digitization) .  The  resulting  array  values  are  referred  to  as 
"pixels",  short  for  "picture  elements." 

The  equipment  used  in  the  process  include  a  standard 
video  monitor,  a  Cohu  6150  vidicon  camera  with  6950  camera 
controls,  and  3  Tecmar  digitizer  boards  (A/D  converter, 
direct  memory  access,  D/A) .  The  Tecmar  digitizer  is 
interfaced  with  a  DG  NOVA  2  processor  via  a  CROMEMCO  Z-80 
based  microcomputer.  The  NOVA  terminal  can  be  used  to 
communicate  with  the  A/D/A  converter  with  the  Fortran 
callable  subroutine  CHANNEL,  developed  in  an  earlier  AFIT 
thesis  [8] .  High  speed  processing  can  be  performed  on  the 
digitized  images  with  the  powerful  DG  ECLIPSE  S/250 
minicomputer.  Both  DG  machines  are  16-bit  processors.  See 
Figure  7  for  a  schematic  of  the  equipment  layout. 

The  program  which  controls  image  input  (digitization) 
and  output  (display)  through  CHANNEL  is  VIDEO? .  When  running 
VIDE07  and  the  input  option  is  chosen,  seven  digitized 
versions  of  an  input  image  are  stored  in  files  named  "AO" 
through  "A6",  which  can  be  averaged  together  later.  When  the 
output  option  is  chosen,  the  user  is  given  the  choice  to 
display  from  one  to  ten  files  named  "AO"  through  "A(n-l)", 
where  n  is  the  number  of  files  to  be  displayed.  In  addition 
to  the  main  purpose  of  checking  images  before  averaging,  n 
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Figure  7:  Equipment  Layout  of  the  AFIT  Signal  Processing  Laboratory 
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images  of  interest  may  be  easily  presented  in  sequence  for 
demonstrations  by  re-naming  the  images  "AO"  through  "A(n-l) ." 
The  third  mode  of  VIDE07  allows  the  user  to  display  an 
existing  file  any  number  of  times  consecutively.  This  mode 
is  used  when  the  D/A  converter  malfunctions,  and  usually 
several  attempts  to  view  an  image  must  be  made  before  a 
satisfactory  image  is  displayed. 

Files  created  by  VIDE07  are  written  to  disk  in  what  is 
referred  to  as  "packed  video  form."  The  file  is  "packed" 
because  4  pixels  are  stored  in  each  16-bit  word.  Packed  form 
is  ideal  for  minimizing  storage  requirements,  while  posing 
only  minor  processing  inconveniences.  In  packed  form,  video 
files  will  be  64  blocks  long,  where  one  block  is  256  16-bit 
words.  Thus  one  image  requires  only  32K  bytes  of  memory. 
Note  that  each  block  holds  4  packed  video  rows.  As  a  result, 
the  processing  programs  operate  on  multiples  of  4  rows  at  a 
time  between  RDBLK  calls.  The  fastest  way  to  transfer  data 
from  disk  storage  to  core  memory  is  by  the  RDBLK  call,  which 
passes  data  in  the  data  channel  mode.  The  data  channel  mode 
of  moving  data  does  not  require  program  control  once  a 
transfer  is  initiated.  Figure  8  shows  more  clearly  the 
relation  between  packed  and  unpacked  forms. 

The  seven  digitized  images  created  by  VIDE07  can  be 
averaged  to  produce  an  output  image  that  has  an  improved 
signal  to  noise  ratio.  The  program  that  does  this  averaging 
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is  called  QUICKAVE7,  which  creates  an  output  file  named 
"AVERAGE7 .VD",  where  the  .VD  extension  indicates  a  video 
file.  The  seven  files  to  be  averaged  are  assumed  to  be  named 
"AO"  through  "A6." 

It  is  often  desirable  to  produce  a  hard  copy  of  a 
digitized  image.  One  of  the  ways  this  can  be  done  is  by 
displaying  an  image  through  VIDE07,  then  activating  the 
Tektronix  4632  Video  Hard  Copy  Unit,  which  makes  a  photocopy 
of  the  image  being  sent  to  the  video  monitor.  This  method  is 
acceptable  most  of  the  time,  but  for  numerous  reasons  it  is 
necessary  to  produce  hard  copies  of  stored  images  with  the 
Printronix  P-300  lineprinter.  The  program  originally  written 
to  do  this  is  DISPLAY,  by  Robin  Simmons  in  an  earlier  thesis 
[9],  DISPLAY  used  3x3  dot  patterns  to  simulate  the  16 
gray-levels,  which  resulted  in  two  shortcomings.  Distortion 
occurs  in  the  picture  because  the  P-300  horizontal  dot 
density  is  less  than  the  vertical  density,  resulting  in  a 
1.2:1  aspect  ratio.  Also  the  3x3  patterns  do  not  fully  take 
advantage  of  the  16  gray-levels  available.  DISPLAY  was 
modified  to  solve  these  problems  by  using  a  combination  of 
3x3  and  3x4  dot  patterns.  Up  to  four  different  dot  patterns 
are  used  per  gray-level,  instead  of  just  one,  and  the  aspect 
ratio  is  very  nearly  1:1.  Other  modifications  include 
allowing  the  user  to  choose  the  number  of  rows  to  be 
displayed  along  with  the  starting  row.  The  run  time  for  an 
11x13  image  hard  copy  was  reduced  to  less  than  90  seconds. 
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down  from  5-6  minutes.  The  modified  program  is  called 
PICTURE. 

SCENE  and  TEMPLATE  SYNTHESIS 

After  an  image  is  digitized  and  stored ,  the  next  step  is 
to  create  a  template  or  scene  to  be  used  in  the  correlation 
process.  To  do  this  programs  were  written  to  improve  the 
image,  create  a  scene  or  template  by  combining  images,  and  to 
put  the  images  into  the  correct  format  for  the  correlation 
process. 

The  program  REMOVE  was  written  to  perform  a  3x3  pixel 
mask  processing  of  an  image  for  the  purpose  of  noise  re¬ 
moval.  The  main  program  handles  the  bookkeeping  of  passing 
the  three  rows  to  be  operated  on  to  the  subroutine  TEST3 , 
which  produces  the  noise-removed  output  row.  The  subroutine 
UNPACK 2  is  used  to  unpack  the  video  rows  from  four  pixels  per 
word  to  one  pixel  per  word.  REMOVE  was  not  used  extensively, 
but  is  included  to  demonstrate  an  efficient  method  to 
perform  mask  processing.  The  mask  function  can  be  changed  by 
modifying  TEST3.  TEST3  presently  computes  the  difference 
between  the  center  pixel  value  and  the  average  value  of  the 
surrounding  pixels.  If  the  difference  is  greater  than  some 
threshold,  the  center  pixel  value  is  modified  accordingly. 

If  REMOVE  and  QUIKAVE7  fail  to  produce  a  satisfactory 
image,  a  histogram  can  be  generated,  and  then  modified  to 
enhance  the  image.  The  program  to  produce  the  histogram  is 


called  EVIDHIST  (the  "E"  Indicates  an  Eclipse  only  program) . 
TOMER  modifies  the  histogram  by  a  mapping  function  of  the 
type:  0  a,  1  b,  .  .  .  ,  15  p,  where  the  user  defines 
the  new  values  of  "a"  through  "p".  TONER  is  used  to  increase 
the  contrast  or  raise  the  average  brightness  level  when  de¬ 
ficiencies  occur  due  to  A/D  or  camera  gain  misadjustment. 

To  create  a  template  from  a  scene  with  a  target  in  it,  a 
"window"  is  placed  over  the  target  information,  and  the 
background  is  set  to  some  constant  value  (usually  0  or  15) . 

Program  NMOVE  allows  the  user  to  specify  a  template 
scene  file,  a  background  file,  and  a  combined  filename.  The 
template  window  size  and  position  are  variable,  as  is  the 
combined  window  position.  Figure  9  demonstrates  the 
capability  of  NMOVE. 
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FIGURE  9:  LEFT:  MOVE  INPUT  SCENE.  RIGHT:  SCENE  CREATED  BY  NHOVE. 
THE  BACKGROUND  SCENE  IS  THE  U  GRAY-LEVEL  BARPATTERN. 


The  next  step  requires  that  the  "negative"  image  be 
formed,  using  the  equation  NEGATIVE  =  15  -  POSITIVE.  This 
forces  the  expected  high  energy  background  to  become  a  low 
energy  background,  improving  the  correlation  results.  TONER 
can  also  be  used  for  this  purpose,  but  NEGATE  is  used  in 
macro  files,  as  it  requires  no  user  input. 

The  last  step  before  the  correlation  sequence  is  to  do  a 
four-to-one  reduction  to  imbed  the  256x256  scene  into  the 
upper  left  quadrant  (lower  right  for  template)  so  that  linear 
correlation  will  be  obtained.  Program  REDUCE  does  this  by 
averaging  four  pixels  to  create  one  output  value.  This  re¬ 
duction  is  also  effectively  a  low-pass  filtering  operation. 
See  Figures  10  and  11  for  flowcharts  of  the  scene  and 
template  synthesis  processes. 

CORRELATION  IMPLEMENTATION 

The  correlation  method  used  is  to  perform  a  multi¬ 
plication  in  the  frequency  domain,  then  inverse  transform  the 
product,  resulting  in  a  correlation  in  the  spatial  domain. 
The  complex  arrays  multiplied  bave  been  obtained  from 
Fourier  transforming  normalized  image  arrays. 

The  normalization  scheme  used  is  the  rectangle  grid 
normalization  described  in  Chapter  4.  The  program  which 
carries  out  the  normalization  is  suitably  named  NORMALIZE. 
The  number  of  grid  rows  and  grid  columns  are  chosen  by  the 
user.  Only  the  upper  left  or  lower  right  quadrants  are 
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FIGURE  11:  FLOWCHART  OF  TEHPLATE  SYNTHESIS. 
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normalized;  the  remainder  of  the  array  is  set  to  zero.  The 
input  file  is  a  256x256  packed  video,  and  the  output  is  a 
256x256  complex  file. 

After  the  normalized  scene  and  template  complex  files 
have  been  created,  they  are  Fourier  transformed  by  the 
program  DIRECT.  The  resultant  arrays  are  then  complex 
multiplied  by  CMULTIPLY,  and  the  product  is  inverse 
transformed  by  INVERSE.  INVERSE  and  DIRECT  are  Eklundh  FFT 
algorithm  based  programs  written  by  Ron  Schafer. 

The  file  created  by  INVERSE  is  a  complex  array  with  the 
imaginary  part  zero  (since  the  scene  and  template  arrays  are 
np  always  real) ,  and  with  the  real  part  having  values  between 

zero  and  two.  CTOI  conserves  file  space  by  converting  the 
complex  array  to  integer  form  by  multiplying  each  real  number 
by  16384  to  take  advantage  of  the  16  bit  word.  The  integer 
file  uses  only  one-fourth  as  much  disk  memory  as  a  complex 
file. 


The  program  CTOV  can  be  used  to  convert  a  complex  file 
(imaginary  part  assumed  zero)  into  a  video  file.  It  performs 
a  linear  scaling  of  the  input  file  into  a  0-15  range.  CTOV 
was  used  (along  with  PICTURE)  to  display  the  complex 
normalized  scene  of  Figure  5. 
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The  last  step  in  the  correlation  process  is  to  combine 
the  results  of  several  correlations  between  a  scene  and  a  set 
of  templates.  Program  IMULTIPLY  computes  the  geometric  mean 
of  two  correlation  functions  by  determining  the  square  root 
of  the  product  of  the  input  arrays.  An  example  of  when 
IMULTIPLY  may  be  used  is  when  the  correlation  functions 
created  from  the  left-half  and  the  right-half  of  a  template 
are  to  be  combined.  IMULTIPLY  can  also  be  used  to  combine 
the  positive  correlation  (positive  scene  with  positive 
template)  with  the  negative  correlation  (negative  scene  with 
negative  template) .  See  Figure  12  for  a  flowchart  of  the 
correlation  implementation. 

PROCESS  EVALUATION 

The  function  resulting  from  the  correlation  process  next 
needs  to  be  evaluated.  The  correlation  function  can  be 
evaluated  by  viewing  it,  or  by  a  numerical  analysis  of  its 
peaks. 

PLTTRNS,  by  Schafer,  enables  the  user  to  view  3-D, 
contour,  and  row  plots  on  the  Tektronix  4010  graphics 
terminal.  The  capabilities  of  PLTTRNS  are  enhanced  by  ITOC, 
which,  among  other  tasks,  converts  integer  files  to  complex 
files  usable  by  PLTTRNS.  See  the  source  code  listing  in  the 
appendix  for  further  information  on  the  use  of  ITOC. 

The  plots  generated  by  ITOC  and  PLTTRNS  give  a  rough 
idea  of  the  success  of  the  global  search.  PEAK  gives 
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FIGURE  12:  FLOWCHART  OF  CORRELATION  INPLEHENTATION. 


quantitative  information  on  up  to  ten  supra-threshold  corre¬ 
lation  plane  peaks:  location,  width,  length,  and  value. 
More  peak  information  is  generated  internally  by  PEAR  if 
further  peak  discrimination  is  desired.  As  the  correlation 
functions  generally  are  not  monotonically  non-decreasing 
functions  toward  the  absolute  peak  value,  the  selection  of  a 
threshold  value  is  not  a  straight-forward  task. 

DISTANCE  uses  the  peak  locations  found  by  PEAK  to 
calculate  a  set  of  distance  factors  between  the  template  and 
the  scene  window.  The  distance  factors  are  based  on  the  LI 
and  L2  norms,  and  take  into  account  the  scene  and  template 
energies,  the  template  average  pixel  value,  and  the  window 
size.  The  factors  can  be  computed  for  the  256x256  pixel 
original  images,  or  the  128x128  pixel  reduced  images  used  in 
the  correlation  process.  The  factors  are  scaled  into  a  0-100 
range.  See  Figure  13  for  a  flowchart  of  the  evaluation 
process. 

SUPPORT  SUBROUTINES 

Several  subroutines  are  common  to  many  programs;  they 
will  be  briefly  described  below. 

TIMER  is  used  to  measure  the  execution  time  of  a  program 
in  hours,  minutes,  and  seconds.  The  first  call  to  TIMER 
starts  the  "stopwatch,"  and  the  second  call  stops  it.  The 
run  time  is  printed  on  the  console. 
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FIGURE  13’  FLOWCHART  OF  EVALUATION  PROCESS. 


IOF  is  used  to  allow  switches  and  filenames  to  be  read 
in  the  execution  command  line.  IOF  is  a  powerful  device  when 
macro  (or  automatic)  programs  are  run. 

REPACK  is  the  routine  to  convert  unpacked  arrays  into 
packed  arrays  before  they  are  written  to  disk.  UNPACK 
converts  packed  arrays  read  from  disk  into  unpacked  to  be 
processed. 

XROBLK  is  used  to  read  a  packed  block  (256  words)  from 
disk,  unpack  it,  and  return  the  unpacked  array  (1024  words) 
to  the  calling  program. 

In  the  next  chapter,  observations  will  be  made  on  some 
results  of  the  correlation  process. 


VI.  OBSERVATIONS 


In  this  chapter,  the  behavior  of  the  detection  process 
is  briefly  discussed,  and  demonstrated  with  several 
appropriate  examples.  In  particular,  the  problems  of  dealing 
with  the  clutter  energy,  variable  target  illumination,  window 
positioning  errors,  and  window  classification  are  considered. 
(See  the  appendix  for  more  test  results.) 

CLUTTER  ENERGY 

The  results  of  some  early  testing  led  to  several 
modifications  of  the  correlation  process.  Consider  the 
cross-correlation  functions  (CCF)  between  a  template  and  a 
scene  with  a  target  in  it,  shown  in  Figure  14,  and  the  CCF 
between  a  template  and  a  scene  without  a  target,  shown  in 
Figure  15.  The  template  and  both  scenes  have  energies  of 
unity  (globally  normalized) .  In  neither  case  are  there 
clearly  defined  peaks  in  the  CCF  to  suggest  possible  target 
locations.  The  maximum  function  values  in  both  cases 
corresponded  to  windows  of  light  background  only,  rather  than 
to  "interesting"  objects  such  as  trees  or  tanks.  This 
behavior  is  attributed  to  the  dominance  of  the  high-energy 
background.  Relatively  low  energy  objects  will  not  correlate 
well  with  the  template,  regardless  of  their  forms.  This 
failure  is  a  classic  weakness  of  the  non-normal ized  CCF. 
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To  take  advantage  of  the  a  priori  knowledge  that  the 
object  of  interest  will  have  low-energy  content,  the  CCF's 
were  computed  using  the  negatives  of  the  scene  and  template. 
These  are  shown  in  Figures  16  and  17.  The  peaks  of  the 
functions  are  much  more  distinct  than  those  of  the  positive 
CCF's.  Notice  especially  the  center  peak  of  Figure  16,  which 
corresponds  to  a  tank  in  the  scene. 

The  negative  CCF's  appeared  to  be  an  improvement,  but 
two  problems  still  remained.  First,  the  digitizer  noise  of 
columns  1-8  resulted  in  the  line  of  peaks  seen  at  the  right 
hand  side  of  the  negative  CCF's.  This  line  can  clearly  be 
seen  in  the  top  20%  contour  plot  in  Figure  18.  The  second 
problem  was  that  the  false  peaks  in  Figure  16  had  values  as 
large  as  the  true  peak  (the  true  peak  is  the  single  peak 
corresponding  to  the  target) .  One  method  to  overcome  these 
problems  is  to  investigate  each  candidate  peak  by  recomputing 
the  CCF  for  a  smaller  sized  scene  window  corresponding  to 
that  peak.  Each  new  CCF  is  then  investigated  for  the  occur¬ 
rence  of  distinct  peaks,  and  the  process  continues  until  one 
window  in  the  scene  is  chosen  as  being  most  likely  to  con¬ 
tain  a  target.  Then  a  decision  is  made  as  to  the  content  of 
the  window. 

The  process  of  iteratively  "windowing-in"  on  a  target 
can  be  successful  because  the  target-to-clutter  area  ratio 
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FIGURE  16:  NEGATIVE  CCF  FOR  TANK  SCENE  TEST.VD  AND  TEMPLATE  H3 
TARGET-TO-CLUTTER  AREA  RATIO  IS  0.12. 


FIGURE  18:  TOP  20Z  CONTOUR  PLOT  OF  FIGURE  17. 


increases  with  each  successive  CCF  computed.  Consider  Figure 
16,  in  which  the  target-to-clutter  area  ratio  was  0.12. 
Notice  (in  Figure  19)  the  dramatic  improvement  obtained  by 
increasing  the  target-to-clutter  area  ratio  to  0.40. 

An  alternate  method  of  choosing  candidate  windows  is 
desired,  as  the  windowing  process  can  become  computationally 
prohibitive  for  large  scene  areas.  For  this  reason,  the  scene 
was  locally  normalized  by  grid  blocks,  in  hopes  of  obtaining 
a  distinct  true  peak  with  just  one  CCF  calculation.  In  all 
test  cases  a  4x6  normalization  grid  was  used.  Rows  121-128 
and  columns  1-8  of  the  reduced  scenes  were  set  to  zero  prior 
to  normalization  to  remove  the  possibility  of  digitizer  noise 
dominating  the  correlation  function.  The  effect  of  these  two 
changes  can  be  seen  by  comparing  Figure  20  with  Figure  17. 


of  the  score  factor  as  the  search  window  is  shifted.  The 
scene  tested  did  not  have  a  target  present  in  this  case. 


TARGET  ILLUMINATION 

The  grid  method  of  normalization  was  next  tested  on 
scene  PTANKH3.  This  is  the  same  scene  that  was  used  to 
create  template  H3  (PTEMPH3 ) .  The  scene  was  re-digitized  to 
test  the  effect  of  the  noise  added  in  the  digitization 
process.  See  Figure  21  for  the  digitized  version  of  PTANKH3. 


FIGURE  21:  DIGITIZED  TANK  SCENE  H3. 


The  CCF  between  PTANKH3  and  PTEMPH3  is  shown  in  Figure 
Unfortunately,  there  are  no  distinct  peaks  present  to 
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suggest  further  investigation.  Apparently  the  grid  rectangle 
size  was  not  small  enough  to  sufficiently  normalize  the  high 
energy  non-target  areas.  For  a  4x6  grid,  the  grid  rectangle 
size  is  30x20  pixels,  for  an  area  of  600  pixels.  A  finer 
grid  may  produce  a  more  meaningful  CCF  for  the  case  when  the 
target  average  energy  per  pixel  is  less  than  the  clutter 
average  energy  per  pixel.  The  effect  of  varying  the  grid 
size  was  not  studied  in  this  thesis. 


FIGURE  22:  CCF  OF  PTANIH3  AND  PTEHPH3. 


The  CCF  of  the  negative  of  PTANKH3  (NTANKH3)  and  NT EM PH 3 
did  exhibit  a  sharp  peak,  with  a  maximum  value  at  (129,129) . 
This  coordinate  pair  corresponds  to  a  registration  shift 
error  of  only  one  pixel  in  each  shift  direction.  See  Figures 
23  and  24  for  the  3-D  and  contour  plots.  The  peaks  found, 
with  corresponding  distances,  are  listed  in  Table  IV.  The 
distance  computed  for  Peak  #1  from  the  reduced  scene 


comparison  indicated  the  best  match  for  a  shift  of  (128,128) , 
as  was  expected.  The  original  256x256  scenes  were  compared 
rising  (128,128)  as  the  peak  location,  and  a  distance  score  of 
35  was  computed. 


TABLE  IV.  PEAK  EVALUATION  OP  NPRQIH3,  KITH  CORRESPONDING  DISTANCES 
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One  of  the  reasons  for  this  low  score  can  be  attributed 

to  the  method  of  normalization  used  prior  to  the  distance 

* 

calculation.  As  the  template  window  and  the  target  have  been 
given  to  be  the  same  size,  the  search  window  will  contain 
only  target  information  when  properly  placed.  The  clutter 
energy  will  be  eliminated,  so  the  only  normalization  thought 
to  be  needed  was  search  window  normalization  (in  other  words, 
a  lxl  grid  normalization) .  This  type  of  normalization  will 
not  take  into  account  energy  changes  over  the  target  due  to 
illumination.  Recall  that  the  brightness  function  is  a 
product  of  the  reflectance  and  illumination  functions.  It  is 
an  implicit  assumption  that  the  target  reflectance  function 
is  expected  to  be  about  the  same  as  that  of  the  template.  To 
take  into  account  the  variability  of  the  illumination,  both 
the  scene  and  the  template  should  be  normalized  by  a  method 
similar  to  that  used  in  computing  the  CCF.  Grid  normali¬ 
zation  to  improve  the  confidence  of  the  distance  factors  was 
not  implemented  due  to  time  constraints.  However,  a  simple 
demonstration  of  the  advantage  of  using  search  window  grid 
normalization  is  given  in  Table  V.  Note  the  increase  in  the 
best  distance  score  factor  from  35  to  79. 

POSITIONING  ERRORS 

The  correlation  process  was  next  tested  using  NTANKD2  as 
the  input  scene,  and  NTEMPH3  as  the  template.  The  results  of 
the  peak  evaluation  of  NPR0DD2  are  given  in  Table  VI.  Tank 
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scene  D2  and  the  window  choices  of  PEAK  are  shown  in  Figures 
25  and  26.  The  distances  corresponding  to  the  windows  chosen 
by  PEAK  are  given  in  Table  VII.  Notice  that  the  first  window 
choice  clearly  corresponds  to  the  target,  but  that  the 
distance  factors  favor  the  second  window  choice.  As 
expected,  the  distance  factors  are  very  sensitive  to  window 
positioning  errors,  since  they  are  based  only  upon  a 
pixel-by-pixel  comparison. 


To  achieve  smaller  window  positioning  errors,  a  second 
CCF  can  be  computed  between  the  template  and  an  intermediate 
size  window.  The  window  size  would  be  larger  than  the 
template,  but  much  smaller  than  the  original  scene  area.  An 
improvement  of  this  type  would  be  relatively  computation¬ 
ally  inexpensive,  as  the  FFT  based  correlation  time  is 
proportional  to  NxN  log  N,  where  N  is  the  window  width. 


FIGURE  25:  DIGITIZED  TAME  SCENE  D2. 


FIGURE  24:  FIRST  TWO  WINDOW  CHOICES  FROH  PEAI 


EVALUATION  OF  NPRODD2. 
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TABLE  VI.  PEAK  EVALUATION  OF  HPR0BB2 


TABLE  VII.  BISTANCES  CORRESPONBINC  TO  UINBOUS  OF  FIGURE  26 
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Thus,  the  correlation  computation  time  quickly  diminishes  as 
N  is  decreased. 

CLASSIFICATION  ERRORS 

The  ability  of  the  distance  factors  to  classify  windows 
was  tested  using  PTANKD2 ,  and  the  standard  template,  PTEMPH3 . 
A  window  was  selected  in  PTANKD2  to  correspond  to  the  tank 
only  (the  ideal  choice  of  PEAK  evaluating  NPRODD2 ) . 

None  of  the  distance  factors  of  the  nine  windows 
investigated  by  DISTANCE  were  greater  than  zero,  a  somewhat 
disappointing  result.  There  are  several  possible  explana¬ 
tions  for  this  result.  The  first  is  that  scene  D2  and 
template  H3  were  digitized  in  separate  sessions.  As  a 
result,  the  size  and  orientation  of  the  tanks  will  slightly 
differ.  The  second  explanation  is  that  the  normalization 
approach  used  in  the  distance  calculation  admittedly  does 
not  account  for  illumination  changes  over  the  target,  as 
mentioned  previously.  Finally,  the  16  level  quantization 
cannot  accurately  represent  the  continuous  image  unless  the 
brightness  function  histogram  is  reasonably  "spread  out." 

The  third  explanation  requires  further  discussion. 
Consider  the  histograms  of  template  H3,  and  of  the  tank  of 
scene  D2  (Figures  27  and  28) .  The  histogram  of  template  H3 
is  almost  evenly  distributed  from  levels  0  to  9.  On  the 
other  hand,  the  nistogram  of  the  tank  of  scene  D2  shows  63% 
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FIGURE  2A:  HISTOGRAM  OF  TANK  D2 

(BACKGROUND  PIXEL  VALUES  SUPPRESSED). 

portions  of  the  tanks,  and  also  the  bottom  portions.  Unless 
an  elegant  method  of  normalization  is  used,  any 

point-by-point  based  distance  measure  may  give  misleading 
results. 


There  are  several  methods  of  computing  distances  that 
take  into  account  the  spatial  relationships  of  the  pixels  in 
an  image.  Although  none  were  implemented  in  this  study,  they 
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FIGURE  30:  TONER  TANK  02  (LEFT)  MB  TEMPLATE  H3  (RIGHT). 


are  suggested  as  logical  improvements  to  the  detection 
process.  One  method  to  be  tested  is  to  compare  in  some  way 
the  filtered  Fourier  components  of  the  normalized  scene 
window  with  those  of  the  normalized  template  window.  Another 


possibility  to  explore  is  to  compute 
correlation  measure  of  the  form 


statistical 


Ng(mfn) 


f {s(xry) }f {t(x-m,y-n) } 


[f{s(x,y)}J 


x  y 


[f {t(x-mry-n) }] 


x  y 


where  f{  •  }  indicates  some  3x3  mask  operation  of  the  array , 
most  likely  an  edge  enhancer.  Linear  edge  enhancers  include 
the  Laplacian  mask,  while  non-linear  enhancers  include  the 
Kirsch  and  Roberts  operators  [5:  482-491].  A  3x3  mask 
operator  can  easily  be  implemented  by  modifying  the 
subroutine  TEST3  of  the  program  REMOVE. 


VII.  CONCLUSION 


SUMMARY 

In  this  thesis,  a  method  for  two-dimensional  pattern 
recognition  was  developed  and  tested.  The  method  included  a 
global  search  scheme  for  candidate  targets,  based  on  high 
speed  cross-correlation  between  a  normalized  scene  and 
template.  A  target  classification  measure  dependent  on  the 
normalized  LI  and  Euclidean  distances  was  also  presented. 

Several  computer  programs  were  written  to  carry  out  the 
process,  from  image  input  to  target  classification. 
Especially  significant  was  the  program  PEAK,  which  performs  a 
search  of  two-dimensional  cross-correlation  functions  for 
isolated  supra- threshold  peaks.  Another  program,  DISTANCE, 
computes  a  classification,  or  similarity  score  between 
template  and  scene  windows  corresponding  to  the  correlation 
peaks.  Simple  image  procesing  programs  were  also  written  and 
described. 

CONCLUSIONS 

1.  The  global  search  for  candidate  targets  by  cross¬ 
correlation  template  matching  can  be  counted  on  to  find 
candidate  targets  to  only  a  limited  degree.  The  search 
scheme  presented  in  this  study  is  best  used  to  determine  the 
approximate  location  of  amorphous  "blobs." 

2.  The  reliability  of  cross-correlation  in  finding  targets 
increases  as  the  target-to-clutter  area  ratio  increases. 


This  property  can  allow  the  target  to  be  found  by  an 
iterative  process  of  correlationss  using  smaller  and  smaller 
Window  sizes. 

3.  Similarity  measures  based  on  pixel-by-pixel  comparisons 
are  sensitive  to  slight  mis-registration  errors ,  and  to  minor 
degradations  in  the  brightness  function  due  to  illumination 
and  digitization  noise. 

4.  The  grid  normalization  method  of  chapter  four  for  use  in 
the  correlation  process  is  not  a  suitable  approximation  of 
search  window  normalization.  The  method  has  difficulties  when 
a  grid  rectangle  covers  an  area  which  contains  both  high  and 
low  energy  objects.  However,  it  works  very  well  for 
normalizing  scenes  and  templates  in  the  classification 
process,  after  a  search  window  has  been  chosen. 

RECOMMENDATIONS 

1.  An  algorithm  to  automatically  evaluate  the  information 
generated  by  the  program  PEAK  needs  to  be  developed. 
Specifically,  a  simple  rule  to  determine  the  selection  of  a 
threshold  value  (or  possibly  several  values  for  a  given  CCF) 
needs  to  be  determined.  Also  to  be  determined  is  a  method  to 
recognize  a  false  peak  (one  not  corresponding  to  a  target) , 
given  its  width,  length,  and  percentage  of  maximum  peak 
value.  Perhaps  the  cross-section  of  the  peak  could  be 
compared  with  the  rross-section  of  the  autocorrelation 
function.  Elimination  of  false  peaks  from  further 
consideration  would  result  in  computational  savings  in  the 


classification  process. 

2.  A  more  elegant  normalization  scheme  is  needed  for  use  in 
t'he  global  search  process.  One  possibility  would  be  to 
compute  a  finer  array  of  normalization  constants  that  take 
into  account  the  energies  of  the  surrounding  grid  rectangles 
(see  Figure  31) .  The  normalization  coefficients  computed 
would  essentially  be  the  result  of  overlapping  the  grid 
rectangles. 

3.  An  intermediate  correlation  using  a  smaller  window  area 
should  be  computed,  using  the  peak  information  from  the 
global  search  CCF  to  determine  the  window  size  and  center. 
This  intermediate  correlation  will  allow  for  a  precise 
positioning  of  the  window  used  in  the  computation  of  the 
similarity  measure. 

4.  The  use  of  grid  normalization  of  scene  and  template 
windows  for  the  classification  process  needs  to  be  further 
developed  and  tested.  Preliminary  results  indicate  that  grid 
normalization  can  be  used  successfully  to  account  for 
differences  in  energy  between  the  template  and  scene  tanks, 
as  long  as  most  of  the  clutter  energy  is  eliminated  (which  is 
the  case  when  the  search  window  is  accurately  positioned) . 

5.  The  choice  of  features  to  be  compared  needs  to  be 
studied.  Instead  of  comparing  the  pixel  values  of  the  scene 
candidate  window  directly  with  those  of  template  window,  the 
comparison  of  the  high  or  low-frequency  Fourier  components, 
for  example,  may  lead  to  more  promising  results.  It  might 
also  be  interesting  to  compute  the  distances  between  scenes 
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FIGURE  31:  SUGGESTED  NORMALIZATION  SCHEME 
Consider  Che  case  where  Che  energies  of  9  local  recCangular 
"whole”  regions  A-I  are  compuced. 
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From  Chese  9  energy  Ceras,  we  wane  Co  compuCe  81  normalizaClon 
coefficients.  The  coefflclenCs  for  any  "sub"-reglon  can  be  approx- 
imaced  by  a  welghCed  average  of  Che  whole  region  coefflclenCs  cor¬ 
responding  Co  Che  center  sub-region  and  its  eight  nearest  neighbors; 
where  Che  weights  are  determined  by  Che  distance  from  Che  center 
sub-region. 
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For  s(x,y)  in  sub-region  A',  one  possible  method  of  computing  < 
normalization  coefficient  is  by.  Eq.  50  (let  N(A’)  represent  the  nor 
malization  coefficient  for  region  A'): 


N(A') 


y  KlxA  +  K2x[(A  +  A  +  B  +  D)  +  (A+B  +  E  +  D  ] 


(50) 


This  normalization  method  may  avoid  some  of  the  severe  discontinuities 
sometimes  created  by  whole  region  normalization. 


and  templates  that  have  been  mask  processed  (for  example  edge 
enhanced) ,  as  mentioned  in  chapter  six. 
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APPENDIX  Ai 
USE.  QL  macro  piles 


Under  the  RDOS  operating  system,  macro  (or  indirect 
files)  can  be  run  to  execute  a  series  of  CLI  commands.  Thus, 
a  process  consisting  of  a  series  of  program  can  be  run 
automatically  as  long  as  none  of  the  programs  require 
interactive  user  input.  In  order  that  file  names  and  program 
options  can  be  specified  in  the  execution  line,  all  of  the 
programs  in  the  correlation  process  use  the  COMARG  call  to 
read  the  command  line  argument  string.  This  appendix 
describes  the  use  of  several  macro  programs. 

NTEMPCP.MC 

Macro  NTEMPCP.MC  creates  NTEMP.CP,  the  Fourier  transform 
of  the  reduced  template,  RNTEMP.VD.  RNTEMP.VD  is  the 
negative  of  the  256x256  PTEMP.VD,  reduced  to  the  lower  right 
quadrant  by  REDUCE  and  NMOVE.  Links  to  NORMALIZE  and  DIRECT 
must  exist  before  NTEMPCP.MC  is  executed.  The  program  lines 
are  as  follows: 

NORMALIZE/L  RNTEMP.VD  NTEMP.CP 
DELETE/V  RNTEMP.  VD 
DIRECT  NTEMP.CP/I  256/N 


Execution  of  NPRODIN.MC  results  in  the  CCP  between  the 
template  and  SCENE. VD  to  be  computed  and  stored  in  NPROD.IN. 
Piles  HOLD1.CP  and  H0LD2.CP  should  exist  as  1024-block 
contiguous  files;  they  will  be  created  if  they  do  not  exist. 
The  use  of  contiguous  (as  opposed  random)  files  decreases  run 
time  by  a  factor  of  3.  The  following  files  must  be  linked  to 
before  NPRODIN.MC  can  be  run  (the  links  in  CROMER. DR  are 
given  as  an  example) : 

CROMER 


: 26. 26 

HOLD1.  CP 

STROUPE: HOLD1.  CP 

H0LD2.  CP 

STROUPE:  HOLD2.  CP 

NTEMP .  CP 

STROUPE:  NTEMP.  CP 

SCENE.  VD 

STROUPE: SCENE.  VD 

NEGATE.  SV 

STROUPE:  NEGATE.  SV 

REDUCE.  SV 

STROUPE: REDUCE.  SV 

NORMALIZE.  SV 

STROUPE: NORMALIZE.  SV 

DIRECT.  SV 

DP4F: DIRECT.  SV 

CMULTIPLY.  SV 

STROUPE: CMULTIPLY.  SV 

INVERSE.  SV 

DP4F:  INVERSE.  SV 

CTO  I.  SV 

STROUPE.  CTO I.  SV 

Run  time  will  be 

less  than  10  minutes.  The 

interest  should  be  renamed  SCENE. VD  before  executing; 
NPROD.IN  should  be  renamed  after  the  macro  program  execution. 
See  COMPUTE. MC  for  further  clarification  on  the  use  of 
NPRODIN.MC.  One  version  of  NPRODIN.MC  is  as  follows: 


RENAME  HOLD1.  CP  NSCENE.  CP 
RENAME  H0LD2.  CP  NPROD.  CP 
NEOATE  SCENE.  VD  NSCENE.  VD 
REDUCE  NSCENE.  VD  RNSCENE.  VD 
DELETE/V  NSCENE.  VD 


NORMALIZE/U  RNSCENE.  VD  NSCENE.  CP 

DELETE/V  RNSCENE.  VD 

DIRECT  NSCENE.  CP/ 1  256/N 

CMULTIPLY  NTEMP.  CP  NSCENE.  CP  NPRQD.  CP 

RENAME  NSCENE.  CP  HOLD1.  CP 

INVERSE  NPROD.  CP/ 1  256/N 

CTQI  NPROD.  CP  NPROD.  IN 

RENAME  NPROD.  CP  H0LD2.  CP 


COMPUTE. MC  is  used  to  automatically,  without  user 
interaction,  compute  the  CCF's  corresponding  to  six  different 
scenes  and  one  template.  The  total  run  time  is  less  than  one 
hour.  The  macro  controls  the  re-naming  (or  re-linking  in  this 
case)  of  the  dummy  files  SCENE. VD  and  NPROD. IN.  The  CCF's 
created  are  moved  to  other  disks  to  allow  room  to  run  the 
rest  of  the  correlations.  The  macro  program  requires  256K 
bytes  of  free  disk  (assuming  NPROD. IN  exists).  If  for  any 
reason  NPROD. IN  cannot  be  moved,  it  is  simply  overwritten. 
The  program  may  be  aborted  at  any  time  the  user  desires; 
results  up  to  that  time  will  be  saved.  The  program  is  given 
as  follows: 

LINK  SCENE.  VD  CROMER: PTANK04.  VD 
NPROD IN 

RENAME  NPROD.  IN  NPR0DG4.  IN 
UNLINK  SCENE.  VD 

LINK  SCENE. VD  CROMER: PTANKD4.  VD 
NPROD IN 

MOVE/V/D/R  DP5F  NPROD.  IN/S  NPR0DD4.  IN 
LINK  NPR0DD4.  IN  DP5F: NPR0DD4.  IN 
UNLINK  SCENE.  VD 

LINK  SCENE.  VD  CROMER:  PTANKC3.  VD 
NPROD IN 

MOVE/V/D/R  DP5F  NPROD.  IN/S  NPR0DC3.  IN 
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LINK  NPR0DC3.  IN  DP5F:  NPR0DC3.  IN 
UNLINK  SCENE.  VD 

LINK  SCENE.  VD  CROMER: PTANKB2.  VD 
NPRODIN 

MOVE/V/D/R  DP5F  NPROD.  IN/S  NPR0DB2.  IN 
LINK  NPR0DB2.  IN  DP5F:  NPR0DB2.  IN 
UNLINK  SCENE.  VD 

LINK  SCENE.  VD  CROMER: PTANKD2.  VD 
NPRODIN 

MOVE/V/D/R  DP 5  NPROD.  IN/S  NPR0DD2.  IN 
LINK  NPR0DD2.  IN  DP5: NPR0DD2.  IN 
UNLINK  SCENE.  VD 

LINK  SCENE.  VD  CROMER: PTANKE2.  VD 
NPRODIN 

MOVE/V/D/R  DP 5  NPROD.  IN/S  NPR0DE2.  IN 
LINK  NPR0DE2.  IN  DP5:  NPR0DE2.  IN 
UNLINK  SCENE.  VD 
DIR  DP4 

MESSAGE  RECOGNITION  PROGRAM  COMPLETED 

GDIR 

GTOD 


The  scene  files  correlated  with  the  template  are  as  follows: 

CROMER 


PTANKB2.  VD 

32768 

PC 

PTANKC3.  VD 

32768 

PC 

PTANKD2.  VD 

32768 

PC 

PTANKD4.  VD 

32768 

PC 

PTANKE2.  VD 

32768 

PC 

PTANKG4.  VD 

32768 

PC 

The  files  and  links  created  by  this  version  of  COMPUTE. MC  are 
as  follows: 


CROMER 


NPR0DB2.  IN  DP5F:  NPR0DB2.  IN 

NPR0DC3.  IN  DP5F:  NPRODC3.  IN 

NPR0DD2.  IN  DP5:  NPR0DD2.  IN 

NPR0DD4. IN  DP5F: NPR0DD4. IN 

NPR0DE2.  IN  DP5:  NPR0DE2.  IN 

NPR0DG4.  IN  131072  C 


>  „• 

r.- 

.’i  ,* 
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APPENDIX  B: 

r 

r,' 

/v’ 

SUMMARY  OF  PROGRAM  USAGE 

•T-  /**. 

i  ^ 

INPUT  FILE 

OUTEOT  FILE 

1 

TYPE 

EXECUTION  LINE  FOFMAT 

TYPE 

[  .VD] 

VIDEO? 

[.VD] 

-> 

.VD 

QUICKAVE7 

AVERAGE?  .VD 

i 

.VD 

REMOVE  inf  lie  outl  exit  2 

.VD 

•: 

.VD 

EVIDBIST 

— 

.VD 

TONER 

.VD 

i 

.VD 

MOVE 

.VD 

.VD 

NEGME[/F]  infile  outfile 

.VD 

.VD 

REDUCE  infile  outfile 

.VD 

i 

•  * 

.VD 

NORMALIZE l/u  or  /L]  infile  outfile 

.CP 

.CP 

DIRECT  infile/I  256/N  [outfile/O] 

.CP 

.,  ■ 

.CP 

CMJLTIHjY  infilel  infile2  outfile 

.CP 

1  & 

.CP 

INVERSE  infile/I  256/N  [outfile/O] 

.CP 

.CP 

CPOI  infile  outfile 

.IN 

.IN 

IffJLTIHiY  infilel  infile2  outfile 

.IN 

I 

.IN 

IT0C/[A,N,E,H  or  0]  infile [/C]  [A]  outfile  [/M] 

.CP 

y 

.CP 

PLTTRNS  infile/I  256/N 

— 

* 

a 

.IN 

PEAK 

— 

.VD 

DISTANCE 

— 

> 

.CP 

CTOV 

.VD 

1 

.VD 

PICTURE 

— 

where 

.VD  —  32K  bytes  packed  video 

,n 

1 

.IN  —  128K  bytes  integer 

.CP  —  512K  bytes  oonplex 

»• 

,* 

[  ]  —  optional  input  or  output 

i 
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This  appendix  contains  the  following  programs  (subroutines 
given  in  parenthesis) : 

1.  VIDE07  (CHAN7,  VABORT,  ERCHK) 

2.  QU1CKAVE7 

3.  PICTURE  (OUT3X3,  OUT4X3) 


c******  NOVA  VIDEO?  INPUT/OUTPUT  ROUTINE  VIA  CRCMEMQD  COMPUTER  ***** 


12  Aug  1982 


This  program  allows  the  user  to  display  a  video  file 
repeatedly  any  number  of  times.  It  also  allows  the 
user  to  input  or  output  video  files  named  AD,  Al,  A2, 
A3,  A4,  A5,  etc.  automatically.  These  files  nay 
then  be  averaged  by  QUICKAVE7. 

Execution  Line  Format: 

VTDE07 

Load  Line  Format: 

RLDR  VIDEO?  CHAN7  ERCHK  VABORT  CHANNEL 
DCHRX  DCHTX  SANDS  CANDR  PORT.  LB 


q  ***** 

C  WRITTEN  BY  Lt.  Jim  Cromer  12  Aug  1982  ***** 

C  Fortran  IV  ***** 

c  ***** 

C'  This  program  allows  the  user  to  display  a  video  file  ***** 

C  repeatedly  any  number  of  times.  It  also  allows  the  ***** 

C  user  to  input  or  output  video  files  named  AD,  Al,  A2,  ***** 

C  A3,  A4,  A5,  etc.  automatically.  These  files  may  ***** 

C  then  be  averaged  by  QUIO&VE7.  ***** 

C  ***** 

C  Execution  Line  Format:  ***** 

C  VTDE07  ***** 

Q  ***** 

C  Load  Line  Fornat:  ***** 

C  RLDR  VIDEO?  CHAN7  ERCHK  VABORT  CHANNEL  ***** 

C  DCHRX  DCHTX  SANDS  CANDR  FORT.  LB  ***** 

Q  ***** 

(]******************************************************************** 

DIMENSION  IRAR(2)  ,IHQID(7) 

INTEGER  FILE (7) 

IPAR(2)=0 

C 

C 

C********  nsf!R  PARAMETER  INPUT  ************************************ 

C 

TYPE  "NOTICE:  CHOPS  must  be  running  K12>" 

ACCEET'Input  or  output  (IN-0/0UT-1) ?  ",IDIR 
IF(IDIR.IE.0.AND.IDIR.I®.1)G0  TO  8  ; Error  checking 

1  ACCEPT"Enter  time  (SEC.):  ",ITIME  jMonitor  display  time 
IF(IDIR.EQ.O)QO  TO  4 
ACCEPT  "Type  1  to  output  AO-An:  ",IDIR 
IF  (IDIR.EQ.l)QO  TO  4 

ACCEPT"What  is  the  name  of  the  data  file  (13  Char  max)?  " 

READ  (11 ,8)  FILE  (1)  ;Video  frame  filename 

ACCEPT"Enter  4  times  to  be  displayed:  ",11 
IF  (FILE  (1)  .IE.  10752)  GO  TO  42  ;An  ***  means  to 

DO  2  1*1,7  ;run  with  the  last 

2  FILE(I)*nm>(I)  jfile  used 

GO  TO  42 


(}****«**«  INPUT  OR  OUTPUT  FILES  A0_A6  ******************************* 
C 

4  IPAR(1)=ITIME 

CALL  CHAN7(IDIR,IPAR) 

GO  TO  6 


C********  DISPIAY  VIDEO  FILE  OF  USER'S  CHOICE  *********************** 
C 

42  DO  43  KK«1,II 

""YPE’CHECK  MONITOR  -  -  ",KR 
CALL  VABORT  (NISK,  IM,  IPCNT,  IDCNT) 


n  o 


IDUR=1 

IPAR(l)  “ITD® 

TEppafl 

CALL  CHAJ®CL(OTSK,IDIR,IM,IPQn?fIDCNT,FILE,64f0,IPARfIERR,ISYS) 
43  CALL  EICHK(IERRf  IDIR,  IDCNT,  IPQ/T,  ISYS)  ; identify  errors 


C********  USER  INPUT  *********************************************** 

C 

6  TXPE"<12>" 

DO  7  >1,7  ;Save  filename 

7  IHOLD ( I ) =FILE ( I )  ;for  next  loop 

ACCEPT  "What  next  ( INPUT-0, OUTPUT-1, STOP-2) ?  ",  IDIR 

IF  (IDIR.BQ.O)  00  TO  1 
IF(IDIR.BQ.l)  GO  TO  1 

IF  (IDIR. BQ. 2)  STOP  "Type  VIDE07  to  rerun." 

8  FORMAT (SI 3)  ; Filename  input  format 

TYPE"<7>Irput  error;  try  again." 

GO  TO  6 
END 
C 

C************  program  VIDE07  **************************************** 


oonoooonnnoooo 


SUBROUTINE  CHAN7 (IB, IPAR)  ;by  Lt  Jin  Cromer 

q********************************************************************** 


Subroutine  CHAN7  will  digitize  a  picture  seven  tines,  and 
will  store  the  video  data  in  files  A0-A6  (i.e.  in  a  format 
usable  by  QUICKWE7) .  It  will  also  display  the  digitized 
pictures  consecutively  on  the  video  monitor. 

(Called  by  VIDE07) 

PARAMETERS  PASSED: 

IB=0 - >  input  pictures  from  camera 

IB=1 - >  output  pictures  to  the  video  monitor 

IPAR(l)  — >  display  time  in  seconds 
IPAR(2)  — >  unused 

C********************************************************************** 

INTEGER  FILE  (7)  ,IBAR(2) 

DO  100  1-1,7 
100  FILE  (I)  *0 
IERR-0 

IF(IB.BQ.l)QO  TO  300 

DO  200  1-16688,16694  ;if  digitizing,  delete  A0-A6 
FILE  (1) -I 

200  CALL  DELETE  (FILE) 

INU*=7 
GO  TO  400 

300  TYPE "Enter  the  nunber  of  pictures  to  be  displayed." 

TYPE  "They  must  be  named  AO-A(N-l)  to  be  displayed." 

ACCEPT"  Enter  N  - >  ",INIM 

IF  (INUM.GT.  10)  TYPE  "Sorry.  The  maximum  nunber  of  pictures 
$  is  iO • * 

IFdNWI.GT.10)  IN»^10 

400  DO  500  I=1,INUM  ;if  1=1,  then  FILE="A0" 

; if  1=2,  then  FILE="A1" 

FILE (1) =16687+1  ;and  so  on 
CALL  VABOKT(IA,IC,ID,IE) 

WRITE  (10, 600) FILE (1) 

CALL  CHANNEL (IA, IB, IC, ID, IE, FILE, 64, 0,IPAR,IERR,ISYS) 
500  CALL  ERCHK(IERR,IB,IE,ID,ISYS) 

600  FORMAT ("  check  monitor  —  Picture  being  displayed  — >  ",S13) 
RETURN 
END 
C 

C***********  subroutine  CHAN7  ***************************************** 


onnon 


SUBROUTINE  VABOKT  (NTSK,  IM,  IPCNT,  IDCNT) 

C*  **************************************************************** 

This  subroutine  sends  an  abort  call  to  CHANNEL, 
and  resets  the  calling  parameters  for  the  next 
VIDEO?  call. 

C***************************************************************** 

J=0  ;dumy  variable 

K=3  ;mode  3  — >  Abort 

CALL  CHANIGL(J,J,K,J,J,”A",J,J,J,IE,IS) 

IDCNT=4 

IPCNT=1  ; parameter  count 

NESK=3 

IM=2  ;mode  2  — >  video 

RETURN 
END 
C 

C*************  subroutine  VABORT  ********************************* 


n  nnnn 


SUBROUTINE  ERCHK(IERR,IDIR,IDCNrf  IPCNT,ISYS) 

C********************************************************************* 


This  subroutine  checks  for  errors  made  during  an  attempted 
call  to  CHANNEL,  and  prints  them  to  the  screen. 

(Called  by  VIDEO?) . 

C********************************************************************* 

INTEGER  CROERRfHXDNT 
LOGICAL  BTEST 

IF(IERR.EQ.O.OR*  (IERR.EQ. 13323  .AND.IDIR.EQ.O)  )G0  TO  6  ;ok 
IF  (IERR.EQ. 13323  .OR.  IERR.BQ.-24832  .OR.  IERR.  EQ. -24064. OR.  IERR 
*.EQ.-22528)GO  TO  5  ;Specific  error  messages  will  be  given 

IF (BTEST (IERR, 15)) GO  TO  9  ; Abnormal  return 

GO  TO  10  ;Error  without  abort 

5  IF  ( IERR.EQ. -24832) TYPE" <7>ABORT— FILE  DOES  NOT  EXIST" 

IF  (IERR.EQ. 13323. OR.  IERR.EQ. -240 64)  TYPE  "<7>ABORT — FILE  DOES  NOT 
*  CONTAIN  VIDEO" 

IF  (IERR. EQ.-22528)TYPE"<7>ABORT— FILE  CANNOT 
6  TYPE"<12>Channel  cleared" 

GO  TO  20 

9  TYPE"<7><12>  ABORT  INITIATE)  111  <12>" 

10  CROERR=15  .AND.  IERR 

FDO0NP=ISHFT<240  .AND.  IERR, -4) 

NOVERR=ISHFT  (-256  .AND.  IERR,-8) 

TYPE"  CRCMEMOO  ERROR  RETURNED:  ",CROERR 
TYPE"  PARAMETER/DATA  (DUNT  RETURNED: ",HXDNT 
CALL  BCLR(NOVERR,7) 

TYPE"  NOVA  ERROR  RETURNED: ",NOVERR 
TYPE"  ERROR  CODE  RETURNED: ",  IERR 
TYPE"  DATA  COUNT: ",  IDCNT 
TYPE"  PARAMETER  COUNT: ",1PCNT 
TYPE"  SYSERR  RETURNED: ", ISYS 
PAUSE 

20  RETURN 
END 
C 

0**************  subroutine  ERCHK  ************************************* 


BE  CREATED" 


;There  are 
;two  error 
;  codes  in  the 
;variahle 
;'IERR' 


;Error 

; messages  are 
;printed  for 
;user  information 
;and  correction 


C********************************************************************** 

c 

C  Program  QUICKAVE7  Written  by  Lt.  Jim  Cromer 

C  Fortran  5  (by  DATA  GENERAL) 

C 

C' 

C  This  program  will  average  the  packed  video  files 

C  "AO"  -  "A6"  pixel  by  pixel#  and  output  the  averaged  picture 

C  to  the  packed  video  array  "AVERAGE7 . VD" .  [A  packed  video 

C  array  contains  a  256x256  4-bit/pixel  picture  in  a  64  block  file 

C  on  disk#  where  1  block=256  16-bit  words  (i.e.  1  blocks  video 

C  rows) .  Each  16-bit  word  holds  4  pixels.  Video  files  are 

C  stored  on  disk  in  packed  form  to  conserve  disk  space.  File 

C  size  is  32K  bytes.]  Total  run  time  is  less  than  one 

C  minute  clock  time. 

C 

C  Execution  Line  Format:  (run  SL0WAVE7  on  the  NOVA) 

C  QUICKAVE7 

C 

C  Loader  Command  Line  Format: 

C  KLDR  QUICKAVE7  TIMER  UNPACK  REPACK  20/C  @FUB@ 

C  The  20/C  opens  up  enough  channels  to  run  the 

C  program  (by  default#  only  8  channels  (designated  and  pre- 

C  designated)  are  normally  available] .  Links  to  AD-Afi  in 

C  DPOF  should  be  created. 

C 

C 

c********  array  management  ******************************************** 

c 

C  Arrays  A0-A6  hold  2  blocks  each  (i.e.  8  packed  video  rows) 

C 

C  Arrays  A0U-A6U  hold  8  blocks  each  (i.e.  8  unpacked  video  rows) 

C 

C  Array  AVE  is  used  for  the  unpacked  averaged  picture  (8  rows) 

C 

C  Array  AVEP  holds  the  packed  averaged  picture  (8  rows) 

C 

C  [If  array  sizes  are  modified#  note  that  arrays  A*U  must  be  diraen- 

C  sioned  to  be  4  times  larger  than  arrays  A*.] 

C 

C  The  EQUIVALENCE  statement  is  used  to  reduce  memory 

C  requirements.  Once  an  array  is  unpacked#  the  packed  array  is 

C  no  longer  used.  Therefore  unneeded  packed  arrays  can  be  used 

C  to  hold  unpacked  arrays  (but  not  at  the  same  time  l ) . 

C 

INTEGER  FILE (2) 

INTEGER  AO (512)  #A1(512)  #A2(512)  #A3(512)  #A4(512)  ,A5(512)  # 

$  A6(512),AVE(2048)#A#B#C,D 

INTEGER  A0U(2048)  #A1U(2048)  #A2U(2048)  #A3U(2048)  ,A4U(2048) , 
$A5U(2048) #A6U(2048) #AVEP(512) 

COMMON  A0U,A1U,A2U#A3U#A4U#A5U,A6U#AVEP  /COMMON  must  be  declare 

/before  EQUIVALENCE 
/can  be  used  (Fortran  IV,' 
EQUIVALENCE  (AVE#  AOU)  #  (A0#A1U)  #  (A1#A2U)  #  (A2#A3U)  #  (A3#A4U)  #  (A4#A5U) 
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$(A5,A6U),(A6,AVEP) 


TYPE”  Program  QUIC^VE7  now  executing  . 
CALL  TUdR(O)  ;start  timer 


,<15><7>" 


C********  I/O  FILE  MANAGEMENT  ****************************************** 
C 

C  Seven  channels  to  "AO"  thru  "A6"  must  be  OPENed  for 

C  reading.  An  eighth  channel  to  "AVERA3E7.VD"  must  be  OPENed  for 

C  writing  (after  "AVERAGE7  .VD"  is  created) . 


FILE  (2)  =0 
DO  998  1=0,6 

FILE (1) =16688+1  ; ASCI I  for  A_ 

CALL  OPEN(I, FILE, 1, IER) 

IF  (IER. Id.  1)  WRITE  (10, 999)  FILE  (1) ,IER 
EORMAT("  OPEN  ",S4,"  error  #",I5) 

CONTINUE 

DELETE  "AVERAGE7.VD" 

CALL  CFILW( "AVERAGE? .VD", 3, 64, IER)  ;create  a  contiguous 

IF  (IER.EQ.l)QO  TO  997  ;file,  if  possible 

CALL  CFILW ( "AVERAGE7 . VD" , 2 , IER)  jcreate  a  random  file 

;of  variable  size 

IF  (IER.  Id.  1)  TYPE  "  AVERAGE?  .VD  create  error:", IER  STOP 

CALL  OPEN  (7 ,  "AVERAGE7 . VD" ,  3 ,  IER)  jOPEN  a  channel  to 

IF  (IER.  Id.  1)  TYPE  "  OPEN  7  error  #",IER  ;  "AVERAGE7 .  VD"  for  writi-.a 


DO  4  1=0,62,2  ;average  64  blocks,  2  at  a  time 

C*********************************************************************** 

C 

C  Read  2  blocks  from  each  picture  to  be  averaged  (8  video  rows) 

C 

CALL  RDBLK(O,I,A0,2,IER) 

CALL  RDBLK(1,I,A1,2,JER) 

CALL  RCBLK(2,I,A2,2,KER) 

CALL  FDBLK(3,I,A3,2,LER) 

CALL  RDBLK(4,I,A4,2,MER) 

CALL  RDBLK(5,I,A5,2,NER) 

CALL  BDBLK(6,I,A6,2,I6ER) 

IF  (IER. Id.  1 .  .OR.iJER. Id.  1 .  .OR.KER.Id.l .  .*•  R. IER. Id.  1 .  .OR. 

$  MER.Id.l. .OR.NER.td.l. .OR.I6ER.!d.l)GO  TO  5 
C 

C  Unpack  the  arrays 
C 

CALL  UNPACK ( 51 2, AD, AO U) 

CALL  UNPACK(512,A1,A1U) 

CALL  UNPACK (512, A2,A2U) 

CALL  UNPACK (512, A3 ,A3U) 

GALL  UNPACK (512, A4 ,A4U) 

CALL  UNPACK (512, A5,A5U) 


nnnn  x  non  non 


CALL  UNEACK (512,A6,A6U) 

Perform  pixel  by  pixel  averaging 
DO  1  K=l,2048 

AVE  (K)  =IFIX  ( FLOAT  (AOU  (K)  +A1U (K)  4A2U  (K)  +A3U  (K)  +A4U  (K)  +A5U 
$+A6U(K))/7.0+0.5)  ; round  to  nearest  integer 

1  CONTINUE 

Pack  array  AVE  into  AVEP  and  write  it  to  the  disk 

CALL  REPACK  (51 2, AVE, AVEP) 

CALL  HRBLK(7,I,AVEP,2,IER) 

IF(IER.NE.1)TYPE  *  WRBLK  #"fIf*  error:  ",IER 
TYPE"  Block  averaged. <15>”  ;can  be  deleted  to 

; decrease  run  time 

4  CONTINUE 

Print  total  run  time  message  to  the  console  CRT, 
and  EXIT  the  program. 

CALL  TUGR(l)  5  stop  timer 

GO  TO  6 

5  TYPE"  <7><15>RDBLK  AO  error : " , IER, " <15>RDBLK  A1  error: ",JER, 
$"<15>RDBLK  A2  error : \KER,  "<15>RDBLK  A3  error  :",LER, 

$"<15>RDBLK  A4  error :",MER,"<15>RDBLK  A5  error :",NER, 

$"<15>RDBLK  A6  error  :\I6ER,"<7X7><15>" 

6  CALL  RESET 

TYPE"  <15>  Program  QUICKAVE7  finished." 

TYPE"  The  averaged  file  is  named  AVFRAGE7.VD<15X7X7>" 

STOP 
END 

C*************  program  QUICKAVE7  ************************************ 


nono  ooo  onn  non  nonoonnnnnononnnnn 


C********************************************************************* 


Program  PICTURE  Written  by  Lt.  Simmons  14  Oct  1981 
Fortran  5  Revised  by  Lt.  Cromer  12  July  1982 

This  program  will  convert  video  pixels  to  1  in epr inter  pixels, 
and  will  put  the  picture  in  a  file  or  to  the  Printronix  300 
linepr inter.  This  program  prints  either  the  complete  256  by 
256  pixel  picture,  or  a  smaller  picture  that  does  not  contain 
the  noise  created  by  the  video  digitizer  (the  last  five  blocks 
in  the  video  file  are  noise) .  Odd  nunbered  Video  lines  (rows) 
are  represented  by  3x3  pixels,  while  even  numbered  Video  lines 
are  represented  by  4x3  pixels.  Run  time  for  a  shortened 
picture  of  length  222  lines  is  about  1.5  minutes;  the  file 
size  will  be  about  204Kbytes  (or  102,000  words) .  For  a  256x256 
picture,  run  time  should  be  2  minutes,  with  239Kbytes  used  on 
disk. 


********************************************************************** 
INTEGER  ATITLE(40) 

DIMENSION  I1ARRAY (268)  ,I2ARRAY(264)  ,I3ARRAY(201)  ,I4ARRAY(198) 
EQUIVALENCE  (IlARRAY , I2ARRAY , I3ARRAY , I4ARRAY) 

DIMENSION  ILP(4,67)  ,IFILE(7)  ,IRBC(64)  ,ISAV(4) 

LOGICAL  SHORT, TITLE 


Set  up  solid  line,  space,  and  line  feed/plot-on  characters 


IL=177777K 

NC=40100K 

LF=012K 

LFP02412K 


;Solid  line 
;  Space 
;Line  feed 
;Line  feed/plot  on 


Set  up  parameters  for  complete  picture  display. 


SH0RT=.  FALSE. 

m=66 

N2=256 

N3»l 

N4=66 

N5=67 

N6=l 


Open  the  video  file  for  input 


;Short  picture  test 
;Top  and  bottom  border  length 
;Number  of  lines  displayed 
;Location  of  left  border 
;Location  of  right  border 
/•Location  of  line  feed 
; Length  of  lines  displayed 


ACCEPT"  What  is  the  name  of  the  input  file?  " 

READ (11, 17) IFILE(l)  ;Read  video  file  name 

CALL  OPEN (1,IFILE,1, IER)  ;Open  the  video  file 

CALL  CHECK (IER) 

IF(IER.NE.1)TYPE"  Input  open  error:", IER 


Ask  for  an  output  file,  either  the  linepr  inter  or  a  disk  file, 
and  open  the  disk  file  if  necessary. 

ACCEPT"  Do  you  want  a  disk  file  created  (Y  or  N)?  " 


1  READ  (11, 19)  N  ;Read  one  ASCII  character 

IF(N.EQ.19968)G0  TO  2  ;File  was  not  selected 

IF  (N.  IE. 22784)  GO  TO  20  ;  Input  error 

ACCEPT"  A  file  was  selected;  what  should  its  name  be 
*  (13  char  max)?  " 

READ(11,17) IFILE(I)  ;Read  output  file  name 

CALL  DFILW(IFILE,JER)  ;Make  sure  that  the 

IF(JER.fE.l .AND. JER.NE.13)TYPE"  Output  delete  error : ",JER 
CALL  CFILW(IFILE,2,KER)  ;file  does  not  exist 

IF(KER.fE.l)TYPE"  Output  create  error :"fKER 
CALL  OPEN (12,IFILE,3 ,LER)  ;Open  the  output  file 

CALL  CHECK (LER) 

IF  (LER.  IE.  1)  TYPE"  Output  file  open  error:", LER 
GO  TO  3 

2  TYPE"  The  picture  will  only  go  to  the  linepr inter." 

C 

C  Choose  between  complete  picture  and  noiseless  picture. 

C 

3  CONTINUE 

ACCEPT  "Do  you  want  to  title  the  picture  (Y  or  N)?  " 
READ(11,19)N 

IF(N.EQ.22784)TI1LE=.7HJE. 

DO  300  KT=1,40 

ATITLE(KT)  =0 
300  CONTINUE 

IF (TITLE) TYPE "Enter  title  below  (up  to  80  characters)." 

IF  (TITLE)  READ  (11 , 3000 )  ATITLE  (1 ) 

3000  FOFMAT(S80) 

ACCEPT  "Do  you  want  a  complete  picture  (Y  or  N)?  " 

4  READ (11, 19) N  ;Read  one  ASCII  character 

IF (N.BQ. 19968) GO  TO  22  ; Response  was  "NO" 

IF(N.NE.22784)GO  TO  21  ; Input  error 

TYPE"  A  complete  256  by  256  pixel  picture  was  chosen." 

C 

C  Put  a  border  at  the  top  of  the  picture. 

C 


5  DO  7  1=1,3 

IF  (SHORT)  WRITE  BINARY  (12)  NC 
IF (SHORT) WRITE  BINARY (1 2) NC 
DO  6  J=1,N1 

6  WRITE  BINARY  (12)  IL 

7  WRITE  BINARY  (12)  LFPC 
C 


; Space  right 
; Space  right  again 

;Print  a  line 
; Terminate  the  line 


C****************************************************************** 


C 

C 

C 

C 

C 

C 


Each  line  of  the  picture  will  have  a  border  on  each  end.  A 
DO-LOOP  loops  233  (or  256  for  whole  frame)  times  around  the 
next  three  program  parts 


JTEST  ■  -1 
DO  13  JA-1,N2 
JTEST  ■  JTEST  *  -1 
J j  -  4 
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f)D-fll24  788 
UNCLASSIFIED 


SCENE  ANALYSIS:  NON-LINEAR  SPATIAL  FILTERING  FOR 
AUTONATIC  TARGET  DETECTIONS)  AIR  FORCE  INST  OF  TECH 
HRIGHT-PATTERSON  AFB  OH  SCHOOL  OF  ENG I.  .  J  H  CROMER 
DEC  82  AFI T 7GE/EE782D-26  F/G  1778 


nnnnn  _  ..  non  non  onrs 


IF(JTEST.Gr.O)  JJ*3  /Odd  iteration 

Put  a  border  down  the  left  hand  side 


DO  8  K=1,JJ 

IF(SHDRT)ILP(K,1)=NC 
IF  (SHORT)  ILP(K,2)=NC 

8  ILP(K,N3)=43500K 

Put  a  border  down  the  right  hand  side 

DO  9  L=1,JJ 

ILP(L,N4)=40170K 

9  ILP(L,N5)=LFPC 


;Put  a  left  border 
;two  spaces  in 
; ("1/4"  space) 

/for  short  picture 


; insert  border  and 
;line  feed  after 

/the  picture 


Convert  the  video  picture  pixels  to  $LPT  pixels 


READ  BINARY (1)  EREC  ;Read  one  video  line 

DO  12  N=N6,64 

IWR=BYTE(IREC(N) ,2)  /Right  two  pixels 

IWL=BYTE(IRBC(N)  ,1)  /Left  two  pixels 

IF(JJ.GT.3.5)CALL  0UT4X3(IWL,ISAV) 
IF(JJ.LT.3.5)CALL  0UT3X3(IWL,ISAV) 

N7=Ntl 

IF  (SHORT)  N7=N 
DO  11  J&*1,JJ 

1  BXTE(ILP(JB,N7)  ,1)=ISAV(JB)  /Move  to  high  byte 
IF(JJ.LT.3.5)CALL  0UT3X3(IWR,ISAV) 

IF(JJ.GT.3.5)CALL  OUT4X3 (IWR, ISAV) 

DO  12  J01,JJ 

2  BXTE(ILP(JC,N7)  ,2)=ISAV(JC)  /Store  low  byte 

L=0 

DO  130  JE»1,JJ 
DO  130  JI>*1,N5 
]>L+1 

IlARRAY  (L)  “ILP  ( JE,  JD) 

30  CONTINUE 

IF(16.BQ.66..AND.JJ.BQ.3)WRITE  BINARY (12) I4ARRAY 
IF(N5.EQ.67..AND.JJ.EQ.3)WRITE  BINARY (12) I3ARRAY 
IF(N5.EQ.66. .AND. JJ.EQ.4)WRITE  BINARY (12) I2ARRAY 
IF(N5.BQ.67. .AND.JJ.EQ.4)WRITE  BINARY (12) IlARRAY 
13  CONTINUE 


Put  a  border  and  title  at  the  bottom  of  the  picture 
DO  15  JF»1,3 

IF  (SHORT)  WRITE  BINARY  (12)  NC  /Space  right  twice 

IF  (SHORT)  WRITE  BINARY  (12)  NC  /for  short  picture 

DO  14  JOlfNl 

14  WRITE  BINARY (12) IL  /Print  a  line 

15  WRITE  BINARY  (12)  IFPC  /Terminate  the  line 

Write  BINARY (12) If  /End  with  a  line  feed 


nan  nnnnn 


16 


WRITE(12,16)ATnLE(l)  /Title  picture 
F0RMAT(15X,S80) 

CALL  RESET  ; Close  all  channels 

STOP 

Format  statements. 

16  FORMATC  '  ,15X, 'Signal  Processing  Laboratory,  Air 

1  Institute  of  Technology,  Wright-Patterson  AFB,  OH 

17  FORMAT (S13)  filename 

18  FQRMAT('O')  ; Double  space  for 

19  FORMAT(Sl)  ; Query  format 

20  ACCEPT"  Input  error.  Try  again  (YES/NO)  >  ■ 

GO  TO  1 

21  ACCEPT"  Input  error.  Try  again  (YES/NO)  >  " 

GO  TO  4 

Set  up  parameters  for  shortened  picture. 

22  TYPE”  Hie  shortened  (noise  removed)  picture  was  chosen." 

220  ACCEPT*Enter  starting  row  (1-256) :  ",JSTART 

IF(JSEART.LT.1.0R.JSTART.GT.256)GD  TO  220 
ACCEPT"  Number  of  lines  to  be  displayed  (255  Max)?  ",N2 
N2=MIN(N2,  (257-JSTART) ) 

SHOFa^.lSUE. 

NL=63 

IF(N2.GT.232) N2=228 
N3=3 
N4=65 
N5=66 
N6*4 

IF(JSTAKT.EQ.1)G0  TO  5 
DO  23  J  =  1,(JSTART-1) 

23  READ  BINARY (1) IREC 
GO  TO  5 
END 

C 

C*********  program  PICTURE  ***************************************** 


;Turn  on  short  test 
;Top  and  bottan  border  length 
;Number  of  lines  displayed 
/Location  of  left  border 
/Location  of  right  border 
/Location  of  line  feed 
/Length  of  lines  displayed 

/Skip  first  (JSTART-1) 

/Video  lines 


Foroe 

45433<14>') 

format 

short  picture 


nnn  r>  nnnnnnnnnnnn  nnn  nnnn 


SUBROUTINE  OUT3X3  (VIDPIX, UNEPRIOT) 
£******************************************************************** 

Written  by  Lt.  J.H.  Cromer 

This  subroutine  converts  the  video  pixel  values 
(i.e.  an  integer  value  from  0-15)  to  linepr inter 
dot  matrix  form.  A  3x3  array  pattern 
is  formed  by  this  subroutine.  Dot  pattern  texture 
(distribution  of  dots)  and  average  brightness  are 
varied  to  create  16  pseudo-gray  levels.  Odd  nunbered 
rows  of  the  picture  created  by  PICTURE  use 
these  3x3  patterns. 

NOTE:  The  six  least  significant  bits  of  each  byte 
sent  to  the  P-300  represent  print  hammer 
switches  (i.e.  a  1  turns  the  hartmer  on 
to  print  a  dot,  a  0  leaves  it  off) .  Bit 
seven  oust  have  a  value  of  1. 

(See  the  Printronix  manual  for  further  discussion) 

******************************************************************* 
INTEGER  VIDPIX, IJNEPRINT (3)  ,RIGHT,PRTTERN(3,16,2) 

Note  that  right  and  left  pixel  patterns  are 
not  necessarily  the  same. 


DATA  EATTERty/4*7,5,7,6,7,3,7,5,2*2,7,2,5,2,5,2,7 ,2*2, 

$5, 2, 5,0, 5,2, 5,2, 5,0, 2*2,0 ,2*2, 0,2, 1,3*0, 2, 4*0, 

$4*170K,150K,170K,160K,170K,130K,120K,150R,170K,4*150K, 

$120K,150K,120K,170K,2*120K,150K,120K,150K,100K,150K, 

$140K,120K,110K,120K,100K,150K,110K,140K,110K,100K,150K, 

$2*100K,120K,7*100K/ 

Riany  (  vidpix.  and.  15) +1 
LEFT=(ISHfT(VIDPIX,-4)  ,AND.15)+1 
DO  10  1-1,3 

LINEPRINT  ( I )  -PATTERN  ( I ,  UEFT  ,  1 ) +BATTERN  ( I ,  RIGHT ,  2 ) 

10  OOKTINUE 
RETURN 
END 
C 

q*********  subroutine  0UT3X3  ************************************ 


o  o  o'o  o  o  u  o  o  o 


SUBROUTINE  CUT4X3  (VIDPIX,LINEPRINr) 

C******************************************************************** 

Written  by  Lt.  Cromer 

This  subroutine  returns  linepr  inter  pixels  to  the 
calling  program  PICTURE,  which  sends  video 
pixels  (an  integer  from  0-15) .  The  pixel  pattern 
returned  is  a  4x3  dot  matrix  array,  to  be  used  for 
the  even  rows  of  pictures  created  by  PICTURE. 

(See  0UT3X3.FR  for  more  explanation) . 

C******************************************************************* 

INTEGER  VIDPIX,LINEPRINT(4)  ,RIGHT,BAT1ERN(4,16,2) 

DATA  EATIER1^5*7,5,7,7,6,7,7,3,3,6,7,5,5,6,3,2,5,2, 

$3*5, 2, 2, 5, 5,0 ,2, 5,2, 2, 5, 2, 2,1, 4, 2, 2,4 ,1,2, 1,4 ,2,0 ,0,1, 

$4,0,0,2,10*0, 5*170K,150K,2*170K,160K,2*170K,2*130K, 

$160K,170K,150K,120K,2*170K,120K,150K,120K,3*150R,2*120K, 

$150K,120K,2*150K,3*120K,150K,2*120K,110K,140K,120K,150K, 

$2*100K,150K,120K,100K,110K,140K,120K,2*100K, 

$120K,2*100K,140K,3*100K,120K,5*100K/ 

RIGHTS  (VH3PIX.AI03.15)  +1 
LEFT=(ISHFT(VIDPIX,-4)  .AND.15)+1 
DO  10  1*1,4 

LINEPRIMr(I)=EATTERN(I,LEPr,l)+PATTERN(I,RIGHT,2) 

10  CONTINUE 
RETURN 
END 

C 

C*********  subroutine  0OT4X3  ************************************ 


c******************************************************************** 

c 

C  Program  REMOVE  by  Lt  Jim  Cromer 

C  Fortran  5 

C. 

C  This  program  utilizes  a  3x3  pixel  mask  function 

C  to  remove  noise  from  packed  video  files.  Hie  mask 

C  can  be  changed  by  modifying  the  Subroutine  TEST3. 

C 

C  Execution  Line  Format: 

C  REMOVE  infile  outfilel  outfile2 

C  All  files  are  64  block  video  files. 

C 

C  Load  Line  Format: 

C  RLDR  REMOVE  TEST3  UNPACK2  IOF  @FLIB@ 

C 

C*******************************-**************************************** 

INTEGER  A(256)  ,B(256)  ,C(256)  ,D(256)  ,VH3BO(256)  rFILEl(7)  , 
$FILE2(7)  ,FILE3(7)  ,THRESH(4)  ,aOUNT(16,4)  ,0UT2(256)  ,DUM,MAIN(7)  , 
$0011(256) 

COMMON  /G0M2/  A,B,C,D,VIDEO 
COMMON  /OOM3/  THRESH,  COUNT, OJTl,OUT2 
C 

C  I/O  FILE  MANAGEMENT 

C 

CALL  IOF  (3  ,MAIN,FILE1  ,FILE2  ,FILE3  ,DUM,DUM,DUM,DUM,DUM) 

CALL  OP*M(l,FILEl,2,IER) 

DELETE  FILE2 
DELETE  FILE3 
CALL  CFII^(FILE2,2,IER) 

CALL  OPEN(2,FILE2,2,IER) 

CALL  CFILW(FILE3,2,IER) 

CALL  OPEN(3,FILE3,2,IER) 

C 

C  Enter  4  threshold  values 

C 

DO  101  1-1,4 

ACCEPT-ENTER  THRESH:  ",  THRESH  (I) 

TOTE*THRESH(",I,")-",THRESH(I) 

101  CONTINUE 
C 

C  Set  top  border  row  values 

C 

GALL  FDBLK(1,0,VIDEO,1,IER) 

IF(IER.NE.1)STQP  "  RDBLK  tO  error:", IER 
DO  1  1-1,64 

OOT2(I)»VIDEO(I) 

1  OUTl(I)-VIDEO(I) 

CALL  UNRACK2(1) 

CALL  UNPACK2(2) 

C 

C  Set  pixel  value  change  counters 

C 

DO  10  1-1,16 
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ono  non  o  noooo 


DO  10  J-1,4 

COUNT  (I,J)=0 

10  CONTINUE 

Operate  on  64  blocks.  Each  call  to  TESF3 
operates  on  3  rows  to  output  1  row.  Four  calls  to 
TEST3  are  made  for  each  packed  block  read. 

DO  2  1=1,63 
CALL  TEST3  (A,B,C,2) 

CALL  TEST3  (B,C,D,3) 

CALL  REBLK(1,I, VIDEO, 1, IER) 

IF(IER.IE.l)SIOP  *  RDBLK  error  **,I,"  error:", IER 
CALL  UNEACK2(1) 

CALL  TEST3(C,D,A,4) 

CALL  WRBLK (2,1-1, CUT1, 1, IER) 

IF(IER.IE.1)SP0P  "  WRBLK  #",I-1,"  error:", IER 
CALL  WRBLK(3,I-l,OUT2,l,IER) 

CALL  TESI3  (D,A,B,1) 

CALL  UNPACK2(2) 

2  TYPE"  <l5>Block  #",I,"  tested  .  .  ." 

CALL  TEST3(A,B,C,2) 

CALL  TEST3(B,C,D,3) 

Set  bottom  border  row  values 

DO  3  1=193,256 
OOTl(I)-VIDBO(I) 

3  OOT2(I)=VID0O(I) 

CALL  W®LK(2,63,0UT1,1,IER) 

IF ( IER. IE.  1) TYPE"  WRBLK  «63  error:”, IER 

CALL  NRBLX(3,63,OUI2,l,IER) 

Write  out  #  of  pixels  changed 

DO  20  J-1,3,2 
DO  20  1=0,15 

WRITE (12, 1000)1+1, 1, OOUNT (I, J)  ,1, 1+1, COUNT  (I+1,J+1) 
1000  FORMAT*"  #",I2,"S  reduoed  to", 13,":", 15,"  #",I2, 
$"s  increased  to  *,I2,":",I5) 

20  CONTINUE 

TYPE"  <15><7><7>Program  REMOVE  finished. <7>" 

CALL  RESET 

STOP 

EM) 


C 

C* 


Program  REMOVE 


SUBROUTINE  TEST3 (ROWl,RCW2,ROW3,I) 

C************************************************************** 

C 

C  (Called  by  REMOVE)  by  Lt  Jim  Cromer 

C. 

C  This  subroutine  determines  the  noise  removed 

C  output  value  of  an  input  pixel  by  calculating 

C  some  function  of  its  8  nearest  neighbors,  and 

C  comparing  it  with  a  user  input  threshold.  Two 

C  methods  of  noise  removal  are  used;  the  output 

C  arrays  are  OOT1  and  0UT2. 

C 

C  Backed  video  blocks  hold  4  rows.  The  parameter  I 

C  determines  which  part  of  the  block  that  CUT 

C  will  be  packed  into. 


c 

1-1  — > 

pack  OUT  into  elements  1-64 

c 

1-2  — > 

65-128 

c 

1=3  — > 

129-192 

c 

1=4  — > 

192-256 

C 

C***************************************************************** 

INTEGER  RGW1(256)  ,FCW2(256)  ,RCW3(256)  , TEMPI  (256) 

INTEGER  TEMP2(256)  ,0UT1(256)  ,0UT2(256)  ,THRESH(4) 

INTEGER  COUNT (16 ,4)  , COLUMN, SURROJtC 
LOGICAL  TEST 

COMMON  /C0M3/  THRESH,  OOUNTrOUTl,OUT2 
C 

C  Set  border  values 

C 

TEMPI (1)-B0W2(1) 

TEMP2 (1) =RQW2 (1) 

TEMPI (256) =R0W2 (256) 

THff2(256)«R0W2(256) 

C 

C  Check  each  of  the  neighbors  for  a  value  of 

C  the  interior  pixel  +/-1. 

C 

DO  20  CQLUMN-2,255 
N-l 

TEMPI  (COLUMN)  -R0W2  (COLUMN) 

TEST-R0W2  (OOLUMN) . EQ.O. OR. RCW2  (COLUMN)  .EQ.15 
IF  (TEST)  GO  TO  20 
I2-R0W2 (COLUMN) -3 
1  J«0 

12*12+2 

DO  5  M6  (COLUMN-1)  ,  (OOLUMN+1) 

IF  (R0W1  (M)  .EQ.I2)  J-J+l 
IF  (RGW2  (M)  .EQ.  12)  J-J+l 
IF (RCW3 (M)  .EQ.I2)  J*J+1 

5  CONTINUE 

TEST-J.GE.THRESH (N) 

IF  (TEST)TBtPl  (COLUMN)  -12 
IF  (TEST)  COUNT  (12  ,N)<XXJNT(I2  ,N)  +1 
IF  (TEST)  GO  TO  20 
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N»Nfl 

IF(N.BQ.2)Q0  TO  1 
20  CONTINUE 

Corapare  the  average  value  of  the  surround  less 
the  interior  value  with  a  threshold. 

DO  100  COLUMN  ,255 
TEMP2  (OOUMN)  =R0W2  (OOLIMN) 

TEST-BDW2  (OOLUMN)  .EQ.O . OR.  RGW2  (COLUMN)  .BQ.15 
IF  (TEST)  00  TO  100 
SURR0UN>-1*R0W2  (OOLIMN) 

DO  50  J*(00UMN-1) ,  (OQLIMNH) 

SURROUlt>=SURROlJlI>*-ROWl  (OOLIMN)  +BOW2  (OOLIMN) 

$+HOW3  (OOLIMN) 

0  CONTINUE 

AVERAGE>=FLQAT(  SURROUND) /8.0 
DIFF-AVERAGE-FLCAT  (RCW2  (OOLIMN) ) 

I2=ROW2  (OOLIMN)  +1 
I3=mr  (AVERAGE+0 .5) 

14=  (RCW2  (OOLIMN)  +13)  /2 
I5=MRX(I4,I2) 

I6=MIN(I4fI2-2) 

TEST-DIFF.GE.THRESH(3) 

IF  (TEST)  TEMP2  (OOLIMN)  =15 
IF  (TEST)  COUNT  (15 ,3)  KXXJOT (15  f  3)  +1 
IF (TEST) GO  TO  100 
TEST=  (-1*DIFF)  ,GB.THRESH(4) 

IF  (TEST)TEMP2  (OOLIMN)  =16 
IF  (TEST)  OOUNT  (12-2 ,4)  KXXJNT  (12-2 ,4) +1 
00  CONTINUE 

Kick  the  noise  removed  pixels  to  CtTTl,  OOT2 
M— 3 

JL0W= ( 1-1 ) *64+1 
JHIGK=JL0W+63 
DO  200  J*JD3W,JHIGH 
M=M+4 

OUT1  ( J)  =ISHET  (TEMPI  (M)  f  12)  +ISHFT  (TEMPI  (M+l)  ,8)  + 

$  ISHFT  (TEMPI  (M+2)  ,4)+TEMPl  (M+3) 

OUT2(J)=ISHFT(TEMP2(M)  ,12)+ISHFT(TEMP2(M+1)  ,8)  + 

$  ISHFT (TEMP2  (M+2)  ,4)-HEEMP2(M+3) 

200  CONTINUE 
RETORN 
END 
C 

€«#******  subroutine  'lEbT3  **************************************** 
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SUBROUTINE  UNPACK2(I) 

q*************************************************************** 

c 

C  (Called  by  REMOVE)  by  Lt.  Jim  Craner 

C. 

C  If  I»1  — >  The  first  two  rows  of  a  packed 

C  block  are  unpacked  into  arrays 

C  A  and  B. 

C  Else  — >  The  last  two  rows  are  unpacked 

C  into  arrays  C  and  D. 

C 

C  (See  the  unpacking  subroutines  for 

C  further  explanation). 

C 

C************************************************************** 

INTEGER  A(256)  rB(256)  ,C(256)  ,D(256)  ,VIDEO(256) 

COMMON  /OOM2/A,B,C,D, VIDEO 
C 

IF(I.EQ.1)G0  TO  3 
C 

L=-3 

DO  2  K=129,192 
L=D+4 

C  (L)  =15 .AND.  ISHFT  (VIDEO (K)  r-12) 

C  (1+1)  “15  .AND.  ISHFT  (VIDEO  (K)  ,-8) 

C  (L+2)  =15  .AND.  ISHPT  (VIDEO  (K)  ,-4) 
C(I+3)“15.AND.VIDEO(K) 

D  (L)  =15  .AND.  ISHPT  (VIDEO  (K+64)  ,-12) 

D  (L+l)  *15 . AND. ISHFT  (VIDEO (K+64)  ,-8) 

D(Df2) =15. AND.  ISHPT  (VIDEO (K+64) ,-4) 

D(W-3)  =15. AND. VIDEO  (K+64) 

2  OONTINCE 
GO  TO  6 

C 

3  D—3 

DO  5  K=l,64 
IrtH-4 

A(L)«15.AND.ISHFT(VIDBO(K)  ,-12) 

A(L+1)«15.AM).  ISHFT  (VIDEO (K)  ,-8) 

A  (L+2)  -15  .AND.  ISHFT  (VIDEO  (K)  ,-4) 

A  (L+3)  *15  .AND.  VIDEO  (K) 

B(L)  *15. AND.  ISHFT  (VIDEO  (K+64)  ,-12) 

B(Dfl)-15.AM>.  ISHPT  (VIDEO (K+64)  ,-8) 

B(L+2)  *15. AMD.  ISHPT  (VIDEO (K+64)  ,-4) 

B(L+3)  *15.AM).VIDEO(K+64) 

5  CONTINUE 

6  RETORN 


Subroutine  UNPA0K2  *********************************** 
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c ********************************************************************* 

c 

C  Program  EVIDHIST  Written  by  Lt.  Jim  Cromer 

C  Fortran  5 

C. 

C  This  program  calculates  a  histogram  of  the  16  gray  levels 

C  of  the  input  video  picture,  and  types  the  results  on  the 

C  CRT  and/or  linepr inter  in  tabular  form.  A  plot  can  be  sent 

C  to  the  linepr inter  if  requested. 

C 

C  Execution  Line  Format:  (run  VIM  1ST  on  the  NOVA) 

C  EVIDHIST 

C 

C  Load  Line  Format: 

C  RLDR  EVIDHIST  UNPACK  HISTHOT  HDT5.IB  @FLIB@ 

C  A  link  to  PL0T5.IB  in  DP5F:CALC0M5  should  exist  or 

C  be  created  before  loading. 

C 

C ********************************************************************* 

REAL  VALUE (0:15) 

INTEGER  INFILE  (7)  ,IARRAY  (2048)  ,IUNPACK(8192) 

C 

C********  I/O  FILE  MANAGEMENT  **************************************** 

c 

1  ACCEPT"Enter  name  of  video  file  to  be  evaluated 
$  — >  ■ 

READ (11 ,1000) INFILE (1) 

1000  P0RMAT(S13) 

CALL  OPEN(l, INFILE, 1,IER) 

IF (IER.NE.l) STOP"INFILE  OPEN  error  #",IER 
C 

£********  initialize  eata  ******************************************** 

c 

DO  10  1=0,15 

VALUE(I)=0.0 
10  00NTINIE 
C 

q********  mOCESS  TOE  PICTURE  **************************************** 

C 

KK=8  ;This  variable  determines  the  #  of  blocks 

;read  at  a  time.  (Maxinuro*16 ,  if  array 
; sizes  are  increased.) 

LL=64-KK 

DO  100  J=0,LL,KK 

CALL  BDBLK(1,J,IARRAY,KK,IER) 

IF(IER.JE.1)TYPE"1RDBLK  error:", IER 

ISIZE=256*KK  ;256  words/block 

CALL  UNPACK  (ISIZE,IARRAY,IUNPACK) 

IMAX=ISIZE*4  ;4  pixels/word 

DO  50  I=1,IMAX 

IPIX”IUNEACK (I) 

VALUE  (IPIX)  -VALUE  (IPIX)  +1 .0 
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C******  Print  the  histogram  ********* 

C 

100  CONTINUE 
TYPE"<15>" 

TYPE "Where  should  the  histogram  table  be  printed?" 

TYPE "Enter  1  to  send  to  linepr inter,  2  to  send  to  CRT," 

TYFE"or  3  for  both;  any  other  integer  to  continue." 
ACCEPT"<ll><ll>Enter  integer  option  — >  ",IOPT 
TYPE"<15>" 

IF(I0PT.LT.1.0R.  I0ET.GT.3)Q0  TO  160 
IF  (IOPT.EQ.l  .OR.  IOPT.EQ.3)  ICH=12 
110  IF(IOPT.EQ.2)ICH=10 

IF (ICH.EQ.10) WRITE (10,5000) INFILE (1) 

IF (ICH.EQ.12)  WRITE (ICH, 2000) INFILE (1) 

DO  120  J=0,7 
JJ=J+8 

WRITE  (ICH, 3000)  J, VALLE  (J) ,  JJ,  VALUE  (JJ) 

120  CONTINUE 

IF (ICH.EQ.12) WRITE (12, 4000) 

IF(IOPT.EQ.3)  IOPP=2 
IF (IOPT.EQ.2. AND. ICH.EQ.12) GO  TO  110 
160  TYPE"<15>" 

C 

C ******  plot  the  histogram  ********** 

C 

TYPE "Enter  1  to  plot  histogram  on  linepr inter,  any  " 
ACCEPT"other  integer  to  continue:  ",I 
IF(I.NE.l)GO  TO  200 
CALL  HISTPLOT (VALUE,  INFILE) 

C 

C*****  EXIT  program  or  start  over  ********* 

C 

200  CALL  RESET 

2000  FORMAT (///////27X, "  <10>HISTOGRAM  OF  ",S13,/25X, 

$"  - - "//) 

3000  PORMAT(18X,"  I  of  level ",I2,":\  15, 5X,"  #  of  level", 13,":", 15) 
4000  FORMAT ("0") 

5000  FORMAT (//27X,  "HISTOGRAM  OF  " , S13 ,/) 

TYFE"<15>" 

ACCEPT"Enter  1  to  evaluate  another  file,  any  other  integer 
$  to  stop:  ",I 
IF(I.EQ.1)G0  TO  1 
STOP 
END 
C 

C************  program  EVIDHIST  *************************************** 
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SUBROUTINE  HISTPDOT(VALU:  .  JJFILE) 

C***************************** *************************************** 

C 

C  Written  by  Lt.  Jim  Cramer 

C  This  subroutine  plots  the  16  element  array  VALUE 

C  passed  to  it  as  a  histogram.  It  is  called  by  the 

C  program  EVIDBIST. 

C 

c******************************************************************** 

REAL  VALUE (0:15), INFILE (7) 

REAL  X(481) rY(481) 

Y(1)=0.0 

K=1 

X(l)=*-0.25  ;X-axis  starting  point 

C 

C********  create  the  arrays  to  be  plotted  *************************** 
C 

VNORM=655.36  ;#  pixels/100 

1  ACCEPT"Enter  histogram  level  to  be  suppressed 
$  (999  to  continue):  ",ISUPH*ESS 
IF(ISUPPRESS.EQ.999)G0  TO  3 

IF  (ISUPPRESS.LT. 0  .OR.  ISUPFRESS. CT. 15) TYIE’INPUT  ERROR  1  <7><15>" 
IF  (ISUPERESS.LT.O  .OR.  ISUPFRESS. CT.  15)  GO  TO  1 
VNORM=VNDRM- VALUE  ( ISUPFRESS)  /100.0 
VALUE (ISUPFRESS) =0 .0 
GO  TO  1 

3  DO  20  1=0,15  ;do  16  values 

DO  5  J*l,15  ;determines  width  of  bars 

K=K+1 

X(K)=X(K-1) +1.0/30.0 
Y(K)=VALUE(I)/VN0FM 
5  CONTINUE 

X  (K-15)  =X  (K-14) 

DO  10  J=16,30  ;determines  spacing  between  bars 
K=K+1 

X(K) =X (K-l) +1.0/30 .0 
Y(K)=0.0 
10  CONTINUE 

20  CONTINUE 

C 

C  Find  end  of  input  filename,  then  insert  blanks  after  it 

C 

DO  30  J=2,13 

IF{  (BYTE  (INFILE,  J) )  .EQ.0)  BYTE  (INFILE,  J)  =32 
IF  ((BYTE  (INFILE,  (J-l) ) )  ,BQ.32)BYTE(INFILE,J)=32 

30  CONTINUE 

Y  (480)  =0.0 
X (480) =16.0 

XAX=5.0  ;set  X-  and  Y-axes  length 

YAX-4.0 
C 

C********  plot  the  arrays  ******************************************* 
C 

TYPE  "Please  wait  while  a  plot  is  generated  (50  secs)." 
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CALL  PL0TS(0,0,6) 

CALL  PL0T(1.5,5.0,-3) 

CALL  ASCALE(X,XAX,480,1,FX,DX)  ;  scale  the  arrays 
CALL  ASCALE  (Y,  YAX,480  /I  ,FY,DY) 

Title  the  axes 

CALL  AXIS (0.0 ,0.0, "PIXEL  VALUE", -11, XAX,0.0,FX,DK) 
CALL  AXIS (0.0, 0.0, "PERCENTAGE  OF  EVALUATED 
$  PIXELS", 30  ,YAX,90.0,FY,DY)  ;vertical  axis  title 

CALL  ALINE (X,Y,480,1,0,1,FX,DX,FY,DY) 


Title  the  plot 

CALL  SYMBOL (0.1, 4. 0,0. 2, "HISTOGRAM  OF  ",0.0,13) 

CALL  SYMBOL (3. 1,4. 0,0. 2, INFILE, 0.0, 13) 

CALL  PLOT  (0.0 ,0.0 ,999)  ;send  to  1  inept inter 
RETURN 
END 
C 

C**** *********  subroutine  HISTPLOT  ********************************** 
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q*************************************************************** 

Program  TONER  Written  by  Lt.  Jim  Cromer 

This  program  will  convert  individual  pixel  values 
of  an  input  video  file  into  new  values  assigned  by 
the  user.  It  can  be  used  to  adjust  gray-levels 
(i.e.  histogram  equalization),  increase  contrast, 
display  selected  pixel  values,  create  "negative  im¬ 
ages”,  create  files  of  a  constant  gray-level,  and  for 
many  other  purposes.  (NOTE:  The  program  EVIDHIST  can 
be  run  to  generate  a  histogram  of  any  VIDEO  file.) 

Execution  Line  Format: 

TONER 

The  program  will  ask  for  the  input  and  output 
file  names,  and  the  new  pixel  values. 

Load  Line  Format: 

RLDR  TONER  TIMER  UNPACK  REPACK  @FLIB@ 

C************************************************************** 

INTEGER  NEWVALUE  (0:15) , INFILE  (7)  ,0UTFILE(7) 

INTEGER  PACKED  (4096)  , UNPACKED  (16384) 

ISTART=0  ; start  timer 

ISTOP=1  ;stop  timer 


C 

£**«**★**  user  INPUT  OF  VARIABLES  ***************************** 

C 

ACCEPT  "Enter  name  of  video  file  to  be  toned  - >  " 

READ (11 ,1000) INFILE (1) 

ACCEPT"Enter  name  of  output  file  - >  " 

READ  (11, 1000)  OUTTTLE(l) 

TYPE"<15>" , "Enter  new  pixel  values<15>" 

TYPE”  NOTE:  Leading  zeros  are  significant,  i.e.  alter” 
TYPE"  a  '1'  as  '01',  enter  a  '2'  as  '02',  etc.<15>" 

DO  1  J=0,15  ;16  gray-levels 

WRITE (10, 2000)  J 
READ  (11 ,3000) NEWVALUE (J) 

IF (&SWVALUE ( J) .LT. 0. OR. NEWVALUE (J) .CT.15) 
$NEWVALUE(J)»15 
1  CONTINUE 

CALL  TBGR(ISTART) 

C 

C* *******  i/o  FILE  MANAGEMENT  ********************************* 

C 

CALL  OPEN (1, INFILE, 1,IER) 

IF(IER.IC.1)TYPE"INFILE  OPEN  error  #",IER 
CALL  DFILW(OUTFILE,IER) 

IF(IER.IE.l  •AND.IER.NE.13)TYI!E“CUTFILE  DFILW 
$  error  #",IER 

CALL  CFILW(0UTTTLE,3,64,IER)  /create  a  contiguous  file 
IF(IER.BQ.41)CALL  CFILW(0UTFILE,2,IER) 
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IF(IER.ie.l)TYIE"GUTFILE  CFILW  error  #\IER 
CALL  OPEN(2,OUTFILE,3,IER) 

IFdER.NE.DTYFE'OU'IFILE  OPEN  error  #",IER 
C 

0********  re-tq^E  tok  PICTURE  ******************************** 

C' 

DO  3  J=0 ,48,16 

CALL  REBLK(1,J,  PACKED,  16,  IER) 

IF(IER.M3.1)TYFE"RDBLK  #*,J,"  error:", IER 
CALL  UNPACK  (4096, PACKED, UNPACKED) 

DO  2  1=1,16384  jdo  1/4  of  picture 

UNPACKED (I)  =NEWVALUE  (UNPACKED (I) ) 

2  CONTINUE 

CALL  REPACK  (4096  ,UM>ACKED, PACKED) 

CALL  WRBLK(2,J, PACKED,  16, IER) 

IF(IER.NB.1)TYFE"WRBLK  #",J,"  error:", IER 

3  CONTINUE 

Send  message  to  CRT  terminal 

CALL  TDCR(ISTOP) 

WRITE(10,4000)OUTFILE  (1) 

TXPE"<15>" ,  "Have  a  nice  dayl<7><15>" 

C 

^* ********  write  NEW  TONE  VALUES  TO  THE  LINEPRINTER  *************** 

C 

WRITE (12,5000) 

WRITE  (12,6000)  INFILE (1)  ,0UTFILE(1) 

DO  5  1=0,15 

WRITE (12 ,7000) I ,NEWVALUE (I ) 

5  CONTINUE 
1000  F0RMAT(S13) 

2000  FORMAT {"  Change  old  pixel  value", 13,"  to  ?") 

3000  P0RMAT(I2) 

4000  FORMAT  ("  The  toned  picture  is  in  the  file  — >  ",S13) 

5000  FORMAT (//////26X , "  RESULTS  OF  TONER<10>"/ 

$26X,  "  - "////) 

6000  FORMAT (10X,"  Input  file  — >  ",S13,A0X,"  Output  file 

$ - >  " ,S13 ,//20X, "OLD  PIXEL",  10X,  "NEW  PIXEL",/20X, 

■  20X, "  "/) 

7000  F0PMAT(23X,  12 ,17X,  12) 

CALL  RESET 
STOP 
END 
C 

q***********  Program  TONER  ****************************************** 


Q********************************************************************* 

c 

C  Program  MOVE  written  by  Lt.  Jim  Cromer 

C  Fortran  IV  16  Aug  1982 

C. 

C  This  program  will  place  the  video  information  in  the  win- 

C  dow  given  for  the  template  (inputf  ilel)  inside  of  the 

C  window  given  for  the  background  (inputf  ile2) ,  and  write 

C  the  combined  picture  to  the  outputfile.  The  window  may  be 

C  placed  anyvhere  within  the  background,  and  may  be  taken 

C  from  anywhere  within  the  template.  Window  width,  length, 

C  and  position  are  input  by  the  user. 

C 

C  Execution  line  format:  (on  the  NOVA  only) 

C  MOVE  (run  EMJVE  on  the  ECLIPSE) 

C 

C  Loader  command  line  format  (NOVA  only) : 

C  KLER  MOVE  TEST  BLOCK  CHANGE  XWKBLK  XRDBLK 

C  UNPACK  REPACK  FORT.  LB 

C 

C********************************************************************* 

c 

INTEGER  IPAR(2)  ,INFILE1  (7)  ,INFH£2(7)  ,CUTFILE(7) , 

$CB1,  (BLOCKS,  CCDL,CIi3,CSTQP,CTOP,Q£,  CLEFT,  TTC^,TB1,TBL0CKS,CH3, 
$TOOL, HjS,TSTOP, TIB, HjEFT, 00MB (1024)  ,TOIP(1024)  , BACK (1024)  , WIDTH, TB 
C0MM3N/LIST1/G0MB,TEMP,CLS,TLS 
OOMMW/LIST2/LENGTH ,  WIDTH 
E9QUIVALQO  (COMB,  BACK) 

C 

C 

C*******  I/O  FILE  MANAGEMENT  *********************************** 

C 

TYFE"<15>" ,  "Program  MOVE  is  to  be  run  on  the  NOVA  only  I" 

99  TYPE"<15>" , "*******************************************<15>" 

ACCEPT"  Enter  template  file  name:  " 

READ(11,1000) INFILE1 (1) 

ACCEPT" <15>*,"  Enter  background  file  name:  " 

READ(11,1000) INFILE2 (1) 

DO  999  J»l,7 

999  OUTYILE  (J)  =INFILE2  ( J) 

ACCEPT" <15>","  Enter  combined  output  file  name:  ” 

READ (11, 1000) OUTFILE(l) 

1000  F0BMAT(S13) 

CALL  OPEN  (1 , INFILE1 ,2 ,  IER) 

IF  ( IER.  MI.  1}  TYPE"  Channel  1  OPEN  error:”, IER 
CALL  OPEN (2 , INFILE2 ,2 , IER) 

IF ( IER. MI.  1) TYPE"  Channel  2  OPEN  error:", IER 

CH3-2 

ICOUNT-O 

DO  1002  J-1,7  /check  for  BACRGR0UM>OGMBINED 

1002  IF  (OUTFILE ( J)  .EQ.INFILE2  ( J) )  IOOUNT-IGOUMPfl 

IF  (IC0UNT.BQ.7)  GO  TO  1 
CH3-3 

CALL  DFILW(OUTFILE,IER) 
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IF  (IER.EE.l  .AND.  IER.JE.13)TYFE"0UTFILE  DFILW  error  :",IER 
CALL  CFILW(0UTFILE,2,IER) 

IF(IER.NE.1)TOPE"  CFILW  error :",IER 
CALL  OPEN(CH3,OUTFILE,2,IER) 

IF(IER.IB.1)TOPE"  Channel  3  OPEN  error :",IER 


*******  enter  window  parameters  ******************************* 

1  ACCEPT"  <15>" , "  Enter  top  row  of  template  window  (1-256) 

ACCEPT”  Enter  left  oolum  of  template  window  (1-256)  :”,TLE1 
ACCEPT"  Enter  width  of  window  (1-256) WIDTH 
ACCEPT”  Enter  length  of  window  (1-256)  :*, LENGTH 

The  calls  to  TEST  check  to  see  if  the  input  parameters  are 
legal,  and  modifies  them  if  necessary: 

0  <  TOP  <  257,  (TOP  +  LENGTH)  <  258 

0  <  LEFT  <  257,  (LEET  +  WIDTH)  <  258 

CALL  TEST  (TTOP,TLEET) 

AOCEPT"<15>", "  Enter  top  row  of  background  window  (1-256) 
ACCEPT"  Enter  left  column  of  background  window  (1-256)  :",Q 
CALL  TEST  (CTOP, CLEFT) 

CALL  HUXX(TBL0CKS,TB1,TES,T0DL,TTOP,TLEET) 

CALL  BLOCK  (CHL0CKS,CB1 , CIS, COOL, CPOP, CLEFT) 


Determine  oolum  number  of  the  last  video  row  (0-3) 

J1=MQD(UNG1H,4) 

TSTOMCD(  (T00L+J1)  ,4)-l 
CSTOP=MOD{  (CCOL+Jl)  ,4)-l 
IF  (CSTOP.EQ.-l)  CSTOM 
IF  (TSTOP.EQ.-l)TSTOM 

Determine  the  last  significant  block  of  window 

CLBKB1+CBL0CKS-1 

TLB=TB1-«BL0CKS-1 

User  check  of  window  parameters 

TOPE  "a5>","WIEfIH«",  WIDTH,"  LENGTH-" ,  LENGTH 

TOPETEMLATE  TOP  RCW*",TTOP,"  BACKGROUND  TOP  RCW-",CTOP 
TOPE  "TEMPLATE  LEET  OOLCW-",TLS, "  BACNGROUM)  LEET  OOLIMN 
$»",OS,"<15>" 

ACCEPT "Enter  l  to  see  expanded  set  of  variables,  any 
$  other  integer  to  continue:  ",I 
IF(I.M2.1)QO  TO  5 

TYPE"************************************************" 

TOPE”  PARAEETER  TEMRATE  BACHGROUM)" 


nrtfloortrtftfio  o  on  n  non  noon 


TYPE"<15>", "  FIRST  BLOCK" ,1B1  ,CB1 
TYFE"<15>","  LAST  BLOCK  ",TU3,CX£ 

TYFE*<15>","  *  OF  BLOCKS", TBLOCKS,CBLXKS 
TYFE"<15>", "  LEFT  OOL  ",TLS,CLS 
TYPE"<15>","  WIDTH-  ",  WIDTH 
TYPE"<15>","  LaCTH-MiQKJffl 

TYPE"<15>", "***********************************************■ 
ACCEPT"  Enter  1  to  try  another  set,  any  other  integer 
$  to  continue:  ",I 
IF(I.BQ.1)G0  TO  1 


******  create  the  contained  picture  ************************** 


5  IO0UNT=0 

IF(CH3.EQ.2)G0  TO  20  ; If  combined  picture  file 

?is  the  same  as  the  back- 
; ground  picture  file,  then  no 
;need  to  write  to  itself 
Write  background  only  blocks  (before  window) 
to  the  combined  picture  file. 

JMAX=CB1-1 

IF(JMAX.LT.O)QO  TO  20 
DO  10  J=0,JMAX 

CALL  RDBLK(2,J,BACK,1,IER) 

IF(IER.1B.1)TYPE"  2FDBLK",J,"  error :",IER 
CALL  WRBLK  (CH3,J,  COMB,  1,  IER) 

IF(IER.NE.1)TYFE"  WRBLK",J,"  error: ",IER 
ICOUNMCOUNW-l 

10  OONTINUE 

20  TYPE"  Background  before  window  completed." 

TYPE"  #  Blocks  written: ",ICCUNT 


.  .  .  Overlay  template  window  onto  background . 

CALL  XRDBLK(1,TB1,  TEMP,  1,  IER) 

IF(IER.NE.1)TYFE"1KDBLK  #",0B1,"  error:", IER 
CALL  XRDBLK(2,CB1, BACK, 1, IER) 

IF(IER.NE.1)TYPE”2FDBLK  I*,CB1,"  error:", IER 
Nl-TOOL  ; 4-MAX (Nl, N2)  gives  the  nutter  of  rows 

N2-C00L  ;to  change  before  the  next  RDBLK 

IF(TOX.GT.CQX)GO  TO  100 


There  are  four  colums  in  the  packed  video  array  (64x4) , 
designated  0,  1,  2,  and  3.  If  the  template  starting 
column  nutter  is  less  than  or  equal  to  the  background  (combined) 
starting  column  nutter,  then  the  background  block  will  be  "used 
up"  before  the  tenplate  block.  When  the  background  block  is 
finished,  a  WRBLK  is  done,  and  the  next  background  block  is  read 
When  the  template  block  is  finished,  the  next  tenplate  block  is 
read,  but  no  WRBUC  needs  to  be  performed.  Note  that  the  back- 


ground  and  combined  files  are  always  at  the  same  block  number. 

CALL  CHANGE (N2 ,N2,NL) 

CALL  XWRBLK  (CH3  ,CB1 ,00MB,  1  ,  IER) 

IF(IER.NE.1)TYPE"  WRBLK  #",CB1,"  error:", IER 

Write  the  template  window  into  the  background 

TB=TB1+1 

IMDKB1+1 

ioounm. 

IMAX=OB-l 

IF (IMIN.GT. IMAX)QO  TO  60 
DO  50  I=IMIN,IMAX 

CALL  XRDBLK(2,I,BACK,1,IER) 

IF(IER.IE.1)TYPE"  2RDBLK  #",I,"  error:", IER 
CALL  CHANGE (N1,N2,N1) 

CALL  XRDBLK(1,TB,T01P,1,IER) 

IF ( IER.  NE . 1 ) TYPE  *  1RDBLK  #",TB,"  error:", IER 
CALL  CHANGE (N2,N2, Nl) 

CALL  XWRBLK(CH3, 1,OOMB,1,IER) 

IF ( IER. l®.  1) TYPE"  WRBLK  #",I,"  error:", IER 
ICOONMGOUOTtl 

TB*TB+1 

TYPE”  TOOL. LT. COX — Window  portion  complete. " 

TYPE"  #  blocks  written: ",IC0UNT 
GO  TO  250 

In  this  case  the  template  starting  column  number 
is  greater  than  the  background  starting  oolum  number. 

The  template  block  must  be  "finished”  first. 

CALL  CHANGE (N1,N2,NL)  ; finish  TEMP  block 

TB-TB1+1 

IMAX<LB-1 

IOOUNT-O 

IF  (CB1  .GT.  IMAX)  GO  TO  225 
DO  200  I=CBl,IMAX 

CALL  XRDBLK(1,TB,TEMP,1,IER) 

IF  (IER.  I®.  1}  TYPE"  1RDBLK  #*,1B,"  error:", IER 
CALL  CHANGE  (N2,N2,N1)  /finish  BACK  block 

CALL  XHRBLK(CH3,I,G0MB,1,IER) 

IF ( IER. I®.  1) TYPE"  WRBLK  *",I,"  error:"fIER 

ioouN>iaxjNm 

IBLK-I+1 

CALL  XRDBLK(2,IBLK, BACK, 1, IER) 

IF ( IER. I®. 1) TYPE"  2RDBLK  #",IBLK,"  error:", IER 
CALL  CHANGE (N1,N2,N1)  /finish  TEMP  block 

TB^B+1 

TYPE"  TOOL. (TT. COOL — Window  portion  complete." 

TYPE”  #  blocks  written: ",ICOU NT 


If  the  combined  (background)  stopping  oolum  nunber 
is  greater  than  the  template  stop  col  inn  nunber,  then  the 
second  last  template  block  (I2IB*TIE-1)  must  be  read  (i.e. 
there  are  more  video  rows  to  be  changed  in  the  last  back¬ 
ground  block  than  there  are  available  in  the  last  template 
block  to  change  them  to) .  If  TSTCP  is  greater  than  or  equal 
to  CSTCP,  then  there  are  sufficient  video  rows  available  in 
the  last  template  block  to  complete  the  last 
background  block  to  be  changed. 

IF  (CSTOP. CT. TSTCP)  00  TO  400 
M1=TSTOP-CSTOP 

CALL  XRDBUC(1, TIB, TEMP,  1,  IER) 

IF (IER.  183.1) TYPE"  1FDBLK  #",HB,"  errors", IER 
CALL  XRDBLK(2,CXB,BACK,1,IER) 

IF(IER.1E.1)TYPE"  2RDBLK  #",OB,"  error:",IER 

N1*3-CST0P 

N2=0 

CALL  CHANGE (N1,N2, Ml)  ; finish  BACK  block  to  CSTCP 
CALL  XWKBLK(CH3,CXB,OOMB,1,IER) 

IF(IER.NE.1)TYIE"  WRBLK  #",CLB,"  errors", IER 
TYPE"  CSTOP.LT. TSTOP— Last  block  of  window  complete." 

GO  TO  500 

Complete  the  last  block  of  the  window 
NOTEs  CSTCP  is  greater  than  TSTCP.  Therefore  finish 
TEMP  before  BACK. 

Ml=CSTOP-TSTOP 

CALL  XRDBLK(2,OB,BACK,l,IER) 

IF  (IER.  IE.  1)  TYPE"  2RDBLK  #",CXB,"  errors  ",IER 
I2IB=H£-1 

CALL  XRDBLK(1,I2LB, TEMP,  1, IER) 

IF ( IER. IE. 1) TYPE"  1RDBLK  #",I2IB,"  errors", IER 

N1=4-M1 

N2**0 

CALL  CHANGE (NL,N2,N1)  ; finish  TEMP  block 

CALL  XRDBUC (1, TIE, TEMP,  1, IER) 

IF(IER.1E.1)TYH:"  lRDBLK  #",TIB,"  errors", IER 
M1-3-TST0P 

CALL  CHANGE  (Ml  ,N2,N1)  ;  finish  BACK  block  to  CSTCP 

CALL  XWFBLK(CH3,CIB,OCMB,l,IER) 

IF ( IER, IE.  1) TYPE"  WRBLK  #",OB,"  errors", IER 

TYPE"  CSTOP.CT. TSTOP — Last  block  of  window  complete." 

IGOONM 

Finish  the  combined  file  (background  only  portion) 

JMHK2B+1 

IF(JMIN.OT.63)G0  TO  601  ;if  finished  ST CP 

IF (CH3.EQ.2)G0  TO  601  ;if  CX3MBINEI>>BACKGR0U1D,  STO 

DO  600  J-JMIN,63 

CALL  FCBLK(2, J, BACK, 1, IER) 


IF(IER.re.l)THE"  2REBLK  #",J,"  error:", IER 
CALL  WRBLK(CH3,J,OOMB,1,IER) 

IF  ( IER.  EE.  1)  TYPE"  HRBLK  *",J,"  error :",IER 
iGourayiGouMPfi 

600  OONTINLE 

TYPE"  Finished  background  only  portion." 

601  TYPE"  *  blocks  written*  MCOUNT 

TYPE"<15>",  "<7>"f  "<15>","  Program  NMOVE  execution  completed.  <7>" 
WRITE  (10, 2000)(XrTFILE(l) 

2000  FORMAT ("  The  combined  picture  is  in  the  file  — >  ",S13) 

C 

C 

C******  present  option  Msnu  ********************************* 

C 

GO  TO  2010 

2002  TYPE"<15>", "Input  error.  Try  again." 

2010  TYPE"<15>" , "******************************************■ 

TYPE " <15 >", "What  next?<15>","Here  are  the  options:" 

TYPE*<15>" ,  "<11>1  -  Try  another  set  of  window  values" 

TYPE" <15 >","<11  >2  -  Start  over  with  new  input  pictures” 

C  TYPE"<15>", "<11>3  -  Display  combined  picture  on  the  video  monitor" 

TYPE"<15>" ,  "<11>3  -  Save  combined  picture  and  STCP<15>" 

ACCEPT" <ll>Enter  option  — >  ",IOPT 
IF(IOPT.UP.1.0R.IOPT.<IP.3)GO  TO  2002 
IF(IOPT.EQ.1)GO  TO  1 
CALL  RESET 

IF(IOPT.EQ.2)GO  TO  99 
IF  (IOPT.BQ.3)  STOP 

TYPE"<15>" , "Check  monitor  —  Press  green  CHOPS  control 
$button  to  aontinue." 

IDCOT=4 

IPAR(l)-9999 

IEAR(2)-0 

WRITE  (10, 3000)(XnFILE(l) 

3000  FORMAT ("0", "Picture  being  displayed  — >  ",S13) 

C  CALL  CHANNEL  (0,0,3 ,0,0,  "A", 0,0,0,  IE,  IS)  ;call  abort 

C  CALL  CHANN&(3, 1,2,1, IDCNT,OUTFILE, 64,0, IPAR,IERR,ISYS) 

C  CALL  ERCHK(IERR,1  ,IDQW,1,ISYS) 

TYPE"<15>","CHANIEL  currently  not  loaded." 

TYPE  "Use  VIDEO  to  display  combined  pictures.  <15>" 

CALL  OPEN(l,INFILEl  ,2,IER)  ;re-OPEN  channels 

IF(IER.NE.1)TYPE"CH1  RE-OPEN  ERROR:", IER 
CALL  OFEN(2,INFILE2,2,IER) 

IF ( IER. IE .  1 ) TYPE "CH2  RE-OPEN  ERROR:", IER 
IF(CH3.EQ.3)GALL  OPEN(3  ,CUTFILE,2,IER) 

IF ( IER. NE . 1 ) TYPE "CH3  RE-OPEN  ERROR:”, IER 
GO  TO  2010 
END 
C 

C*********  Program  (MOVE  ******************************************** 


nnnnn 


SUBROUTINE  TEST  (TCP, LEFT) 

C******************************************************************* 

Subroutine  TEST  checks  to  see  if  the  input  parameters 
to  program  MOVE  are  legal,  and  modifies  them  if 
necessary.  (It  is  also  called  by  DISTANCE.) 

C******************************************************************* 

INTEGER  TOP, WIDTH 
OOMMDN/LIST2/LENGIH ,  WIDTH 
IF(LEET.LT.l  .QR.I£FT.Gr.256)  LEET=1 

MAXWIOTH=257 -LEFT  /picture  has  256  colunns 

if  (wnxra .  or .  maxwioth  .  or.  width .  lt  .  1 )  wmiH=MAXwiniH 

IF(T0P.LT.1.C«.TOP.GT.256)T0P=1 

MAXL£NGTH=257-TOP  /picture  has  256  rows 

IF  (LENoIH.CTr.MAXLENGTH.  OR.  LENGTH.  LT.1)LENGTH=MAXLE2«3TH 
RETORN 
END 
C 

C*********  subroutine  TEST  ***************************************** 


nnnnoo 


SUBROUTINE  BLOCK  (NllBIOCKS,BliXKl,  LEFTS  IDE,  COLUMN,  TOP,  LEFT) 
eft********************************************************************* 


Subroutine  BLOCK  determines  the  total  nunber  of  blocks  to 
be  read  into  the  window,  the  first  block  to  be  read, 
and  the  first  video  row  "column"  number.  This  subroutine 
is  called  by  WOVE. 

C********************************************************************** 

INTEGER  BLOCKl/COLLMN, TCP, REMAINDER, WIETH 
OMON/LIST2/LE1CTH ,  WIDTH 

BLOCK1=IOT((TQP-1)/4.0)  ;4  rows  per  block 

OQLCMWOD(  (TOP-1)  ,4) 

LEFTS  IDE«LEPr 
REMAINDE&40D  (LENGTH,  4) 

Kl=LENGTB+3 

NDMBLOCKS»INT(K1/4.0) 

IF  (RH1AINDER.EQ.2  .AND.OOLUMN.GT.2)  NDMBL0CKS=NUMBL0CKS+1 
IF  (REMAINDER. EQ.3  .AND.OOLUHN.GT.l)  NCWBLOOCS=NUMBLOCKS+1 
IF  (REMAINDER.  EQ.O  .AND.00L111N.GT.  0)  NCMBLOOCS-NUMBLOCKS+ 1 
IF  (NUMBLOCKS.CT.l)  RETURN 

TYPE "WARNING:  t  Blocks  to  be  read  =",NUMBLOCKS 
PAUSE 
RETURN 
END 
C 

C************  subroutine  BLOCK  ************************************* 


n  n  o.  o  o  o 


SUBROUTINE  CHANGE  (JMIN,  (START, TSTART) 

q******************************************************************* 


Written  by  Lt.  Jim  Cromer 

Subroutine  CHANGE  changes  the  corresponding  background 
(i.e.  the  combined  picture)  pixels  to  template  pixels; 
it  is  called  by  the  program  MOVE. 

C******************************************************************* 

INTEGER  OOMB  (1024)  ,T31P(1024)  rCLS  ,TLS ,  (START,  TSTART, WIDTH 
COMMON  /LIST2/  LENGTH,  WIOTH 
COMMON  /LIST1/  OOMB, TEMP, CLS,TLS 
C 

DO  2  J=JMIN,3 

K=TSTABT*256+TLS  ;Set  left  side  of  input (template 

M=CSTART*256+CLS  ;and  output  (combined)  windows. 

KMAX=K+WIEnn-l  ;  Change  values  over  the  width  of  window 
DO  1  L=K,KMAX 

COMB(M)=TEMP(L) 

1  M=M+1 
TSTART=TSTAR1M-1 

2  CSTART=CSTARTtl 

IF(CSEART.EQ.4)CSrART=0  ; reset  row-pointer  if  necessary 

IF(TSTART.EQ.4)TSTART=0 
RETURN 
END 
C 

C*********  subroutine  CHANGE  ***************************************** 


£******************************************************************** 

c 

C  Program  NEGATE  by  Lt.  Jim  Cromer 

C  Fortran  5 

C, 

C  This  program  writes  to  an  output  video  file  the  "negative" 

C  pixels  values  of  the  input  video  file  (i.e.  dark  pixels 

C  become  light,  and  vice-versa) .  It  will  also  "flip"  the 

C  picture  about  its  horizontal  axis  (i.e.  turn  it  upside  down) . 

C 

C  Execution  Line  Format: 

C  NEGATE[/F]  infile  outfile 

C  The  /F  is  chosen  if  the  picture  is  to  be  reversed 

C  horizontally.  The  program  can  be  run  twice  to  produce 

C  a  horizontally  flipped  positive  image. 

C 

C  Load  Line  Format: 

C  RLDR  NEGATE  IOF  UNPACK  REPACK  TIMER  @FLIB@ 

C 

C********************************************************************* 

INTEGER  INFILE (7) , OUTFILE  (7)  ,GS(2)  ,NEW(1024)  ,0LD(1024) 

INTEGER  VH3EO(256)  ,DUM,MAIN(7) 

C 

q********  I/O  FILE  MANAGEMENT  *************************************** 

c 

c 

CALL  IOF  (2, MAIN,  INFILE, OUTFILE,  DUM,GS,IXJMfEUM,IXJM) 

CALL  TlfGR(O)  ; start  timer 

OPEN  1, INFILE 

DELETE  OtTTFILE 

CALL  CFILW (OUTFILE, 3, 64, IER) 

IF  (IER.EQ.l)  TYPE  "Contiguous  file  created" 

IF  ( IER.  EQ.  41)  CALL  CFILW  (OUTFILE,  2,  IER) 

IF(IER.NE.l)STOP  "  Random  file  creation  error  " 

OPEN  2,  OUTFILE 
C 

C********  Test  switch  and  set  variables  **************************** 

C 

IF(GS(1) .EQ. 1024) GO  TO  1  ;test  global  switch 

TYPE"  Creating  negative  video  file  only." 

11=0 
12=256 
13=512 
14=768 
GO  TO  2 

1  TYPE"  Creating  negative  and  horizontally  flipped  video  file." 
11=768 
12=512 
13=256 
14-0 
C 

C********  Loop  around  the  next  section  64  times 

C  to  process  the  entire  picture  ********************* 

C 


noon 


2 


DO  4  1=0,63  jdo  entire  picture  (64  blocks) 

K=I 

IF(I1.BQ.768)K=63-I  ;if  flipped,  start  at  bottom  of  infile 
CALL  RDBLK(l,K,VXDEO,l,IER) 

IF(IER.JC.1)TYFE  ■  RDBLK  error  #s",IER  STOP 
CALL  UNPACK  (256  /VIDEO, CLD) 

Negate  pixels,  and  re-arrange  if  required. 

Work  on  4  rows  at  a  time. 

DO  3  N=l,256 
NEW(N)=15-0LD(N*I1) 

NEW  (Nf 256)  =15-0LD(NtI2) 

NEW(Nf512)  =15-OLD  (N4-I3) 

NEW (Nf 76 8)  =1 5-OLD  (Nf  14) 

3  CONTINUE 

CALL  REPACK (256, NEW, VIDEO) 

CALL  HRBLK(2,I,VIDEO,l,IER) 

IF(IER.NE.1)ST0P"  WRBLK  error  #",IER 

4  OONTINUE 
C 

C*******  Write  aompletion  message  to  CRT  ************ 

C 

CALL  RESET 
CALL  TIMER(IOOO) 

WRITE  (10,5)OUTFILE  (1) 

5  FORMAT ("  The  negative  picture  is  in  the  file  — >  ",S13) 

STOP 

END 

C 

q*************  program  NEGATE  ****************************************** 


c******************************************************************** 

c 

C  Program  REDUCE  Written  by  Lt.  Jim  Cromer 

C  FORTRAN  5 

C. 

C  This  program  reduces  a  256x256  pixel  video  picture  into 

C  a  128x128  array  by  averaging  4  pixels  into  1.  The 

C  reduced  array  is  placed  in  the  upper  left-hand  quadrant; 

C  all  other  pixels  in  the  output  picture  are  made  white  for 

C  display  purposes. 

C 

C  Execution  Line  Format: 

C  REDUCE  infile  outfile 

C  Both  the  input  and  output  files  are  256x256  video  files 

C  (i.e.  packed  integer  form,  64  blocks  per  file) . 

C 

C  Load  Line  Format: 

C  RLDR  REDUCE  IOF  TIMER  UNPACK  REPACK  @FLIB@ 

C 

c*********************************************************************** 
INTEGER  CLDPACKED(4096)  , OLDUNPACK (16384)  ,N0fUM>ACK(8192) 

INTEGER  NEWFACKED  (2048)  , INFILE (7)  , OUTFILE  (7)  ,  MAIN  (7) 

COMMON  QLDPACKED , OLDUNEACK 

EQUIVALENCE  (OLDUNPACK  ,  NEWUNPACK)  f  ( NEWFACKED, CLDPACKED ) 

C 

C 

Q*t*****+  Z./0  FILE  MANAGEMENT  ************************************** 

C 

CALL  IOF(2,MAIN,INFILE,OUTFILE,I1,I2,I3,I4,I5) 

CALL  TDCR(O)  ;start  timer 

CALL  DFILW  (OUTFILE  ,  IER) 

IF (IER. IE. 1. AND. LER.NE.135TYPE "OUTFILE  DFILW  error:", IER 
CALL  CFILW (OUTFILE, 3 ,64, IER) 

IF  (IER. EQ. 41)  CALL  CFILW (OUTFILE ,  2  ,  IER) 

IF  (IER.  IE.  1)  TYPE  "OUTFILE  CFILW  error:", IER 
CALL  OPEN  (1, INFILE,  1, IER) 

IF (IER.NE.l) TYPE" INFILE  OPEN  error:", IER 
CALL  OPEN (2, OUTFILE, 3, IER) 

IF  ( IER.  NB.l)  TYPE  "OUTFILE  OPEN  error:", IER 
C 

C********  process  the  picture  *************************************** 

C 

DO  400  M-0,3 

Ml«41*16  ;RDBLK  counter 

CALL  RDBLK(1, Ml, CLDPACKED, 16, IER)  ;read  64  rows 
IF(IER.IE.1)TYFE"RDBLK  #",M1,"  error:", IER 
CALL  UNPACK  (4096 ,0LDPACKED, OLDUNPACK) 

K*0  ;new  row  counter 

C 

C  This  section  reduoes  64  rows  of  256  elements  each  into 

C  32  rows  of  128  elements  each.  It  executes  this  function 

C  4  times  so  that  a  total  of  256  rows  are  reduced  to  128. 

C 

DO  300  J-1,32  /create  32  new  rows  from  64 


112 
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DO  100  1=1,128  /create  128  eleroents/new  row 
K-K+l 

L=K*2  /old  row  1  counter 

LL=L+256  /old  row  2  counter 

A1*CLDUNPACK(I/-1) 

A2=CLDUNPACK(L) 

A3=CLDUNPACK(LL-1) 

A40LDUNPACK  (LL)  /average  the  4  pixels 

NEWUNPftCK (K) =IFIX ( (Al+A2+A3+A4)/4.0+0.5) 

100  CONTINUE 

IMIN=K+1 

IMAX=K+128 

DO  200  I=IMIN, IMAX  /finish  video  row 

200  NEWUNRACK ( I )  =15  y set  pixels  to  white 

300  K=K+128  /jump  to  next  row 

CALL  REPACK  (2048,NEWUNPACK,NEWPACKED) 

M2=m/2  /write  one-half  the  tof  blocks  read 

CALL  WRBLK(2,M2,NEWPACKED,8,IER) 

400  IF(IER.NE.l)TyPE"WRBIK  #",M2,"  error:", IER 

C 

q*********************************************************************** 

Set  the  rest  of  the  picture  to  white  (15)  for  the  most 
aesthetic  display.  These  pixels  are  later  set  to  zero 
by  the  program  NORMALIZE. 

DO  500  Jsl,4096 
500  CLDPACKED ( J) =177777K 

DO  600  J=32 ,48,16 

CALL  WRBLK (2, J, CLDPACKED,  16 ,IER) 

600  IF(IER.M;.l)T3fPE"WRBLK  #",J,"  error:", IER 

CALL  TIIGR(I)  /stop  timer 

WRITE  (10 ,700)  OOTFILE  (1) 

700  FORMAT ("  The  reduced  picture  is  in  the  file  — >  ",S13) 

CALL  RESET 
STOP 
END 
C 

C*************  Program  REDUCE  **************************************** 


oonnooooooooonnnnnnonnorionnooooooonooooorjnoo 


C******************************************************************** 

Program  NORMALIZE  Written  by  Lt.  Jim  Cromer 

Fortran  5 

If  switch  'U*  is  chosen: 

This  program  normalizes  the  upper  left  quadrant  of  an 
input  packed  video  file  into  twenty-four  (30  oolunn  by 
20  row)  normalized  grid  blocks,  with  the  average  sun 
of  the  squares  of  the  normalized  pixel  values  per  unit 
template  window  area  equal  to  one  (assumes  a  23x47 
reduced  template  window) . 

If  switch  'L'  is  chosen: 

The  entire  lower  right  quadrant  is  normalized  to  give 
it  an  energy  of  unity.  Switch  L  is  usually  used  to 
normalize  templates. 

Otherwise: 

The  program  will  ask  the  user  to  input  the  nunber 
of  horizontal  and  vertical  grid  rectangles. 

Choose  from: 

Horizontal  — )  1,2,3,4,3,6,8,10, 

12,15,20,24,30,40,60,120 
Vertical  - >  1,2,3,5,6,10,15,30 

The  output  file  is  a  256x256  element 
complex  contiguous  file  (or  random,  if  a  contiguous  file 
cannot  be  created) .  The  program  assumes  the  input  file 
is  a  reduced  picture  in  the  upper  left  or  lower  right  hand 
quadrants.  It  is  normally  used  in  the  sequence 

->  REDUCE  NORMALIZE  DIRECT  CttJUTIILY  INVERSE  <- 

Execution  Line  Format: 

NORMALIZE [/U  or  L]  infile  outfile 
One  of  the  switches  must  be  selected  for  use  as  an 
automatic  program,  as  in  a  macro  file;  a  U  indicates 
to  normalize  the  upper  left-hand  quadrant;  an  L 
indicates  that  the  lower  right-hand  quadrant  is  to 
be  normalized.  Pixel  values  outside  of  the  selected 
quadrant  will  be  set  to  zero. 

Load  Line  Format: 

RL ER  NORMALIZE  IOF  TIMER  UNPACK  @FL3B@ 

C********************************************************************** 

REAL  NORM (120, 30) , ENERGY (120,30) 

INTEGER  ARRAY (256) , INFILE (7) , OUTFILE (7)  ,R0W(32,4,2) 

INTEGER  ONEOOUNT,TNOOOUNT,MS(2) , UNPACKED (51 2)  ,MAIN(7) 

INTEGER  VGRID,HGRID,GRDBLK 
INTEGER  WIDTH ,  VNUM0FG,  HNUMOPG,  BEERGRID 
COMPLEX  CNORM(1024)  ,CZERO 
LOGICAL  TEST 

EQUIVALENCE  (NORM, ENERGY) 
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ISEARP-0 

ISPOE^l 

C 

C 

0 ******  i/Q  FILE  management  ****************************************** 

c- 

CALL  RESET 

CALL  IOF(2, MAIN, INFILE, OUTFILE, IDUM, MS, 12, 13, 14) 

CALL  DFILW  ( "TEMP" ,  IER) 

IF(IEK.lE.l.AND.IER.ie.l3)TYIE"TEMP  DFILW  error:", IER 
CALL  CFILW  ( "TEMP" ,  2  ,  IER) 

IF  {IER.  NB.  I)  TYPE  "TEMP  CFILW  errors’1, IER 

IER=1 

JER=1 

CALL  OPEN  (2, OUTFILE,  2,  IER) 

IF ( IER. EQ. 1)00  TO  55 

TYPE "Attempt ing  to  create  a  contiguous  file." 

IF(IER.EQ.13)CALL  CFILW  (OUTFILE,  3 ,1024,  JER) 

IF  ( JER.  41)  TYPE  "Successfully  created  a  contiguous  file." 

IF  (JER.  BQ.  41)  TYPE  "Must  create  a  random  file  instead." 

IF  (JER.  EQ.  41)  CALL  CFILW  (OUTFILE,  2,  JER) 

IF  ( JER.  IC.l)  TYPE  "OUTFILE  CFILW  error:",  IER 
IF  (IER.  EQ.  13)  CALL  OPEN  (2, OUTFILE,  3,  IER) 

IF  ( IER.  IG.l)  TYPE  "OUTFILE  OPEN  ERROR:", IER 
55  CALL  OPEN (0, INFILE,  1, IER) 

IF  ( IER.  NS.  1)  TYPE  "INFILE  OPEN  error:*, IER  STOP 
CALL  OPEN (1, "TEMP", 2, IER) 

IF  (IER.  NS.  1)  TYPE  "TEMP  OPEN  error:", IER 
C 
C 

0******  determine  switches  and  set  variables  ****************** 

c 

IF  (MS  (1)  .E)Q.16)G0  TO  500  yswitch  was  L 

IF  (MS  (2)  .E3Q.2048)GO  TO  1  yswitch  was  U 

GO  TO  556 

5555  TYPE"Input  errorl<7>  0<  input  <121","<15>" 

556  ACGBPT”Enter  horizontal  #  of  grid  rectangles 
$  (1-120):  ",HNUMDFG 

TESFsHNUMOFG.LT. 1 .OR.HNUMDFG.GT.120 
IF  (TEST)  GO  TO  5555 
IHOID-MOD  (120  ,HNUM0FG) 

IF  (MOLD.!®.  0)  TYPE  "Try  again.  Input  must 
$  divide  evenly  into  120." 

IF (IHOLD.I'E.O)GO  TO  556 

558  ACCEPT”Enter  vertical  #  of  grid  rectangles 
$  (1-30):  ",VNUOPG 
MOLD-MOD  (30  ,VNUM0FG) 

TEST-VNUM0FG.LT.1.OR.VNUMDFG.GT.30 
IF  (TEST)  TYPE  "Out  of  range!<7><15>" 

IF  (TEST)  GO  TO  558 

IF(IHQLD.lC.O)TYFE”Try  again.  Input  must  divide 
$  evenly  into  30." 

IF(IHQLD.tG.O)GO  TO  558 
GO  TO  3 


I) 
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500  TYPE  "Lower  right  quadrant  option  on." 

IMDt*32 
IB0W“2 
JMHt-0 
TW0C0UMIV512 
KSTART=128 
00  TO  2 

1  TYPE "Upper  left  quadrant  option  on.* 

HNCM0PG=4 
VNUM0PG=6 
3  IRCW=1 

IMDH) 

JMI»*512 

TWXXXJNT=0 

KSTART=0 


******  qreATE  the  normalized  pile  ******************************** 

2  IMAX=IMIN+31  ;limits  on  infile  KDBLK  (lower) 
JMAX=*JMIN+-511  /limits  on  outfile  zeroed  blocks 
CALL  TBER(ISTART)  /start  timer 

TYPE"<15>" , "Creating  the  normalized  file." 

ONEOOUNF=0  /workfile  WRBLK  counter 

SUMSQ-0.0 

CZERD=CMHjX(0.0,0.0) 

DO  10  J=l,1024 

10  CN0FM  ( J)  =CZERO 

DO  20  J=JMIN,JMAX,16  /zero  appropriate  outfile  rows 

CALL  WRBLK(2,J,CN0RM,16,IER) 

20  IF(IER.bE.l)TYPE"CNC*M  zero  WRBLK  error:", IER 

IF(IMIN.EQ.32)QO  TO  200  /switch  was  L 

***  NORMALIZE  UPPER  QUADRANT  *** 

Create  the  unpacked  workfile  and  determine  the 
energy  content  of  it.  Work  on  4  picture  rows  per  loop. 

HWIDTH-120/HNUMDPG  /width  of  rectangle 

BFERSRU>3 0/VNUMOrc  /RDBLKS  per  grid 

TLE2)GnBs23.0  /length  of  reduced  template 

TWIDTH-47.0  /width  "  "  " 

VLEK7IH“4.0 "FLOAT (BPERGRID)  / 4  rows/block 

TAREA*ILENGTH*IWIDra  /template  area 

GAREAHFLQAT ({WIDTH)  *VLEHGTH  /rectangle  area 

AREAFACTOR»TAREA/GAREA 

DO  19  K»l,30  /initialize  energy  terns 

DO  19  J-1,120 

19  ENERGY(J,K)-0.0. 

Determine  energy  in  30  blocks  (120  rows) 

HGL0CK-IMIN-1  /set  RDBLK  counter 

DO  40  VGRII>1  ,VNUMDFG  /do  5  rows  of  grids 
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DO  39  GRDBLK-l ,BPERGRID  ;of  6  blocks  each 

INBLOCK= BBLOCK+1  ;REBLK  counter 

CALL  RDBLK(0,HBLOCK,ARRAY,1,IER) 

KK-IM3N 

DO  25  11*1,4  ;do  4  video  rows 

DO  21  J-1,2  ;set  first  8  oolurms 

KK-KK+1  jto  zero 

;  (noise  terms) 

CONTINUE 

DO  23  J=3 ,32  ; arrange  non-zero  portion 

KK-KK+1  ;of  picture  for  processing 

ROW(J,II,IROW)=ARRAY(KK) 

CONTINUE 

KK=KK+32 

CONTINUE 

CALL  UNPACK (128, RCW (1, 1, IRCW) , UNPACKED) 

Determine  energy  in  4  rows 


26 

28 

30 


39 

40 
C 

C 

C  3000 

C 

C 

C  2000 
C  100 
C 
C 
C 

666 


45 

50 


MIN00Ir=9 

DO  30  HGRH>1  ,HNCMDPG  ;do  5  oolurms  of  grids, 

MAXO5L=MINa0L+WIDra-l 

DO  28  KCQL*MINCOL,HAXCOL 

DO  26  JRGt^0,3  ;do  4  rows 

ENERGY  (HGRID,VGRID)  -ENERGY  (HGRID,  VGRID)  + 

(FLOAT (UNPACKED (KCOL+JRCW*128)  +1) )  **2 
CONTINUE 
CONTINUE 
MINOOL=MAXCOL+1 
CONTINUE 

CALL  WRBLK(1,ONEOOUNT,UM3ACKED,2,IER) 

IF(IER.I«.1)TYPE"1RDBLK  *",CNEO0UNT,"  error:", IER 
ONEGOUNIXUEOOUNIM-2  jWRBLK  counter 
CONTINUE  ;do  next  block  of  grids 

CONTINUE  ;do  next  grid  row 

IF  (HNUM0FG.GT.5)  00  TO  666 

WRITE  (12,3000)  INFILE  (1) 

FORMAT ("  ENERGIES  OF  ",S13,///) 

DO  100  VQRII>=1  ,VNLMDFG 

WRUE (12,2000)  VGRID,  (ENERGY (HGRID,VGRID)  ,HGRID=1  ,HNUM0FG 
FORMAT ("  GRID  RCW",I2,5(10X,F12.2) ,/) 

CONTINUE  ' ' 

Determine  the  normalization  factors 

DO  50  VGRID»1,VNUM0FG  • vertical 

DO  45  HGRII>1  ,HIW1)FG  ; horizontal 

IF  (ENERGY  (HGRID,  VGRID)  .LE.l  .0)  ENERGY  (HGRID,  VGRID)  *=1.0 
NORM  (HGRID,  VGRID)  -SORT  (ENERGY  (HGRID,  VGRID)  "AREAFACTOR) 
CONTINUE 

CONTINUE 


C 

C 


Normalize  and  create  the  output  file 


INBLOCK*- 2  ;RDBLK  counter 

DO  70  voin>i,vT«Drc 
DO  60  GREBLK=1,BPEHGRID 
HBLOCK-INBLOCK+2 

CALL  REBLK(1,DBL0CK, UNBACKED, 2 fIER) 

MNCDL-9  ; starting  column 

ICNT**0 

DO  58  BSRH>*1  f HNOMOFG  ;do  portion  of  5  grid  blocks 

MAHEL-MNOOL+IWinrH-l  ;width=24 
DO  56  KCOL=MINOOL , MAXCOL  ;do  width  of  1 

;grid  block 

ICNP=ICNPtl 
DO  54  JRCW=0 ,3 

OUTPUT=FLQAT ( UNPACKED ( ICNT+JR0N*128) +1) 

$  /NORM  (HGRID,VGRID) 

CNORM  (KOOL+JPCW*256)  =CMPLX  (OUTPUT,  0 .0) 

54  CONTINUE 

56  CONTINUE 

MINaOWBVXCOLfl 
58  OONTINUE 

CALL  WRBLK(2,THOOOUNT,GGEM,16,IER) 
IF(IER.tE.l)TYPE"2WRBLK  #",TOOCOUNr, "  error:", IER 
TWOCOUN^*IWDCOUNP+-16 
60  CONTINUE 

70  OONTINUE 

DO  80  J»l,1024 

CNORM  ( J)  =CZERO 

80  CONTINUE 

DO  90  J=l,2  ;zero  out  noise  rows 

CALL  HRBLK(2,'MXXXJNT,CNOFH,16,IER) 
IF(IER.NE.1)TYPE"2HRBLK  #",TWOOOUNT,"  error:", IER 
1WOOOUNMWOGOUNrH6 


90  OONTINUE 

GO  TO  555 
C 

C  ****  Normalize  lower  quadrant  (template)  **** 

C 

200  DO  240  I-IMIN,IMAX  ;read  non-zero  portion  of  infile 
CALL  RDBLK (0,1, ARRAY,  1, IER) 
IF(IER.1C.1)TYFE"0RDBLK  #",I,"  error:", IER 
C 

C  Set  nonzero  portion  of  ARRAY  equal  to  RON 

C 

KK-IMIN 
DO  225  11-1,4 

DO  223  J-1,32 
KK-KK+1 

223  R0W(J, II, IRGN)  “ARRAY (KK) 

225  KK-KK+32 

C 

CALL  UNPACK  (128  ,RCN(1,1,IRCN) , UNPACKED) 

DO  230  J- 1,51 2  /determine  energy  in  4  rows 
230  SUMSOSUM9QH1NEACKED  { J)  **2 

CALL  WRBLK(l,ONBOOUNr,UIPACKED,2,IER) 
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240 

C 


IF (IER,bC.l)T5fPE"lWRBLK  #",CNBOOUNT,"  error:", IER 

CNExruNP-a®ixjNrf2 

******************************************************* 

TENEBGY=SQRT  (SUMSQ)  ;the  normalizing  factor 

******************************************************* 


*  Normalize  the  significant  pixels  of  INFILE  **************** 


250 

260 


270 


DO  270  1=0,62,2 

KK=KSTART  ; starting  column  of  nonzero  outfile 

ICNT=0 

CALL  EDBLK(1, 1, UNPACKED, 2, IER) 

IF(IER.EE.1)TYPE"1RDBLK  #",I,"  error:", IER 
DO  260  K=l,4  ; normalize  infile 

DO  250  J=l,128 
KK=KK+1 
ICNT=ICNT+1 

OtrTPUT=UNEACKED  ( ICNT)  /TENERGY 
CNOFM(KK)  KMPLX  (OUTPUT,  0 .0) 

KK=KK+128 

CALL  WRBLK(2,TOOOOUOT,CNOIM,16,IER) 

IF(IER.tC.l)TYPE"2WRBLK  #",™OG0UNT,"  error:",IER 
1WOOOUMMWOQOUNIM-16 


c****************************************************************** 

C 

555  CALL  TIMER  (ISTOP)  ;stop  timer 

WRITE  (10 ,1000)  OUTFILE  (1) 

1000  FORMAT ("  The  normalized  file  is  in  — >  ",S13) 

CALL  RESET 

CALL  DFILW( "TEMP", IER) 

IF  ( IER.  1C.  1)  TYPE  "TEMP  DFILW  error:",IER 
STOP 
END 
C 

C************  program  NORMALIZE  ************************************* 


0******************************************************************** 

c 

C  Program  CMUITIILY  Written  by  Lt.  Jim  Crane r 

C  Fortran  5 

C 

C  This  program  performs  a  point-by-point  complex  multi- 

C  plication  between  INFILE1  and  OONJG(INFILE2) .  The  input 

C  files  must  be  256x256  point  complex  arrays;  the  output 

C  file  is  a  256x256  point  complex  array. 

C 

C  Execution  Line  Format: 

C  CMULTIHjY  infilel  infile2  outfile 

C 

C  Load  Line  Format: 

C  RLDR  CMJUTIHjY  IOF  @FLIB@ 

C 

C******************************************************************** 

INTEGER  FILE (7)  ,INFILE1{7)  ,DJFILE2(7)  ,MAIN(7) 

COMPLEX  MAI (2048) ,MA2(2048) fMA3(2048) 

COMMON  FILE,  INFILEl ,  INFILE2  ,MAIN 
C 

q ********  j/q  ftt.e  MANAGEME24T  **************************************** 

C 

CALL  I0F(3,MAIN,INFILE1,INFILE2,FILE,MS,I2,I3,I4) 

WRITE (10,1999) INFILEl(l) ,INFILE2(1) ,FILE(1) 

1999  FORMAT ("  INl=",S13,"  IN2*",S13,"  0UT*",S13) 

CALL  OPEN(l,FILE,2,IER) 

IF  (IER.EQ.l)  QO  TO  55 
CALL  CFILW(FILE,3,1024,IER) 

IF(IER.EQ.41)CALL  CFILW (FILE, 2 , IER) 

IF  (IER.  IE.  1)  TYPE  "OUTFILE  CFILW  error  *",IER 
CALL  OPEN  (1  ,FH£,2  ,IER) 

IF  (IER.IE.l)  TYPE  "OUTFILE  OPEN  error  #",IER 
55  CALL  0PEN(2,INFILE1,2,IER) 

IF(IER.IE.1)TYFE"INFILE1  OPEN  error  #",IER 
CALL  OPEN(3,INFILE2,2,IER) 

IF(IER.NE.1)TYPE"INFILE2  OPEN  error  #",IER 
C 

C********  Perform  point— by-point  multiplication  *********************** 
C 

DO  30  1*0,992,32  ; process  1024  blocks 

CALL  KEBLK(2,I,MA2,32,IER)  ;read  8  complex  rows 

IF(IER.IE.1)TYPE"2RDBLK  #",I,"  error:", IER 
CALL  KDBLK(3,I,MA3,32,IER) 

IF(IER.NE.1)TYH:"3RDBLK  *",I,"  error:", IER 
DO  20  K-1,2048 

MAI  (K)  =MA2  (K)  *C0NJG(MA3  (K) ) 


20  OONTINUE 

CALL  WRBLK(1,I,MA1,32,IER) 
IF(IER.IE.1)TYFE"WRBLK  #",I,"  error:", IER 

30  OONTINIE 

WRITE (10,40) FILE (1) 

40  FORMAT ("  ",S13,"  created  by  CMUITIILY") 

END 

C************  program  OfJITIILY  *************************' 


c******************************************************************** 

c 

C  Program  CHOI  Written  by  Lt.  Jim  Cromer 

C  Fortran  5 

C. 

C  This  program  converts  a  256x256  complex  file  into  a 

C  256x256  integer  file,  The  real  part  only  of  the  complex 

C  file  is  saved;  the  imaginary  part  is  assumed  to  be  zero. 

C  The  maxinun  value  and  its  position  are  written  to  first 

C  3  words  of  block  255.  [Values  greater  then  2  are  set 

C  to  1.  All  other  values  are  divided  by  two.] . 

C 

C  Execution  Line  Format: 

C  CTOI  complex  infile  integer  outfile 

C 

C  Load  Line  Format: 

C  RLDR  CTOI  IOF  TIMER  @FLIB@ 

C 

c******************************************************************** 
INTEGER  MAIN (7) ,  INFILE  (7)  , OUTFILE  (7)  ,MST' 

INTEGER  OUTINTBGER(4096)  ,OOLUMNNUMBERrRL  iNUMBER 
COMPLEX  INQ0MPLEX(4096) 

C 

C********  i/q  ftt.r  management  ************************************* 

c 

CALL  IOF(2, MAIN,  INFILE, OUTFILE,  IDUM, MS,  ISl,IS2,IS3) 

CALL  TDER(O)  ;start  timer 

CALL  OPEN  (0,  INFILE,  1,  IER) 

IF  ( IER.  1C.  1)  TYPE  "INFILE  OPEN  error  #",IER 
CALL  OPEN (1 , OUTFILE, 3  ,IER) 

IF(IER.EQ.l)GO  TO  1 

CALL  CFILW (OUTFILE, 3, 256, IER) 

IF (IER.EQ.l)TYEE "Created  a  oontigous  output  file." 

IF ( IER. EQ. 41) CALL  CFILW  (OUTFILE, 2, IER) 

IF  ( IER.  NE.l)  TYPE  "OUTFILE  CFILW  error  #",IER 
CALL  OPEN (1, OUTFILE, 3, IER) 

IF (IER. NE.l) TYPE "OUTFILE  error  *",IER 
C 

C*********  CONVERT  COMPLEX  WORDS  TO  INTEGER  ********************** 

C 

1  IMAX=0 

DO  20  J=0, 960,64 

CALL  BDBLK(0,J,INOOMILEXr64,IER) 

IF(IER.IE.1)TYPE"INOOMHjEX  RDBLK  error:", IER 

DO  10  K=l,4096 

AREAL-REAL ( INOOMPLEX (K) )/2.0 
OUTINTEGER  (K)  -INT  (AREAL*32767 . 0) 

IF  (AREAL.GE.l  .O)OUTIMrEGER(K)  =32767 
IF(OUTINTEGER(K)  .LT.IMAX)Q0  TO  10 
OOLUMNNUMBE»*K-(IOT((K-l)/256)  56) 
RCWNUBER-INT(  (K-D/256)  +14J/4 
IMAX^XJTINTEGER(K) 

10  CONTINUE 

17  CALL  WRBLK(1,  (J/4)  ,OUTINrE3ER,16,IER) 
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IF(IER.^.1)TYPE"00TINTBGER  WRBLK 
$(J/2) ,"  error  #",IER 
20  ODNTINIE 

OUTTNTEGER(l)  =IMAX 
£X7TINTBGER(2)  OOLUMNNUMBER 
OUTTNTBGER(3)  =RCWNUMBER 
CALL  WRBLK  (1,255 , OOTINTEXSER,  1 ,  IER) 

IF(IER.lE.l)TYre"IMAX  WRBLK  error:", IER 
CALL  TIEER(l)  ;stop  timer 

TYH)"  The  maximum  integer  value  is:",IMAX 
TYPE"  @  RCWNIMBER=",RCWNUMBER 

TYPE"  @  COLUMN  =",  OOLUMNNUMBER 

RMAX»(FLOAT(IMAX)  )/32767.0 
TYPE"  IMAX/32767  =",RMAX 

WRITE  (10 ,1000)  INFILE  (1)  ,0UTFILE(1) 

1000  P0BMAT("  The  complex  file  ",S8,"  has  been 
$oonverted  to  the  integer  file  ",S8) 

CALL  RESET 
STOP 
END 
C 

C*********  Program  C7POI  ********************************************* 


c**************************************************************** 

c 

C  Program  IMUITIELY  Written  by  Lt.  Jim  Cromer 

C 

C  This  program  calculates  the  point-by-point  geometric 

C  mean  between  INFILEl  and  INFILE2;  the  product  is  output 

C  to  OOTFILE.  All  files  must  be  256x256  integer  files. 

C 

C  Execution  Line  Format: 

C  IMULTIELY  infilel  infile2  outfile 

C 

C  Load  Line  Format: 

C  RLDR  IMJKPIILY  IOF  @FLIB@ 

C 

C************************************************************** 

INTEGER  MAIN  (7)  ,INFILE1(7)  ,INFILE2(7) , OUTFILE (7) 

INTEGER  GS(2)  ,LS1(2)  ,LS2(2)  ,LS3(2) 

INTEGER  FACP0R1(8192)  ,FACTOR2  (8192) ,  PRODUCT  (8192) 

C 

C********  j/o  FILE  MANAGEMENT  ************************ 

C 

CALL  IOF (3 ,MAIN, INFILE1 , INFILE2 , OUTFILE, GS, LSI fLS2 ,LS3) 

JEEfc=l 

KER=1 

CALL  OPEN (1 , INFILE1 ,2, IER) 

IF  (IER.NE.l)  STOP"  INFILE1  OPEN  ERROR" 

CALL  OPEN (2,INFILE2 ,2,IER) 

IF  (IER.NE.l)  STOP" INFILE2  OPEN  ERROR" 

C 

C  First  check  to  see  if  the  file  exists;  if  it 

C  doesn't,  try  to  create  a  contiguous  file. 

C 

CALL  OPEN (3, OUTFILE, 2, IER) 

LER=IER 

IF  (IER.  EQ.  13)  CALL  CF ILW  (OUTFILE,  3 ,256,  JER) 

IF (JER.NE.1)CALL  CFILW (OUTFILE,  2  ,KER) 
IF(KER.NE.1)ST0P"0UTFILE  CFILW  ERROR" 

IF  (KER.EQ.1  .OR.  JER.EQ.l)  JEBrl 
IF  (IER.IC.l)  STOP"OUTFILE  OPEN  ERROR" 

IF(LER.NE.1)CALL  OPEN (3, OUTFILE, 2, IER) 
IF(IER.NE.l)STOP"OUTFILE  OPEN  ERROR" 

C*******  Perform  point-by-point  multiplication 

C 

C 

DO  20  1*0,224,32 

CALL  RDBLK(1,I,FACP0R1,32.TER) 
IF(IER.NE.1)TYPE"1RDBLK  #",I,"  error:", IER 
CALL  PDBLK(2,I,FACT0R2,32,IER) 
IF(IER.NE.1)TYPE"2RDBLK  #",I,"  error:", IER 
DO  10  J-1,8192 

AH0UK9QFT  ( ( FLOAT (FACPOR1  ( J) )  "FLOAT  ( FACTOR2  ( J) ) ) ) 
PRODUCT  (J)*INr(AHQLJ>H). 5) 

10  CONTINUE 

CALL  WRBLK(3,I,JR0DUCT,32,IER) 


f  *  s 


IF(IER.NE.1)TYPE”WRBLK  #",I,"  error:", IER 
TYPE"  ***  DATA  CRUNCH  I  ***<11><11>***  EATA  CRUNCH  1  ***" 

20  CONTINUE 
C 

C******  Write  message  to  the  CRT  terminal  ********* 

C ' 

WRITE  (10, lOOO)OUTFILE(l) 

1000  FORMAT  (///,  "The  integer  product  file  is  named  — >  ",S13,/) 

CALL  RESET 
STOP 
END 
C 

C********  program  MULTIPLY  ******************************************** 
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APPENDIX  F 


PROCESS  EVALUATION 


This  appendix  contains  the  following  programs 

1.  ITOC 

2.  PEAK  (FI) 

3.  CTOV 

4 .  DISTANCE  (EUCLID) 
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C******************************************************************** 

Program  ITOC  Written  by  Lt.  Jim  Cromer 

Fortran  5 

This  program  converts  a  256x256  integer  file  into 
a  256x256  complex  file. 

Execution  Line  Format: 

IT0C/[A,E,H,N,  or  0]  integerf ile  [/C  and  or  T]  complex  [/M] 
Global  switches  — >  A:  convert  all  values  >0 

E:  convert  values  >80%  of  the  maximum 
;values  not  con-  H:  convert  values  >1/2  maxinum 

;verted  will  be  N:  convert  values  >90%  maxinum 

;set  to  zero  0:  accept  other  conversion  value 

Input file  switch  — >  C:  "crunch"  256x256  data  into  256x128 

array  for  use  with  PLTTRNS  contour 
T:  set  output  values  equal  to  input 
value  minus  the  threshold 

Outputf  ile  switch - >  M:  insert  maxinum  value  into  3rd 

column  of  every  4th  row 

Load  Line  Format: 

KLDR  ITOC  IOF  TIMER  @FLIB@ 

C******************************************************************** 

INTEGER  MAIN{7)  , INFILE (7)  ,0UTFILE(7)  ,MS(2)  ,IS2(2)  rISl(2) 

INTEGER  ININTBGER ( 40 96 )  ,00LUMN,RCW 
REAL  INREAL 

COMPLEX  OUTO0MPLEX(4096) 

C 

q ********  j/o  ptt.r  MANAGEMENT  ************************************* 

C 

CALL  IOF(2,MAIN, INFILE, 07TFILE,IIXJM,MS,IS1,IS2,IS3) 

PERCQH^-9999  0 

IF(ITEST(MS(1),15).EQ.1)PERCEN]>0.0  ; switch  was  A 

IF(ITEST (MS  (1)  ,11)  .EQ.l)FERCHn>*0.80  ; switch  was  E 
IF(ITEST(MS(1)  ,8)  .EQ.l)  PERCENTS.  50  ;switch  was  H 
IF(ITEST(MS(1)  ,2)  .BQ.1)PERC0*W).9O  ;switch  was  N 

2  IF  (ITEST  (MS  (1)  ,1)  .EQ.l) ACCEPT"  Enter  the 
$  threshold  percent  (0.0  -  .99):  ", PERCENT 
IF  (PERCENT. LT. 0.0  .OR. PERCENT. GT.1.0)  GO  TO  2 
IF (PERCENT. EQ. -9999.0) STOP "BAD  GLOBAL  SWITCH" 

CALL  TIIER(O)  ;start  timer 

CALL  OPEN (0, INFILE,  1,IER) 

IF (IER.JE.l) TYPE* INFILE  OPEN  error  #",IER 
CALL  OPEN(l,OUTFILE,3,IER) 

IF(IER.BQ.1)G0  TO  1 

CALL  CFILW(OUTFILE,3,1024,IER) 

IF  (IER.EQ.1)  TYPE  "Created  a  contiguous  output  file.” 
IF(IER.EQ.41)CALL  CFILW((XJTFILE,2,IER) 

IF(IER.W!.1)TYPE"00TFILE  CPUM  error  #",IER 


o  o  o  r>  n  n  r>  non 


CALL  OPEN  (1  ,0OTFILE,  3 ,  IER) 
IF(IER.2C.1)TYPE"0UTFILE  error  #",IER 


C 

C*********  CONVERT  INTEGER  WORDS  TO  COMPLEX  ************************** 

c 

1  CALL  RDBLK(0,255,ININrBGER,l,IER) 

IP(IER.IE.1)TYPEBIMAX  RDBLK  ERROR:  ",  IER 

IMAJMNIOTBGERd) 

G0LIM^ININTBGER(2) 

RCW-ININTEGER(3) 

nBRESH»INT(FERCENr*FLQAT(IMAX) ) 

JFACTOR-2 

JBL0CK=€4 

INCREMENTS) 

IF(IS1(1)  .BQ.8192)  JEACT0R=1 
IF(IS1(1) .EQ.8192) JBLOCK*32 
IF(IS1(1)  .EQ.8192)  INCR0ffiNP=256 

Create  the  output  file 

DO  20  J=0,480,32 

CALL  REBLK(0r  (J/2)  ,ININTBGER,16,IER) 

IF  (IER.  NE.l)  TYPE"  IN  INTEGER  RDBLK  #"r(J/2)f"  error:", IE 

LL-0 

KK=0 

DO  10  1*1 ,  ( 8*JEACTOR) 

DO  5  K=l,256 

KK*KK+1 

LLr=LL+-l 

IF (IS1  (2)  .EE. 4096) GO  TO  3 

INREAL=FLQAT  ( ININTEGER  (KK)  -riHRESH)  /327 67 . 0 

GO  TO  4 

1  INREAL=FLQAT(INIMrEGER(KK)  )/32767.0 

1  IF  (ININTBGER(KK)  .LT.  IIHRESH)  TNREAL=0 .0 

OUTOOMHjEX  (LL)  *OUUC  ( INREAL,  0 .0 ) 

5  CONTINUE 

KK*KK+ INCREMENT 
L0  CONTINUE 

Prepare  output  for  PLTTRNS  row  plot 

IF(IS2(1)  .EQ.8.AND.IS1(2)  .EQ.4096)CUTC0MILEX(3) 

$KMPLX ( FLOAT ( IMAX-ITHRESH)  ,0.0) 

IF(IS2(1)  .EQ.8.AND.IS1(2)  .EC. 4096) 
$CXJTOOMH-EX(3)KXEliX(FLQAT(IMAX)  ,0.0) 

CALL  WRBLK(1,  (JMFACTOR)  ,OUTCOMPLEX, JBDXK,  IER) 
IF(IER.IE.1)TYFE"0UT00MFLEX  WRBLK  #",  (JMFACTOR) , 

$"  error  :",IER 
20  CONTINUE 

The  data  will  be  compressed  to  the  front  half  of  the 
output  plane. 


IF (ISl(l).l®. 8192)00  TO  35 


non  o  n o 


DO  23  1-1,3 

OUTOQMHiEX  (1+1792 )  <MHiX  (0 .0 ,0 .0) 

23  CONTINUE 

CALL  WRBLK(1,480,OUT00MHjEX,32,IER) 

IF(IER.NE.1)TYPE"WRBLK  (DATA)  *480  error:", IER 
DO  25  >1,4096 

OOTOOMHiEX  (I)  «OULX  (0 .0 ,0 .0) 

25  CONTINUE 

DO  30  J-512,960,64 

CALL  HBBLK(1,J,0UT00MR£X,64,IER) 
IF(IER.N3.1)TYPE"WRBLK  #",J,"  error:", IER 

30  CONTINUE 

GO  TO  40 

Emended  output  only 

35  DO  38  >1,3 

OUT00MELEX  (1+3840)  -CMPLX  (0 .0 ,0 .0) 

38  CONTINUE 

CALL  HBBLK(l,960,CUT00HHiEX,64,IER) 

IF(IER.NE.1)TYFE"WRBLK  (DATA)  *960  error:", IER 

Send  completion  message  to  CRT 

40  CALL  TIMER(l)  ;stop  timer 

TYPE"IMAX=",IMAX 
TYPE"ITHRESB=",  ITBRESH 
TYPE "PERCENT=", PERCENT 
TEMP=FLQAT(ITHRESH)/32767.0 
TYPE"  Normalized  threshold-" , TEMP 
WRITE (10,1000)  INFILE (1)  ,OUTFILE(l) 

1000  FORMAT ("  The  integer  file  ",S8,"  has  been 

$oonverted  to  the  complex  file  ",S8) 

CALL  RESET 
STOP 
END 
C 

C*********  program  ITOC  ********************************************* 


c******************************************************************* 

c 

C  Program  PEAK  Written  by  Lt.  Jim  Crane r 

C  Fortran  5 

C 

C'  Biis  program  searches  a  256x256  integer  file  for  isolated 

C  regions  of  which  all  values  are  above  a  given  threshold. 

C  If  the  input  array  is  thought  of  as  a  3-dimensional  sur- 

C  face,  then  these  regions  will  be  the  "peaks"  of  the  surface. 

C  The  position  and  value  of  both  local  and  globed  peaks  is 

C  determined. 

C 

C  Execution  Line  Format: 

C  PEAK 


Load  Line  Format: 

RLOR  PEAK  FI  @FLIB@ 

Flag  Values: 

PEAK  - > 

; holds  condition 
; of  global  peaks 


STATUS  - 

; holds  condition 
;of  row  peaks 


-1  -  peak  closed 
0  -  peak  unused 
1  -  peak  open 

-1  -  row  peak  matched 
0  -  unused  in  current  row 
1  -  row  peak  unmatched 


INSIDE  - >  .« 

determines  if 
;pointer  in  .: 

; interior  or  exterior 
;of  a  row  peak 


.TRUE.  -  previous  value 

tested  >  threshold 
.FALSE.  -  else 


C ******************************************************************* 
REAL  NORMALIZE  (10) 

INTEGER  INPILE (7)  , WIDTH  (10)  ,PC£NrMAX(10) , LENGTH  (10) 

INTEGER  VALUE  (256)  ,F1 , THRESH,  ISTART(IO)  ,ISTOP(10) 

INTEGER  PCW,OOLUMN,FVALUE,TEMPreAK,reCW(10) ,IMAX(10) 

INTBGER  MAXOOLtMN(lO)  ,PSTCP  (10,256)  ,PSTART(10,256) 

INTEGER  PWLUE(IO)  ,PCDLtMN(10)  ,STA!IUS(10) , PEAK (10) 

INTBGER  JMAX(10)  ,JMIN(10)  ,RANK(10)  ,IRANK(0:10) 

LOGICAL  ATEST,  INSIDE 
C 

C********  INITIALIZE  VARIABLES  ********************************** 

C 

IOP>0 

100  DO  200  1-1,10 

IRANK(I)-0 
RANK (I) -0 
PEAK (I) -O 
PVALUE(I)-0 
PGDLUfti(I)M) 

PROW  (I)M) 

DO  150  J-1,256 


PSTAKT(I,J)*257 

PST0P(I,J)»0 

150  CONTINUE 

200  CONTINUE 
C 

q********  ggj  I/O  PARAMETERS  ************************************** 

C 

IF(IOFT.EQ.2)GO  TO  250 

ACCEPT  "What  is  the  name  of  the  input  integer  file?  ” 

READ (11 ,1000) INFILE (1) 

1000  FORMAT  (S13) 

CALL  OPEN (0, INFILE, 2, IER) 

IF  (IER. IE.  1)  STOP" INFILE  OPEN  ERROR" 

CALL  RDBLK (0,255, VALUE,  1, IER) 

MAXIMUM=VALUE(1) 

TYPE "Absolute  max=",MAXIMJM 
GO  TO  250 

210  TYPE"Input  error.  Percentage  must  be  between  1-100. <7><15>" 

250  ACCEPT"Enter  integer  percentage  of  absolute  maximum 
$  to  be  included:  ",IPERCENT 

IF (IPERCENT.LT.l. OR. IPERCEMT.CT.  100)00  TO  210 
TORESB^INT  (FLOAT  (IIERCENT}  *FL0AT  (MAXIMUM) /1 00.0) 

TYPE"INTBGER  THRESHOLD*",  THRESH 
C 

C****************************************************************** 

c 

C  Loop  through  the  scan  and  matching  modules  254  times 

C  (test  all  but  first  and  last  rows) 

C 

DO  600  RGW*2,255  ytest  rows  2-255 

CALL  RDBLK (0,  (ROW-1)  , VALUE,  1, IER) 

IF  ( IER.  NE.l)  TYPE  "RDBLK  *",(RCW-1),"  error:", IER 
INSIDE=. FALSE. 

NUMIEAKS=0 
DO  300  1*1,10 

STATUS  (I)  *0 
300  CONTINUE 

C 

C********  SCAN  RCW  ************************************************* 

C 

C 

C  This  module  determines  if  the  scanning  pointer  is  in 

C  the  interior  or  the  exterior  of  a  row  peak.  When  a 

C  row  peak  is  encountered,  the  peak  counter  NUMPEAKS  is 

C  increased,  and  the  flag  STATUS  is  set  to  the  unmatched 

C  condition.  The  maximum  value  and  corresponding  column 

C  number  for  each  row  peak  is  stored. 

C 

DO  350  COLUMNS, 255 

FVALUE-F1  (VALUE  (COLUMN-1) , VALUE  (COLUMN) , VALUE  (ODLUMNU) ) 
IF  (FVALUE.LT. THRESH. AND.  .NOT.  INSIDE)  Q0  TO  350 
IF  (FVALUE.GE.THRESH.AND.  INSIDE)  GO  TO  320 
IF (FVALUE.®. THRESH.  AND.  .NOT.  INSIDE) Q0  TO  310 
IF  (FVALUE.LT. THRESH.  AND.  INSIDE)  OD  TO  330 


uyuouu  uu  ouuou 


310 


NUMPEAKS-NUMPEAKS+1  ; row-peak  counter 

STA3US  (NUMIEAKS)  *1  ; row  peak  opened, unmatched 

IMAX  (NIKEEAKS)  =0 
ISTGP  (NUMIEAKS)  *256 
ISTART (NUMFEAKS)  "COLUMN 
INSIDE-. TRUE. 

320  IF(IMAX(NUMPEAKS)  .GT.FVALUE)00  TO  350 

IMAX (NUMFEAKS) -FVALUE 
MAXCOLUMN  (NUMFEAKS)  OOLUMN 
GO  TO  350 

330  INSIDE=.  FALSE. 

ISTOP  (NUMPEAKS)  =C0LUMN-1 

350  CONTINUE 


If  no  values  above  threshold  were  found  in  the 
last  row  tested,  then  close  all  open  peaks. 

IF  (NUMPEAKS. NE. 0)00  TO  400 
DO  370  1-1,10 

IF  (PEAK  (I)  .IE. l)QO  TO  370 
PEAK  (I)— 1 
JMAX(I) -B0W-1 
370  CONTINUE 

00  TO  600 


100  CONTINUE 

********  attempt  to  match  ****************************************** 
Attempt  to  match  row  peaks  to  open  global  peaks. 

DO  500  1-1,10 

IF  (PEAK  (I)  .IE.  1)Q0  TO  500  ;  check  all  open  peaks 
PEAK(I)— 1  ;will  be  closed  unless  matched 

JMAX  (I)  -B0W-1 
DO  450  TEMFFEAK-1,NUMHSAKS 

IF  (STATUS  (TEMPFEAK)  .EQ.0)QO  TO  450 
ATEST=PSTQP  ( I ,  (RCW-1) ) . LT.  ISTART  (TEMPFEAK ) 

$  .OR.PSTART (I ,  (RCW-1) )  .GT.  ISTOP  (TEMPPEAK) 

IF(ATEST)GO  TO  450  ;did  not  match 
PEAK  (I) -1 

STATUS  (TEMPPEAK)  —1 

PSTART(I,ROW)-MIN(PSTART(I,ROW) ,  ISTART  (TEMPFEAK) 
PSTOP(I,B0W)  -MAX  (PST0P(I,  ROW) ,  ISTOP  (TEMPPEAK) ) 

IF  (PVALUE  (I)  .(21.  IMAX  (TEMPFEAK)  )QO  TO  450 
PVALUE  (I)  -IMAX  (TEMPKAK) 

PROW (I) -BOW 

PQDLUMJ(  I)  ^lAXOOLUMN  (TEMPFEAK) 

450  CONTINUE 


(3********  must  OPEN  A  NEW  GLOBAL  PEAK  ********************************** 
C 

C  Natch  unmatched  cow  peaks  to  unused  global  peaks. 


C 

DO  550  T1MPPEAK=1,NUMPEAKS 

IP  (STATUS  (TEMPIEAK)  .IE.1)  00  TO  550 
DO  510  1-1,10 

IF  (STATUS  (TEMPEEAK)  .MS. l)GO  TO  510 
IF  (PEAK  (I)  .EE. 0)00  TO  510 
STATUS  (TB1PEEAK)  —1 
PEAK  (I)  -1 
JMIN(I)=H0W 
JMAX(I)=H0W 

PSTQP(I,BOW)-ISTOP(TEMPPEAK) 

PSTART  (I  ,RQW)  =  1ST  ART  (TEMPPEAK) 

PVALUE  (I)  -IMAX  (T01PFEAK) 

PROW (I) “ROW 

POOLUMN  ( I )  =MAXOOLUMN  (TEMPFEAK ) 

TYPE "PEAK  #",I,"  START  RGW:",RGW 

510  CONTINUE 

550  OONTINUE 
600  OONTINUE 
C 

C ********  ELIMINATE  mjitiply  defined  peaks  ******************* 

c 

DO  650  1=1,9 

IF(PEAK(I) .EQ.O)GO  TO  650 
DO  620  J«I,10 

IF  (PEAK  (J)  .EQ.O)GO  TO  620 

IF(I.EQ.J)GO  TO  620 

IF (PVALUE (J) .fE.PVALUE(I) )GO  TO  620 

IF  (PROW (J)  .NE.PRCW(I)  )QO  TO  620 

IF  (POOLUMN ( J)  .NE. PCDLUKN (I)  )GO  TO  620 

TYPE"***  MULTIPLY  DEFINED  PEAK  F0U1D  ***" 

PEAK (J) =999 

JNIN(I)-NIN(JMIN(I)  ,JMIN(J) ) 
JMAX(I)«MAX(JMAX(I)  ,JMAX(J)) 

JSTART*JMIN(I) 

JSTQE*JMAX(I) 

DO  610  RGW-JSTART ,  JSTCP 
PSTOP(I,ROW)=MAX(PSTOP(I,RCW)  ,PSTQP(J,RCW) ) 
PSTART ( I , RCW)  =MIN ( PSTAi  T ( I , ROW)  ,PSTART(J,RCW) ) 
610  CONTINUE 

620  CONTINUE 

650  OONTINUE 

DO  700  1=1,10 

IF(PEAK(I) .EQ.999) PVALUE (I) =0 

700  CONTINUE 
C 

C********  SORT  PEAKS  ACCOIDING  TO  THEIR  MAXIMJM  VALLES  ********* 

C 

DO  706  K=1 ,10 

LOOUNT-O 
DO  704  1-1,10 


IF  (PEAK  (I)  .EQ.O)GO  TO  704 

IF  (PVALUE  (I)  .EQ*32767)  IVALUE (I) =32766 

ioouravo 

DO  702  J»l,10 

IF  (PEAK  (J)  .EQ.O)GO  TO  702 
IF(I.BQ.J)GO  TO  702 
ATEST=PVALUE(I) .EQ.PVALUE(J) 
IF(.NOT.ATEST)QO  TO  702 
DOOUNr=imjNis-i 
iaxjNi>ioouNr+i 

PVALUE ( J) =PVALUE ( J) -IOOUNT 

702  CONTINUE 

704  CONTINUE 

IF(LOOUNT.BQ.O)GO  TO  708 

706  CONTINUE 
708  IRANK(O) =32767 
DO  720  1=1,10 

IF(PEAK(I) «EQ.O)GO  TO  720 
DO  710  J=l,10 

IF  (PEAK  (J)  .EQ.O)QO  TO  710 

IF  (PVALUE  (J)  .GE.  IBANK  (1-1))  GO  TO  710 

IRANK (I)  =MAX (IBANK (I)  ,PVALUE(J) ) 

710  CONTINUE 

720  CONTINUE 

DO  800  1=1,10 

IF  (PEAK  (I)  .EQ.O)GO  TO  800 
DO  750  J=l,10 

IF  (IBANK  (J)  .BQ.OJGO  TO  750 
IF (PVALUE (I) .BQ.IRANK(J) )RANK(I)=J 

750  CONTINUE 

800  CONTINUE 
C 

C********  oqmhtpe  OUTPUT  VALUES  ******************************** 

C 

810  DO  830  1=1,10 

IF  (PEAK  (I)  .BQ.999)EEAK(I)=0 

IF(PEAK(I).EQ.O)GO  TO  830 

PCENIMAX  (I )  =INT  (PVALUE  ( I)  /  (0 .01*MAXIMJM)  +0  .5) 

NORMALIZE  (I)  =FLQAT(PVALUE  (I)  )/32767.0 

LENGUH  (I)  =JMAX  (I) -JMIN  (I) +1 

WIOTH(I)=0 

JSTART=JMIN(I) 

JSTOP=JMAX(I) 

DO  820  J=JSTART,JSTCP 

IHOLCNPSTOP  ( I ,  J)  -PSTART  ( I ,  J)  +1 
WIDTH  (I)  -MAX  (IHQLD, WIDTH  (I) ) 

820  CONTINUE 

830  CONTINUE 
850  ICH-10 

870  IF  (ICH.  BQ.  12) WRITE (12, 1500) 

1500  FORMAT  (////) 

WRITE  (ICH,  2000) 

IF  (ICH.BQ.10)  WRITE  (10,3000)  INFILE  (1) 

IF ( ICH. EQ. 12) WRITE (12, 4000) INFILE(l) 


134 


B3HRESH“FLCIAT  (THRESH) /32767.0 
WRITE  ( ICH, 5000)  R3HRESH,  IPERCENT 
WRITE  (ICH,  60 00) 

WRITE  (ICHf 7000) 

WRITE (ICH, 8000) 

2000  FORMAT (//,15X,"  ********************************************** 

§**************** »j 

3000  FORMAT (//,30X, "INTEGER  FILE  EVALUATED  - >  ",S13,//) 

4000  FORMAT (//,30X, "INTEGER  FILE  EVALUATED  - >  <10>",S13,//) 

5000  FORMAT(40X, "THRESHCLE^" ,F5 .3 ,/36X, "%  OF  MAX  PEAK:", 15,//) 

6000  FORMAT  (2IX,  "PEAK  %MAX",33X,  "NORMALIZED") 

7000  FORMAT  (2 IX, "  #  PEAK  ROW  OOLLMN  WIDTH  LENGTH 

8PVALUE") 

8000  FORMAT  (21X,” -  -  - - -  - 

$ - -) 

DO  910  1-1,10 

DO  900  J=l,10 

IF  (RANK  ( J)  .NE.I)QO  TO  900 

WRITE  (ICH,  9000)  I,  PCENTMAX(J)  ,PROW(J)  ,POOLtMN(J) , WIDTH  (J) 
$LENGTH(J)  ,  NORMALIZE  (J) 

9000  F0RMAT(15X,I9,I7,I7,I8,I7,I7,F12.3) 

C  9000  P0RMAT(19X,I5,2X,I5,2X,I5,2X,I6,1X,I6,2X,I5,1X,F11.3) 

900  OONTINUE 

910  OONTINUE 

WRITE (ICH, 2000) 

IF(ICH.EQ.12)GO  TO  950 

ACCEPT"Enter  a  1  bo  send  results  to  the  linepr inter:  ",I 
IF(I.NE.l)GO  TO  950 
ICH-12 
GO  TO  870 

950  T5fPE"<15>" ,  "What  next?" 

TYPE"<15>","  Here  are  the  options:  <15>" 

TXPE"<15>", "<11>1  —  Try  a  new  input  file" 

TYPE"<15>",  "<11>2  —  Try  another  threshold  value" 

TXPE"<15>" ,  "<11>3  —  STOP<15>" 

970  ACCEPT"Enter  option - >  ",IOPT 

I=IOET 

IF(I.LT.1.0R.I.OT.3)T«E"<15>", "Input  error <7><7>1<15>" 

IF  (I.  LT.l.  OR.  I.  GT.  3)  GO  TO  970 
IF(I.EQ.2)G0  TO  100 
CALL  RESET 
IF(I.EQ.l)GO  TO  100 

TXPE"<15>","***  EXITING  PROGRAM  PEAK  ***<15>" 

STOP 

END 

C 

q**********  program  PEAK  ********************************************* 


HUGGER  FUNCTION  Fl(MJEFOREfN,NAFTER) 
c************************************************************ 

C 

C  Function  FI 

C 

C  This  function  is  part  of  the  program  PEAK.  In 

C  PEAK,  a  function  of  N  F1(N)  is  compared  to 

C  a  threshold  of  sane  %  of  the  maximum  value  to 

C  determine  if  a  local  peak  was  found. 

C 

C************************************************************ 

F1=N 

C  Other  possible  functions: 

C 

C  Fl=mr  (FLOAT  (NBEFORE+NtNAFTER)  /3 . 0+0 .5) 

C  F1=IOT  (FLOAT  ( NBEFORE+2  *W-NAETER)  /4 .0+0  .5) 

C  Fl^IKT  (FLOAT  (IBEFORDf  3*NfNAFTER)  /5 . 0+0  .5) 

RETURN 

END 

C 

q*********  Function  Fl  *************************************** 


Q ******************************************************************** 

C 

C  Program  CIW  by  Lt  Jim  Cromer 

C  Fortran  5 

C, 

C  This  program  converts  a  complex  input  file  (imaginary  part 

C  assumed  zero)  into  a  video  output  file.  The  input  file  is 

C  linearly  scaled  to  a  0-15  output  range.  Mininum  and 

C  maximum  values  to  be  included  are  input  by  the  user. 

C 

C  Execution  Line  For  rat: 

C  CTOV 

C 

C  Load  Line  Format: 

C  RLDR  CTOV  XWRBLK  TIMER  @FLIB@ 

C 

C********************************************************************** 

REAL  LOWER 

INTEGER  IARRAY(1024)  fdNFILE(7)  ,VOUTFILE(7) 

COMPLEX  CARRAY (1024) 

C 

C********  j/o  file  MMJAGEMENT  ********************* 

C 

ACCEPT"Enter  name  of  complex  input  file:  * 

READ (11 f 2000) CINFILE (1) 

2000  FORMAT  (S13) 

ACOEPT"Enter  name  of  video  output  file:  " 

READ  (11,2000)  VOUTFILE  (1) 

CALL  TDER(0) 

CALL  DFILW(VCUTFILE,  IER) 

IF(IER.1E!.1.AND.IER.NE.13)STGP"DFIUW  ERROR" 

CALL  CFILW(VCUTFILEf2fIER) 

IF  ( IER.  .  1 )  STQP*CFILW  ERROR" 

CALL  OPEN ( 1 f CINFILE , 2 , IER) 

IF(IER.IE.l)arOP"l  OPEN  ERROR" 

CALL  OPEN(2,VOUTFH£,2,IER) 

IF(IER.IE.l)arOP"2  OPEN  ERROR" 

C 

C*******  Determine  maximum  and  minimum  values  ********* 

C 

RMAX=0.0 
RMIN-99999.99 
DO  2  K=0f63 

CALL  RDBLK(1,(K*16) , CARRAY, 16, IER) 

IF(IER.NE.1)TYFE"1  RDBLK  #\(K*16),"  error:"fIER 
DO  1  J=lr1024 

AfREAL  (CARRAY  (J) ) 

RMAX=AMAX1  (RMAX,  A) 

RMIN=AMIN1  (RMINf  A) 

1  CONTINUE 

IF (MOD ( (K+l) ,4) .BQ.0)  TYPE  "BLOG?", (K*16) ,"  searched." 

2  CONTINUE 
C 

C ********  Determine  the  linear  scale  to  be  used  *********** 


TYPE  "Maximal  "  ,RMAX,  *  Mixuimra=",PMIN 

AOCEPT"Enter  maxinum  to  be  included : " ,  UPPER 
ACCEPT  "Enter  mininum  to  be  included: ",  LOWER 
SCALE=15 .999/ (UPPER-LOWER) 

C* 

C*********  Create  the  output  file  **************** 

C 

DO  20  K=0,63 

CALL  RDBLK(1, (K*16) ,CARRAY,16,IER) 

IF(IER.NE.1)TYPE"1  RDBLK  #",(K*16),"  error:", IER 
DO  10  J=l,1024 

IARRAY ( J) =15 
A=REAL  (CARRAY  (J) ) 

IF  (A.  LT.  LOWER)  IARRAY  ( J)  =0 
IF  (A. GE.  LOWER.  AND. A. LE. UPPER) 

$  IARRAY  ( J)  =INT  ( (A-LOWER)  *SCALE) 

10  CONTINUE 

CALL  XWHBLK  (2, K,  IARRAY,  1,  IER) 

IF(IER.NE.1)TYPE"2  WRBLK  #",K,"  error:", IER 

20  OONTINQE 
C 

C***********  said  cockle tion  message  to  CRT  ********* 

C 

CALL  TIMER(l) 

WRITE  (10 ,1000 ) VDUTFILE  (1) 

1000  FORMAT ("  The  video  file  created  is  called  ",S13) 

CALL  RESET 
STOP 
END 
C 

C**********  program  CTOV  *********************************************** 
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c********************************************************************* 

c 

C  Program  DISTANCE  Written  by  Lt.  Jim  Cromer 

C  Fortran  5  16  Oct  1982 

C 

C  This  program  accepts  as  input  a  template  file  window  and 

C  up  to  10  local  correlation  peak  positions  found  by  the 

C  program  PEAK.  Three  distance  factors  will  be  calculated 

C  between  the  template  window  and  9  scene  windows  (the 

C  center  window  corresponds  to  the  input  correlation  peak, 

C  the  other  8  are  its  nearest  neighbors) .  The  score  is 

C  computed  as  the  cube  root  of  the  product  of  the  factors. 

C  If  any  factor  is  less  than  zero  (corresponding  to  the 

C  measure  for  a  constant  gray  level  input  soene)  it  is 

C  set  to  zero.  Results  are  output  to  the  linepr inter. 

C 

C  FACTORS  USED: 

C  L1FACT=100  (1-NLl/NLLMAX) 

C  L2FACT=100 (1-NL2/NL2MAX) 

C  CFACT=100  (NXY-NXYMIN)  /  (1-NXYMIN) 

C  where  NLl  is  the  normalized  Ll  distance 

C  NL2  "  "  "  L2  " 

C  NXY  "  "  "  Cross-correlation  value 

C  NL1MAX  is  the  NLl  for  a  constant  gray  level  input 

C  NL2MAX  "  "  NL2  "  "  "  "  "  " 

C  NXYMIN  ""NXY""  "  " 

C 

C  Execution  Line  Format: 

C  DISTANCE 

C 

C  Load  Line  Format: 

C  RLDR  DISTANCE  TEST  XRDBLK  EUCLID  @FLIB@ 

C 

C********************************************************************* 

c 

REAL  NL1MAX 

REAL  NL1DIST,L1DIST,NSUMSQ,NDIST(10)  ,NDISTL1(100) 

REAL  QORFEAK(IOO)  ,DISTL1(10)  ,NXYMIN,NL2MAX,11j2 
INTEGER  INFILE1  (7)  ,INFILE2 (7)  , WIDTH, TB,a0MMENT(200) , 
$SB1,SLS,TTOP,TB1,SB,CFACTOR(100)  ,L2FACTQR(100)  ,UFACTOR(100) , 
$TLS,TLEFT, SCENE  (1024)  ,TEHP(1024)  ,SLEFT(100)  ,PGDLIMN(10) , 
$GQLCENT(100)  ,RCWCENT(100)  ,DIST(10)  ,CTOP(100)  ,PRCW(10) , SCORE 
LOGICAL  REDUCED, LTEST,SUPERESS 
COMMON  AIST1/  SCENE, TEMP, SLS,TLS,L1DIST, NUDIST 
COMMON  /LIST2/  LENGTH, WIDTH, IRCWaOUNT,SUMSQ,NSUfiQ 
COMMON  /LIST3/  SL1ENERGY,TL1ENERGY,SENERGY,TENERGY 
COMMON  /LIST4/  S1NOPM,S2NOPM,T1NORM,T2NOFM,CORREL 
LTEST-.TRUE. 

NL1MAX-0.46  ; normalized  Ll  distanoe  between  the 

; template  and  a  constant  gray  level 


C*******  i/q  pile  MANAGEMENT  *********************************** 


tJ 
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c 

c 

99  ACCEPT"Enter  tenplate  file  name  — >  " 

READ(11,1000)INFILE1(1) 

1000  FORMAT(S13) 

ACCEPT"Enter  scene  file  name  — >  " 

READ(11 ,1000) INPILE2 (1) 

CALL  0PEN(1,INFILE1,1,IER) 

IF ( IER. IC .  1 ) TYPE " INFILE1  OPEN  ERROR  #",IER 
CALL  OPEN(2,INFILE2,l,IER) 

IF(IER.tE.l)TYPE"INFILE2  OPEN  ERROR  #",IER 
C 
C 

C*******  enter  window  parameters  ******************************* 

c 

C  The  choice  is  given  to  compare  the  original  256x256 
C  pictures,  or  the  reduced  128x128  versions.  Reduced  saene 

C  files  are  assumed  to  occupy  the  upper  left  quadrant,  temp¬ 
er  lates  are  assumed  to  occupy  the  lower  right  quadrant. 

C 

REDUCEEK. FALSE. 

1  ACCEPT"<15>Enter  a  1  to  compare  original  video, 

$<15><11>  or  a  2  to  canpare  reduced  video:  ",I 
IF(I.LT.1.0R.I.GT.2)TYPE"<7>INFUT  ERRORl  <15>" 
IF(I.LT.1.GR.I.GT.2)G0  TJO  1 
IF  ( I.  BQ,  2)  REDUCED  .TRUE. 

IF  (LTEST)  GO  TO  5 

3  ACCEPT"Enter  a  1  to  change  tenplate  window  parameters:  ",IOPT 
IF(IOPT.ie.l)GO  TO  10 

5  ACCEPT" <15>","  Enter  top  row  of  original  tenplate 

$  window  (1-256)  :",TTOP 

ACCEPT"  Enter  left  column  of  original  tenplate 
$  window  (1-256) :",TLEFT 

ACCEPT"  Enter  width  of  window  (1-256) :", WIDTH 
ACCEPT"  Enter  length  of  window  (1-256) :", LENGTH 
LTEST=. FALSE. 

KDOP^ETOP 
MLEFT=lHiEFr 
MfmnwrajTH 
MLENGUMflJGTH 
GO  TO  10 

9  TYPE  "SORRY  <7  >.  Number  of  peaks  can  be  1-10  oily." 

10  ACCEPT" <15>", "Enter  #  of  candidate  peaks:  ",NUKPEAKS 

IF (NDMPEAKS. GT.lO.OR.IXjMPEAKS.LT.  1)Q0  TO  9 
TTOMCOP 
TLEFT=MLEFT 

wnmj-rwiDTH 
LENGTH«MLENGTH 
DO  20  II=1,NLMPEAKS 

TYPE"<15>","<15>", **********  PEAK", II,"  ************ 

GO  TO  15 

13  TYPE"<15>", "Sorry. <7>  Peak  row  must  be  1-256." 

15  ACCEPT" <15>", "Enter  peak  row  number:  ",HtGW(II) 

IF (FRCW(II)  .LT.l  .OR.FRCW(II)  .GT.256)G0  TO  13 
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CTOP (II) -256+TTGP-2*PRCW  (II) 

IF(RHXXED)CP0P(II)=128+IOT(  (TT0Pfl)/2)-FRCW(II) 

00  TO  19 

17  TYFE"<7>Peak  ooltran  must  be  1-256.  Try  again." 

19  ACCEPT" <15>" , "Enter  peak  oolunn  nunber:  ",P0DL11**(II) 
IF(P00LlMN(II)  .LT.1.0R.PGQLUMN(II)  .GT.  256)  00  TO  17 
SLEET  (II)  =256-WLEPT-2*POJLU1N(II) 
IF(REDUCED)SLEET(II)=128+INr(  (TLEET+l)/2)-FC0LUMN(II) 
RGWCENT  (II)  KTOP  (II)  +IOT  (LQ«7m/2) 

OOLCH7T  ( II)  =SLEET  (II) +INT  (WIEOE/2) 

IF  (REDUCED) RGWONT  (II)  =CK3P  ( II)  +IOT  (LENGTH/4) 

IF  (RHXXED)  OOLCENT  (II)  =SLEET  (II) +IMT  (WmiH/4) 

20  CONTINUE 
1*=-1 

IF  (REDUCED)  LENG1H=INT(  (LENGTH+D/2) 

IF  (REDUCED)  WH?TH=INr  ( (WIOTH+1 )  /2) 

IF  (REDUCED)  TTOP=INT  ( (TTOP+1 )  /2)  +128 
IF  (REDUCED)  TLEFT=INT  ( (TLEFT+1 )  /2)  +128 

ileft=tleet 

ITOMTOP 
iwimH=wiDra 
IU2*HH=IJENGTH 
SUPPRESS*. FALSE. 

ACCEPT  "Enter  a  1  to  suppress  window  messages:  ",  12 
IF  (12  .E9Q.1)  SUPPRESS*  .TRUE. 

TYPE  "Executing  .  .  .  ." 

C 

Q******************************************************************** 


This  section  computes  the  REBLK  and  EUCLID  search 
window  parameters,  then  computes  the  distance  measures 
for  each  of  the  windows  entered. 

11=10 

DO  600  J  J=1 , NUMPEAKS 
DO  600  J=l,3  ;DO  9  windows 

DO  600  K=l,3  ; compute  the  window  shift 

II-II+l 

COLCENT  (II )  KBLCENT  ( JJ)  -2+K 

CTOP(II)=CTQP(JJ) -2+J 

ROWCEWT  ( 1 1 )  =R0WCENT  ( J  J)  -2-hJ 

SLEET  ( 1 1 )  =SLEET  ( J  J)  -2+K 

TL1ENERGY=0.0  ; initialize  energies 

SL1EMERGY-0.0 

SENERGY=0.0 

TENERGY=0.0 

The  calls  to  TEST  check  to  see  if  the  input  parameters  are 
legal,  and  modifies  them  if  necessary: 

0  <  TOP  <  257,  (TOP  +  U2CTH)  <  258 

0  <  LEFT  <  257,  (LEET  +  WIDTH)  <  258 

CALL  TEST(TTOP,TLEFT) 

CALL  TEST(CraP(II)  ,SLEFT(II) ) 


nnnnn  non  , ,  non  non 


Set  RCBLK  and  EUCLID  paraneters 


40  ;if  1M.,  compute  the  energies 

;if  N=-l,  compute  the  distances 

IRO WOOUN3W) 

TB=INT (FLOAT (TTOP-1 ) /4 . 0 )  ;first  template  block  to  be  read 

SB»IOT(PDaAT(CrOP(II)-l)/4.0)  ;first  scene  block 
Nl=4CD(  (TTOP-1)  ,4) 

N2=M0D(  {CTOP(II)-l)  ,4) 

H£=flLEFT 

SLS»SLEET(II) 

User  check  of  window  parameters 

IF  (SUPHlESS)  GO  TO  45 
IF  (N.EQ.l)QO  TO  45 

TYPE"<15>" , ********************************************* 
§**********************« 

TYPE"<15>" , "<11><11>*******  WINDOW", (11-10) *******" 

WRITE (10 ,2000) INFILE1 (1) 

2000  P0FMAT(// ,10X,  "Template  file  name  - >  ",S13) 

WRITE (10 ,3000) INFILE2 (1) 

)000  FORMAT (10X,"  Scene  file  name  - >  ",S13,/) 

TYFE"<15>",  "<ll><ll>WIDrH=",WIDTO,  "<11>  LENGIH=",  LENGTH 

TYPE"  TEMHATE  TOP  RCW=",TTQP, "<11>  SCENE  TOP  RCW=",CTOP(I) 
TYPE"  TEMPLATE  LEFT  COLUMN*  ",TLS,"<11>  SCENE  LEFT  COLUMN 
$=",SLS,"<15>" 

Begin  evaluating  the  energy  (or  distances) 

15  ODRREL=0.0  initialize  distances 
NL1DIST=0.0 

CALL  XRDBLK  (1,18, TEMP,  1,IER) 

IF(IER.1E.1)TYPE"1HDBLK  #",1B,"  error:", IER 
CALL  XRDBLK (2, SB, SCENE,  1, IER) 

IF(IER.NB.1)TYPE"2RDBLK  #",SB,"  error:", IER 

This  module  will  continue  to  loop  until  the 
search  windows  have  been  completed  (i.e.  #  of 
iterations® (length  of  the  window)/4) 

100  CALL  EUCLID (NL,N2,N, $500) 

110  TB-TB+l 

SB-SB+1 

IF(N1.EQ.0)CALL  XRDBLK (1,TB, TEMP, 1 , IER) 
IF(IER.EE.1)TYPE"1RDBLK  #",TB,"  error:", IER 
IF(N2.EQ.0)GALL  XRDBLK (2, SB, SCENE,  1 , IER) 
IF(IER.NE.1)TYH!"2RDBLK  *",SB,"  error:", IER 
CALL  EUCLID (N1,N2,N, $500) 

IF(Nl.EQ.N2)GO  TO  110 
IF(N1.EQ.0)GALL  XRDBLK  (1,TB, TEMP, 1 , IER) 
IF(IER.IB.1)TYPE"1RDBLK  #",TB,"  error:", IER 
IF(N2.EQ.0)Q\LL  XRDBLK ( 2, S8, SCENE,  1, IER) 
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IF(IER.EE.1)TYFE"2REBLK  #",SB,"  error :",IER 
GO  TO  100 

Store  the  distances  computed  for  iteration  II  (or 
temporarily  store  the  energies  if  EM.) 

500  IF(SaJEHGY.LT.1.0)SENERGY»1.0 
IF(TQEBGY.LT.1.0)TENERGY*1.0 
IF  (SL1B4EBGY.LT.1 .0)  SL1ENERGY-1 .0 
IF  (TL1ENEPGY.LT.1 .0)TO1ENERGY=1 .0 
S2N0HMQKT  (SENERGY) 

T2NDREHSQKT  (TENERGY) 

S1N0FE*=SL1  ENERGY 
T1EDFM-TLJ.  ENERGY 
AREA-FLOAT  (LENGTH )  *FLQAT  (WIDTH) 

IF(N.BQ.1)G0  TO  40 

Confute  distance  factors 

NDISTL1  (II)  -NL1DIST 

NXYMIN-TO1  ENERGY/  (SQFT  (AREA*TENERGY) ) 

NL2E1AX-SQFT (2 .0* (1 .0-NXXMIN) ) 

OORPEAK(II)  KDRREL/  (S2NDRM*T2N0RM) 

CFACTOR(II)  »INT  (100 .0*  (GCH®EAK  (II) -NXYMIN) /(l  .0-NXYMIN) +0  .5) 
IF (CFACDOR(II)  .LE.0)CFACrOR(II)=O 
NL2«SQFT(2.0*(1.0-OORPEAK(II) ) ) 
L2FACTCK(II)=INT(100.0*(1.0-rt2/NL2MAX)-K).5) 

IF(L2FACP0R(II)  .LE.0)L2FACTOR(II)  =0 

LIFACTOR(II)  »IOT  (100 .0*  (1 .0-EDISTLl  ( II )  /NL1MAX)  +0  .5) 

IF(LlEftCTQR(II)  .LE.O)L1FACFQR(II)=0 

Reset  the  template  window  parameters 

TTO^rrQP 

tleet-ileft 

wioth-iwu/ih 

LENGTH- ILENGIH 

500  CONTINUE 

TYH3*<7><7>*******************************" 


********  write  RESULTS  TO  t.tmrpp twprr  ************************* 


User  input  of  comment 

DO  625  1-1,200 

C0MIENT(I)»0 

625  CONTINUE 

ACCEPT  "Enter  a  1  to  add  comment  to  the  output:  ",IOPT 
IF(IOFT.NE.1)QO  TO  650 

ACCEFT"Enter  #  of  comment  lines  (max-4):  ",NUMCDM 
IF  (NUMQOM.LT.  1)NUNQ0M-1 
NUMQ0M-MIN  (NUMQ0M,4) 

DO  640  I-lfNUMQOM 


non  _  ^  r>nn  non 


TYPE  "Enter  comment  line  #”,I,"  to  be  printed  with 
$  results  between  the  arrows:” 

TYPE”  |  <11><11><11><11><11><11><11><11><11>  I  ” 
TYPE”V<11><11><11><11><11><11><11><11><11>V* 

READ  (11, 9999)  OOMMSNTf  (50*(I-1)+1) ) 

640  CONTINUE 

Write  output  header 

650  WRITE(12,9005) 

WRITE  (12, 9000) 

WRITE (12,8000) 

IF(RED(JCED)WRITE  (12,7500) 

WRITE  (12,7000)  LENGTH, TTCS5; ,  I’Ll  ENERGY 
WRITE (12,6000)  INFILEl(l) , WIDTH, TLS;  ,TENERGY 
IP  (RHXJCED)  WRITE  (12,5500)  INFILE2  (1) 

IF  ( .N0T.RHXJGED)  WRITE  (12 , 5000)  INFILE2  (1) 

WRITE (12, 4600) 

WRITE  (12 ,4500) 

WRITE  (12,4200) 

WRITE  (12,4100) 

Write  distance  factors 

11=10 

DO  710  JJ=1,NUMPEAKS 
DO  700  KK=1,9 
11=11+1 
A=CPACPOR(II) 

B»L2PACTQR(II) 

OLIFACPOR(II) 

SCOREMMTf  (A*B*C)**(1  *0/3*0) +0.5) 

IF(KK.1E.1)G0  TO  698 

WRITE (12 ,4000)  JJ, PROW (JJ)  ,POXUW(JJ)  ,R0WCENT(II)  ,aOLCEMT(II) , 
$CTQP(II)  ,SLEET(II)  ,CFACTQR(II)  ,L2FACPQR(II)  ,L1FACP0R(II) , 

$SG0RE 

4000  FQRHAT(14X,I2,”:”,I4,”,”,I3,5X,I3,”,”,I3,6X,I3,4X,I3, 

$7X,  13 ,6X,  13 ,5X,  13 , 5X,  13  ,T132 ,  ”  ") 

GO  TO  700 

698  WRITE(12 ,4050) ROWCENT (II)  ,O0LCENT(II)  ,C!EOP(II)  ,SLEET(II) , 

$CFACPOR(II)  ,L2FACTOR(II)  ,L1FACT0R(II) , SCORE 

4050  FQBMAT(30X, 13 , , 13 ,6X, 13 ,4X, 13 , 

$7X,  13 ,6X,  13 ,5X,  13 , 5X,  13  ,T132 ,  ”  ”) 

700  CDNTINIE 

WRITE (12, 40 51) 

4051  F0RMAT(”  ") 

710  CONTINUE 

Write  comments  to  linepr inter 

IF(IOPT.EE.1)GO  TO  800 
DO  790  I=1,NUMQ0M 

WRITE  (12 ,9500)  OOMMHT  ( (50*  (1-1)  +1) ) 

790  CONTINUE 
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800  HRITE(12,9000) 

Format  statements 

9999  FORMAT  (S100) 

$500  F0RMAT(16X,  "COMMENT:  "rSl00) 

9005  FORMAT (///////////) 

9000  P0RMAT(/,11X,81("*"),T132,"  ") 

8000  FORMAT  (////37X,  "RECOGNITION  RESULTS  <10>"f///) 

7500  FORMAT(15X,"  ****RHXXSD****") 

7000  F0RMAT(15X,"  TEMHATE  WINDOW :",7X,"LENG1H=",  13 ,"  ROWS",  12X,  "TOP 
$I3,10X,T132,"  ") 

6000  FORMAT (16X,"  (",S13,")  ",8X,"WiniH=M3,"  COLUMNS", 9Xr"LEFTOOL=", 
$I3,10X,T132,"  ") 

5500  FORMAT(//f30X, "*REEUCED*  SCENE  FILE - >  ",S13,T132,"  "r//) 

5000  F0RMAT(//,35X,  "SCENE  FILE - >  ",SI3,T132,"  *,//) 

4600  FORMAT(14X,"  CORRELATION  WINDOW") 

4500  FORMAT  (14X,"  PEAK  CENTER  TOP  LEFT 

$  CORRELATE  L2  LI") 

4200  FORMAT  (14X, "  (ROW,CDLt*N)  (RCW,  COLUMN)  ROW  CDLIMJ 
$  FACTOR  FACTOR  FACTOR  SCORE”) 

4100  FORMAT  (14X, " -  -  -  - 

§ - ") 

TYFE"<7><15><11>****  CHECK  LINEPRINTER  FOR  RESULTS  *****<15>" 

C 

C******  Present  Option  Menu  ********************************* 

C 

GO  TO  2010 

2002  TXPE"<15>"r "Input  error. <7>  Try  again." 

2010  TYFE"<15>", ******************************************* 

§**********************■ 

TYPE"<15>"f "What  next?<15>"f"Here  are  the  options:" 

TYPE"<15>", "<11>1  -  Try  another  set  of  windows" 

TYFE"<15>" , "<11>2  -  Start  over  with  new  input  pictures” 
tyh:"<i5>","<ii>3  -  srap<i5>" 

ACCEFP"<ll>Enter  option  — >  ",IOFT 
IF(IOFr.Ur.l.CR.IOFT.<7r.3)GO  TO  2002 
TYPE"<15>" 

IF(IOPT.E)Q.1)GO  TO  3 
CALL  RESET 

IF(IOPT.EQ.2)QO  TO  99 
STOP 
BID 
C 

C*********  Program  DISTANCE  ******************************************** 
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SUBROUTINE  EUCLID  (TSTART,  SSTART,  N,$) 
£********************************************************************* 

(Called  by  DISTANCE )  by  Lt  Jim  Cromer 

If  ^1  — >  calculate  LI  and  L2  energies  of 
template  and  soene  windows 
Else  — >  calculate  the  normalized  LI  and  cross- 
correlation  measures  between  the  windows 

TSTART,  S START  are  the  row  position  within  the  packed 
video  block  of  the  first  row  of  the  window.  They 
<u:e  automatically  incremented  after  the  first  call  to 
EUCLID  (TSTART,  S START  between  0-3  inclusive) . 

********************************************************************** 
REAL  NL1DIST ,  LlDIST ,  NSUMSQ 

INTEGER  SCENE  (1024) , TEMP  (1024)  fSLSfTLS, SSTART, TSTART, WIDffl 
COMMON  /LISTl/  SOME, TEMP, SLS,TLS, LlDIST, NLIDIST 
OOMMON  /LIST2/  LENGTH, WII7IB,]ECW00UNr,SUMSQ,NSCMSQ 
OOMMDN  /LIST3/  SUJNEIGY,1L1ENERGY,SENERGY,TE!NERGY 
OOMMDN  /LIST4/  S1NORM,S2NOFM,T1NDRM,T2NDFM,OORREL 

Set  do  loop  parameters 

JMIN=MAX  (SSTART, TSTART) 

DO  2  J=JMIN,3 
K=TSTART*256+TLS 
M=SSTART*256+SLS 
KMAX=K4WIETH-1 
IF(N.BQ.l)GO  TO  3 

Calculate  the  distances 

CO  1  L=K,KMAX 

RSCENE>FLQAT  (SCENE  (M) ) 

RTEMP^ELQAT  (TEMP  (L) ) 

MilDIST-ABS  ( (RSCENE/S1N0PM)  -  (FTEMP/TINOPM) )  +JL1DIST 
OORREL"  (RSCENE*RTEMP)  +CORREL 
M**Mfl 
CONTINUE 
GO  TO  5 

Calculate  the  energies 

3  DO  4  L*K,KMAX 

RSCB&>FLQAT  (SCENE  (M) ) 

BTEM^FLQAT  (TQ(P  (L) ) 

TL1ENERGY-TL1ENERGY+PTEMP 
SL1ENERGY-SL1ENERGY+RSCENE 
SENERGY-  (RSCBE**2)  +SENERGY 
TENERGY-  (FTEMP**2)  +TB1ERGY 
M-Mfl 
OONTINIE 
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Test  for  the  end  of  the  window 


5  IRQWC)CXJNP»IROWaXJNPtl 

IF  (IB0W0CUNT.GE.UNG1H)  RETORN  4 

Increment  block  row  oounters 

TSTAKT=*TSTART+1 
2  SSPARP«SOTAR3M-1 

IF (SSTAKT.BQ.4) SSTART=0 
IF  (TSTART.  EQ.4)  TSTART*0 
RETURN 
END 

C********  subroutine  EUCLID  ******************************************* 
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SUPPORT  SUBROUTINES 


This  appendix  contains  the  following  programs 

1.  IOP 

2 .  TIMER 

3.  UNPACK 

4.  XRDBLK 


148 


non  non  onononoooooooooonooooooooooooo 


SUBROUTINE  K)F(N,MAIN,F1,F2,F3,MS,S1,S2,S3) 

C**************************************************************** 

Written  by  Lt.  Simmons  10  Sep  1981 

Version  2 

This  FORTRAN  5  subroutine  will  read  from  the  file 
COM. CM  (FOQM.CM  in  the  foreground)  the  program  name, 
any  global  switches,  and  up  to  three  local  file 
names  and  corresponding  local  switches. 

Calling  arguments: 

N  is  the  nunber  of  local  files  and  switches  to  be 
read  fran  (F)OOM.CM.  N  must  be  1,  2,  or  3. 

MAIN  is  an  ASCII  array  for  the  main  program  file  name. 

FI,  F2,  and  F3  are  the  three  ASCII  arrays  to  return 
the  local  file  names. 

MS  is  a  two-word  integer  array  that  holds  any  global 
switches. 

SI,  S2,  and  S3  are  two-word  integer  arrays  that 
hold  the  local  switches  corresponding  to  FI  through 
F3  respectively. 

****************************************************************** 
Dimension  the  arrays. 

DIMENSION  MAIN (7)  ,MS(2) 

INTEGER  FI (7)  ,F2(7)  ,F3(7)  ,S1(2)  ,S2(2)  ,S3(2) 

Check  the  bounds  on  N. 

IF(N.I/r.l.0R.N.OT.3)SPCP  "N  out  of  bounds  in  IDF." 

Process  the  data  in  (F)GOM.CM 

CALL  GCUND(I)  ;Find  out  which  ground  program  is  in 
IF(I.BQ.O)OPEN  0, "COM. CM"  /Open  ch.  0  to  (DM. CM 

IF(I.EQ.l)OPEN  0,"FG0M.CM"  ;Open  ch.  0  to  FOOM.OI 

CALL  Q0MARG(0,MAIN,MS,IER)  ;Read  fran  (F)QQM.CM 

IF (IER.NE.l)TYPE"  COMARG  error:", IER 
WRITE (10,1) MAIN (1)  ;Type  program  name 

1  FORMAT ( '  Program  ' ,S13,' running.') 

CALL  COMARG (0,F1, Si ,JER)  /Read  fran  (F)COM.CM 

IF(JER.tE.l)TYPE"  COMARG  error  (F1):",JER 
IF(N.BQ.1)G0  TO  2  /Test  N 

CALL  COMARG (0,F2,S2,KER)  /Read  from  (F)OOM.CM 

IF(KER.ie.l)TYFE”  QOMARG  error  (F2)  :",KER 
IF(N.EQ.2)G0  TO  2  /Test  N 


CALL  CDMARG(0,F3,S3rIZR)  ;Read  from  (F) COM. CM 

IF(LER.fE.l)TYPE"  QOMARG  error  (F3):",LER 
2  CHOSE  0 
RETURN 
END 
C * 

C************  subroutine  IDF  *************************************** 


oooooooorjrsooooo 


SUBROUTINE  TIMER(I) 

C****************************************************************** 

Subroutine  TIMER  Written  by  Lt.  Jim  Cromer 

Fortran  5 

This  subroutine  is  used  to  time  the  real-time  execution 
time  of  the  calling  program.  If  the  parameter  passed,  I, 
is  equal  to  0,  the  timer  is  unconditionally  started. 

If  I  is  not  equal  to  0,  the  timer  is  unconditionally 
stopped,  and  the  total  run  time  is  typed  on  the  console 
CRT. 

Execution  Line  Format 

CALL  TIMER(I)  ;IF(I.EQ.O) ,  start  timing 
jIF(I.ie.O) ,  stop  timing 

C******************************************************************* 

COMMON  /ITIME/  IH1,IM1,IS1 
IF (I. NE. 0)00  TO  100 

CALL  FGTIie(IHl,IKL,ISl)  ;get  starting  time 

WRITE  (10  ,1000)  IH1 ,  IM1 ,  ISl 

1000  FORMAT (//"  START  TIME  - >",I4,":",I3,":",I3) 

RETURN 

100  CALL  FGTIME(IH2,IM2,IS2)  ;get  stopping  time 

WRITE  (10, 2000)  IH2,IM2,IS2 

2000  FORMAT (//"  STOP  TIME  - >",I4,":",I3,"s",I3) 

ITOTAL=3600* (IH2-IH1) +60* (IM2-IM1) +IS2-IS1 
B0URS=INT (ITOTAL/3600) 

TECN=(ITO!EAL-3600*H0URS)  } intermediate  variable 

MINS=INT(TRCN/60) 

ISECS=MQD(TRON,60) 

WRITE  (10 ,3000)  HOURS, MINS,  IS  ECS 

3000  FORMAT  (//"  TOTAL  TIME  - >",I4,":,,,I3,":",I3) 

RETURN 

END 

C 

C**************  subroutine  TIMER  ********************************* 


on  n  nnnnnnnnnnnnnnnnnnnnnnnnnnnnnnn 


C**************************************************************** 

Unpacking  (packing)  routines 
Written  by  Lt.  Simmons  Version  2 

Documented  by  Lt.  Crane r 

These  subroutines  unpack  (repack)  four  4-bit  integers  from 
16-bit  integer  word.  The  pixels  in  a  video  file  have  to 
be  unpacked  if  each  pixel  is  to  be  operated  on  separately. 

Packed  video  (4  pixels/1  word) : 


WORD(Nfl)  I  PIXEL  1 1 PIXEL  2 1 PIXEL  3 1 PIXEL  4| 


Unpacked  video  (4  pixels/4  words) : 


WORD(X  )  | 

<-  unused  -> 

1  PIXEL  1| 

WOHD(X+l)  | 

<-  unused  -> 

|  PIXEL  2| 

WORD(X+2)  | 

<-  unused  -> 

|  PIXEL  3| 

WORD(X+3)  1 

<-  unused  -> 

| PIXEL  4| 

where  N=X  mod  4 

****************************************************************** 
SUBROUTINE  UNPACK  (N,PIXW3RD,  PIXELS) 

INTEGER  PIXWORD(N)  ,PIXELS(4,N) 

DO  1  1=1, N 
DO  1  J=l,4 

PIXELS  ((5-J)  ,I)=15  .AND.PIXWORD(I) 

1  PIXWORD(I)=ISHFT(PD«ORD(I)f-4) 

RETURN 
END 

SUBROUTINE  REPACK (N, PIXELS,  PXWD) 

INTEGER  PIXELS ( 4, N)  ,FXWD(N) 

DO  1  J*1,N 
PXWD(J)-0 
DO  1  1*1,4 

PXWD ( J)  =ISHFT  (PXWD  ( J)  ,4) 

1  PXWD  ( J)  “PIXELS  ( I ,  J)  +FXWD  ( J) 

RETURN 
END 

************  packing  Subroutines  ********************************* 


;Four  pixels  per  word 
;'N'  allows  higher-order 
;arrays  to  be  passed. 

; Pick  off  right  pixel 
;Shift  word  4  bits  right 
;to  pick  off  next  pixel. 


SUBROUTINE  XRDBLK (CH,J, FILE, I, IER) 

C********************************************************************** 

C 

C  by  Lt.  Jim  Cromer 

C  Subroutine  XRDBLK  performs  a  RDBLK  to  the  designated 

C.  channel,  reads  a  packed  video  file  block,  and 
C  returns  an  unpacked  array. 

C 

C********************************************************************** 

INTEGER  CH, FILE (1024)  ,VIDEO(256) 

K=256*I 

IF(J.GE.O  .AND. J.LE.63)  GO  TO  1 

TYPE  "ERROR:  <7 > BLOCK  POINTER  OUT  OF  BOUNDS  IN  XRDBLK" 

TYRE"  J=",J 

STOP 

1  IF  (I. EQ. l)QO  TO  2 

TYPE "ERROR  IN  <7>XRDBLK" 

TYPE"  #  Blocks  to  be  read  =",I 
STOP 

2  CALL  RDBLK (CB,J, VIDEO, I, IER) 

DO  3  L=1,K 

DO  3  M=l,4 

ICXXJNT=5-M+-  (L-l)  *4 

FILE  (ICOUNT)  =15.  AND.  VIDEO  (L) 

VIDEO  (L)  =ISHFT  (VIDEO  (L)  ,-4) 

3  OONTINIE 
RETURN 
END 

C 

C*********  Subroutine  XRDBLK  *********************************** 


APPENDIX  H:  PRELIMINARY  RESULTS 
Summary  of  Tank  SCORES 


Template  and  Tank  windows  both  globally  normalized 


Template 

Scene 

High 

Low 

Average 

PTEMPH3 

PTANKH3 

35 

0 

19 

RPTEMPH3 

RPTANKH3 

36 

0 

5 

PTEMPD4 

PTANKB2 

41 

0 

25 

PTEMPD4 

PTANKC3 

49 

0 

15 

PTEMPD4 

PTANKE2 

56 

0 

31 

PTEMPD4 

PTANKG4 

0 

_0 

_ 0 

Average 

36 

0 

16 

Template  and 

Tank  windows 

both  grid 

normalized  (9x5  grid  used) 

Template 

Tank 

High 

Low 

Average 

N0RMD4 

N0RMB2 

65 

40 

52 

f  » 

N0RMC3 

69 

43 

58 

9  9 

N0RMD2 

64 

43 

54 

»  9 

N0RME2 

67 

45 

57 

9  9 

N0RMG4 

42 

16 

29 

f  1 

N0RMH3 

79 

24 

43 

Average 


64 


35 


49 


Summary  of  Scene  SCORES 

Template  and  Scene  windows  both  globally  normalized 


Template 

Scene 

High 

Low 

Average 

PTEMPH3 

PSCENEI3 

16 

14 

15 

PTEMPH3 

PSCENE04 

24 

18 

21 

RPTEMPH3 

RPSCENE04 

20 

10 

15 

RPTEMPH3 

RPSCENELl 

19 

14 

16 

Average 

19 

14 

17 

Template  and  Scene  windows  both  grid  normalized 

(9x5) 

Template 

Scene 

High 

Low 

Average 

N0RMD4 

WHITE 

0 

0 

0 

N0RMD4 


N0RM7 


15 


2 


8 


VITA 


James  H.  Cromer  was  born  on  10  April  1959  in  Cleveland,  Ohio. 
He  graduated  from  Kenston  High  School,  Chagrin  Palls,  Ohio, 
in  1977.  He  attended  Grove  City  College,  Grove  City,  Penn¬ 
sylvania,  on  a  4-year  AFROTC  scholarship.  He  received  the 
Bachelor  of  Science  degree  in  Electrical  Engineering,  and  was 
commissioned  a  Second  Lieutenant  in  the  United  states  Air 
Force,  in  May  1981.  He  entered  the  School  of  Engineering, 
Air  Force  Institute  of  Technology,  in  June  1981,  and  was 
chosen  to  receive  the  1981-1982  IEEE  Outstanding  Student 
Award  in  May  1982.  He  is  a  member  of  Tau  Beta  Pi,  Sigma  Pi 
Sigma,  and  Eta  Kappa  Nu. 

Permanent  Address:  204  E.  Mill  St. 

Port  Allegany,  PA  16743 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Dim  gnl.r.Q 

I  REPORT  DOCUMENTATION  PAGE 


AFIT/GE/ZE/82D-26 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 
iciwiSri  catalog  number 


4.  TITLE  (mnd  Submit) 

SCENE  ANALYSIS:  NON-LINEAR  SPATIAL  FILTERING 
FOR  AUTOMATIC  TARGET  DETECTION 


S.  TYPE  OF  REPORT  •  PERIOD  COVERED 

MS  Thesis 

4.  PERFORMING  ORO.  REPORT  NUMBER 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  A  WORK  UNIT  NUMBERS 


7.  AUTHORfaJ  0.  CONTRACT  OR  GRANT  NUMBERS) 

Janes  H.  Croner 
2nd  Lt  USAF 

»•  PERFORMING  ORGANIZATION  NAME  AND  AOORESS  10.  PROGRAM  ELEMENT.  PROJECT,  7 

AREA  A  WORK  UNIT  NUMBERS 

Air  Force  Institute  of  Technology  (AFIT/EN) 

Wrlght-Patterson  AFB,  Ohio  45433 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS  12.  REPORT  DATE 

December  1982 _ 

is.  NU^ER  OF  PAGES 

14.  MONlfbhlNG  AGENCY  NAME  A  ADDRESS^  SmNn oi  from  Controlling  office)  IS.  SECURITY  CLASS,  (ei  HU.  report) 

Unclassified 


ISa.  OECLASSIFICATION/DOWNGRAOING 
SCHEDULE 


14.  DISTRIBUTION  STATEMENT  (el  Ihlr  Keporl) 


Approved  for  public  release; 

distribution  unllnltsd 

"*7.  DISTRIBUTION  STATEMENT  (et  Hie  aAalrael  entered  In  Bloc*  20.  II  dltlerent  tram  Report) 


I  14.  SUPPLEMENTARY  NOTES 


An^rovod  loV,  public  rmleasn: 
WOLAVEH 


IAW  APR  190-!?. 


^ean  lot  Reseat  th  and  Professional  Devnlt] 
Per O*  Institute  oi  Uchoology  (ATGJ 

SrWwM-fW cnee  AfS-QH  4WM _ 

IS.  KEY  WORDS  (Continue  on  imni  elde  It  necmerery  end  Identity  b y  block  number) 


JAN  1S82 


Cross-correlation 
Fora  Analysis 
Inage  Processing 
Energy  Nozmallsatlon 


Pattern  Recognition 
Scene  Analysis 
Target  Classification 
Target  Detection 


Fortran 


20.  ABSTRACT  fCenllnue  an  reeeree  tide  If  naeaaaarr  aN  Identity  by  block  number) 

This  work  focuses  on  a  method  for  two-dimensional  pattern  recognition. 
A*  method  includes  a  global  search  scheme  for  candidate  windows  of  interest, 
based  on  Fourier  domain  cross-correlation.  A  method  to  normalise  the  Input 
scene  by  local  rectangular  regions*  In  an  attempt  to  efficiently  approximate 
search  window  normalisation,  is  presented.  Also  developed  la  a  candidate 
window  (potential  target)  similarity  measure,  based  on  the  normalised  LI  and 
Isclldean  distances,  which  Is  Independent  of  the  template  DC  value  and  Its 
energy.  Observations  on  the  performance  of  the  algorithm  applied  to  visual 

DO  ijSTt,  1473  COITION  OP  I  NOV  44 II  OBSOLETE  -  ... 


Unclassified _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  fWheo  note  Fn 


StCUWITV  CLASSIFICATION  OF  THIS  PAOEflWmi  Dmta  Bnfnd) _ 

spectrum  photographs  of  tanks  In  a  realistic  environment  are  Included.  Also 
Included  is  the  software  needed  to  implement  the  algorithm  on  Data  General 
Eclipse  S/250  minicomputer. 


Unclassified 


